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ABSTRACT 


A nonlinear model characterizing two terminal 
microwave active circuits under dynamic and steady state 
conditions is developed and represented by two non- 
linear differential equations. The model is applicable 
to oscillators and amplifiers, and considers the 
admittance dependence on both RF voltage and frequency 
up to the sixth order. The model response to harmoni- 
cally time dependent signals with slowly varying ampli- 
tude and frequency is discussed. Either discrete or 
analytically characterized data can be employed in the 
model. 

The model is applied to Impatt diode oscillators 
and the computed transient and steady state solutions 
are compared with experimental data. A technique for 
measuring large signal admittance parameters is developed 
and is used to obtain the model parameters of an Impatt 
diode and its embedding circuitry. 

In order to measure transient behaviour of Impatt 
diode oscillators, a new oscillator designed for pulsed 
operation and a new frequency discriminator was deve- 
loped. Some of the results measured on a circuit compared 
favourably with those obtained from the model using the 


measured parameters. 
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The transient response is also evaluated for two 
hypothetical oscillators in which a semiconductcr wafer, 
as characterized in the literature, is embedded. The 
results indicate the sensitivity of the transients to 
the embedding circuitry. 

As a practical application of the bias modulated 
oscillator, a digital modulation scheme which operates 


at rates up to 40 Mb/s has been suggested. 
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Illustrating the basic principle of 
the operation of the microwave 
instantaneous frequency measuring 
system, 


Single hybrid tee frequency discrimin- 
ator: 


Polar display of the discriminator 
response (P. = 10 mW, 2% = 1.78 cm). 


Input VSWR of the single hybrid tee 
discriminator. 


Experimental arrangement for the 
calibration of instantaneous power 
and frequency measuring system, 


Detectors dc voltage sensitivities, 


Constant power and frequency response 
of the instantaneous measuring system, 


Practical RF modulating system to check 
the accuracy of the instantaneous power 
and frequency measuring system, 


The two RF envelopes recorded by the 
X- and Y-detectors, 


Enlarged section of the constant power 
and frequency curves of Fig. 6.7 
showing the results extracted from 

the traces in Fig. 6.9, 


Overall iterative scheme for the cal- 
culations of steady state voltage 
amplitude, A, and frequency wt+x, 


Overall iterative scheme for the cal- 
culation of the dynamic solution by 
using forward iterative scheme, 


Overall iterative scheme for the cal- 
culation of the dynamic solution 
using backward iterative method, 


ow 


AS he 


AES fee 


fo: 7 


156 


SS ie 


159 


ca | 


a 


>» 


~~ 


a wesike ae 
we 1s wausst you spe 


~pimbiaeth Yousvpess oer 


| 


2H teak “aa kb oft’ 20 wank. zeke 


¥? ‘mo BY 2 F aia = vin : Os +, > baie: 


“Gen ‘hivxe Vii af aD | te 3 “it ‘36 gwen. fiat 


Paey 
at Hal . salle _psogpnimt 102L8 | 
ve ae : ait yr a = x a oy 
gd 408 J neoorrtatre B atnsiatat wt seat MW 


a 
at 


so 4 


Ry Ee " f Liles ‘ er? , 
move paitvaser peacupeth. bee vite 


‘ PT a i yy , 
oe ' eo oe ee eur . — » 
Bes as in I ~ Lo & ae us ud a ¥ ~ tee a4 ey Glad aoe 4 <2 
» 
- itn ms) i ap os + a) = 
i shania ~ ee 4 date any ‘ne el) Be orien 
- 7% c POAT £3 bits rewotd SnBIen0.. 
= a 
, ¢ ; - we . a * : 
“news. priveeae exostrarnstens oie 2. 
rer a or oe iy 
q : i ‘ Po Ly , o Pereyra ae. r% > as . 
n " J er . ze tee: *F Nie) ars Bria: Se rs P 
F 7 é new wh ee NS . 


hits 2.20 yoo wb)k mee 7 
Hoye, Dabtvanem yousgpea? bas 


‘ ; - 


eis yt Dabtoney aeqoleved i om eat 


r ‘* ne | c er 
bine i i i - 
2 _ Nar a ‘ 4 sy ot a he & 7 ¢ Liste ., A 
Se 7 2 1 
mo ee ae au Be we . Bi t Pa if da 
i ; Ur ihis a5 \ B29 ToD teal ows us a eS 


-. §S OE Se @evayo vortkerees* bre 


r > es 
™ x Py 4 , se 7 ve sae 
f Af Ae cf SR i cot 4 as = os at eet a BS ad 
* : 7 ry 
’ » ne x4 ' 
7 ead * res > taal heme i .§ oar 


: e | : 
~Las. chat: 2m, om pio avis erat . Liexeva- ee : 
bee boy eet wheter to enoetedaie: % 
x @ Yor apasS ee a pghts +x res: et Goo 

~lan ody tod Gpegeine ew bie stedh Hexsvo ia 


baode Weer ost.) sG3-2o onda ative. oe 
ie - ebtatesi brewact- gated. | 


= bso bi!” dsbortiahe gt Le cod. Adeseue: . 
1 af Pe; iy 2 wlnsing) 963-26, hatdadue 
| feat on. etwas Syeenes . 


x 


Figure 


724°, (a) 


(b) 


Impatt dicde admittance as a function of 
RF voltage amplitude and frequency ata 
current density of 200 A/cm. 


Pacotve CitCurt adiiicance, vip asia 
FUNCT OM Ol gL Leguency LO! ClYrcui. A 
(Table 7.1)" 


Passive circuit admittance as a function 
OULperecucncy wfOr Cilrcudtee lable 7 /1).. 


Cross-sectional view of a 7 mm coaxial- 
type resonator embedding the Impatt 
diode. Rz,Ro,R3,XhQ and Alo are,the 
resonator parameters (Table 7.1). . 


Equivalent circuit.of an Impatt diode 
OsGLilaton,. 


Buiiaup 4. Of power in a realistic oscil= 
lator for two, sets of the 7 mm coaxial 
TEecGnaton DalaMmecters. 


Buildup of frequency in a realistic 
oscillator for two.sets of the 7 mm 
coaxial resonator parameters. 


Buildup of power (a) and frequency (b) 
in the Impatt diode oscillator using 
qode 1 2 Ch TV LiL). 


Device line and impedance locus super- 
imposed on a part of Impatt diode 
admittance characteristics. 


Device line and impedance locus super- 
imposed on the measured Impatt diode 
agmittance Characteristics, ydiode ui2 
(Gy = 5 B66 (GHZ). 


Arrangement for the measurement of 
transients in Impatt diode oscillators. 


Cross-sectional view of the Impatt 
diode oscillator used in the measure- 
ment of the transients. 


Cross-sectional view of the 7 mm coaxial 
cavity with end mounted Impatt diode. 
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Bias circuit as seen from pulse 
generator terminals. 


Modulating pulse (500 mV/cm), diode 
current (100 mV/cm) and two RF 
envelopes recorded by the X- and Y- 
detectors (20 mV/cm) for diode #2. 


A reproduced copy of Fig. 8.5 with 
the time scale adjusted for a 10 ns 


delay caused by the measurement system. 


System calibration for instantaneous 
frequency and at different instance 

of time on the RF envelopes of Fig. 

BE 


Measured and calculated power 
waveforms of the Impatt diode 
Oscillator. 


Measured and calculated time depen- 
dence of the frequency of the Impatt 
diode oscillator. 


Diode current (100 mV/cm) and two RF 
envelopes (20 mV/cm) recorded by the 
X- and Y-detectors for a modulating 
pulse with duty cycle of 0.5, 


Diode current (100 mV/cm) and two RF 
envelopes (20 mV/cm) recorded by the 
X- and Y-detectors for a modulating 
pulse with duty cycle of 0.2. 


Frequency variation, for pulses of 30 
ns ON time, as. a Eunction of the 

duty cycle between two instances for 
which the power changes between its 
maximum value and 87% of this value. 


Decay time of the RF envelopes as a 
function of the duty cycle due to a 
36 ns ON pulse. 


Instantaneous frequency at maximum 


power level vs. pulse periods (50% 
ON time pulses). 
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Figure Page 


625 Modulation scheme for utilizing bias 208 
modulated Impatt diode oscillator 


for the transmission of data at nearly 
40 Mbits/s. 
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LIST OF SYMBOLS 


A= Ch. Tit: RF current amplitude flowing through the 
active device. 
Ch. JIVis>- diode area. 
Otherwise = A(t): time dependent applied RF voltage 
amplitude. 


A,B,C and D = the fundamental parameters of the insertion 


network. 

a,b,c and 1 = normalized parameters of the insertion 
network. 

AN,BN,CN and 1 = elements of the error network matrix. 

An = peak amplitude of the driving RF voltage for CW 
excitation. 

Aah’ Soft = constants depending upon the embedding circui- 
try and active device. 

anh’ Doh = coefficients dependent on the active device, 
embedding circuitry and the bias voltage KpaC 

anh, BLES ji eS m-th coefficient of the g-harmonic 
component with an amplitude AG 

Oe simplified notation for oon 

ac = dA/dt. 

Ayres nA = characteristic parameters of the passive 
erro ulti. 

a, = normalized incident wave at the i-th port 
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Direserbo = characteristic parameters of the passive 


circa t. 


oO” 
Il 


normalized emerging wave from the i-th port 


lle 


2 
Go eats ea: Joe 
bi = emerging wave from the i-th port to the j-th port. 
d of as , 

RL os heise ie Daclmcocttvelenmeror the .g—th barmonic 
component with an amplitude ne of the active 
device in the KL subrectangle [A, Ay! x TW) Xx Wray 

zbe = the p-th coefficient of the passive circuit in 

p,90,0 


the L-interval. 


d 


hie, 2hyg Scie ola yer g except, Lome hoe.(2...63) to (2.68). = 


the total oscillator b-coefficients. 
B(AZ, we) = active device susceptance determined at the 


points ae and a 


L 

BAC, wi) = active device susceptance determined at the 
oints AG and ae 
Pp M Wyre 


Fen or bicubic spline interpolate for the active 
device susceptance determined at the points ne 


ae 
and w in the MN subrectangle [Ay Ayyy!) x Loy ey yale 


Ria ep = active device susceptance at the points aC 
and fae 

BL (A7,w*) = B(A’,w-). 

Bia (A ew) = susceptance of the complementary device 


admittance at the points Be and fc in the MN 
subrectangle. 


BL = diode susceptance. 
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B, (4) = circuit susceptance. 


yBo (w”) = spline interpolate for Bo in the N-th interval. 

B = total circuit susceptance of the oscillator. 

Br =», ‘octalwcircuit complementary susceptance of the 
oscillator. 

Ban = device electronic susceptance at the n-th harmonic. 

Be = device electronic susceptance. 

By = load susceptance. 


Bot C,= Ch. Vit: nonlinear coefficients of the equations 
representing the steady state of the oscillator. 
Cy = Chivtilr*passivev-circurt capasicance 
Ch.III: equivalent capacitance of the highly 
SaMmps i rreda  oscritator Circuit: 


CniVill= “bypass” capacitors. 


Ce = circuit capacitance at the n-th harmonic. 

Co = equivalent load capacitance. 

C, = equivalent capacitance of the cavity. 

eo = equivalent capacitance of the device. 

ce = equivalent passive circuit capacitance. 

C(V) = bias voltage dependent capacitance of the Impatt 


diode wafer (V<V)). 
evel =e diode wafer capacitance at bias voltage V5. 


= diode wafer capacitance at bias voltage 


Q 
II 


the average of the CONST EE and C541 at bias 


voitages V. and eee) , respectively. 
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package capacitance. 


heat sink capacitance to ground, 


a ay 
Km KLM 2mt1, 2° KLM 2m, 2 APF x omer 


coefficients in the KL subrectangle [Ay , A, 47! 


= acti ver aevice 


[ w which are functions of A and the 


ee tre 


device d-coefficients. 


~ we a fon = : . 
MNC~2m ¢ mnewom+1, 2! wnC*2m, 2 and MNCYOm+1 active device 


coefficients in the MN subrectangle. 


CW = constant frequency wave. 

dona (Prh)-th coefficient of the general model of the 
active device derived from the ct pK coefficients. 

ao ate Mate coefficient of the simplified model of the 
active device derived from the rsp beara coefficients. 

oe p=th Goefirecient “of ther cixrcurtitderived from the 
Behe. coefficient. 


dBm = dB with .reference to mW. 
D = Cha 510 = a0 


Gh. /V* s*diodestyunction fdvdmeterm Ceni):. 


E = eléctric fieid w(V/em),.. 

E,,, (t) = voltage amplitude of the measured signal. 

Eout = Measured output voltage. 

Eq Ea = de output voltages of diodes #1 and #2. 
1 2 

£N.= ‘Sfrequency in GHz. 


f(t) = driwine@vicorce: 
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Py reeesFy = functional values of either circuit conduc- 
tance or B.(w)/w. 

Fy (u") = spline interpolate for G, ox Bw) /w in the 
N-th interval. 

g = the aj-th,.coefficient,of ithe bicubic spline 
fitting for the active device conductance. 

g = hole and electron generation recombination rates. 


G = Choditis, wassive circuit conductance pp. 6-7. 


total. circuit conductance of the oscillator 
eo. 
Ch.III: equivalent circuit conductance of the 
highly simplified oscillator circuit. 
ate = 
Glw (t)] =conductance operator for non-CW operation. 
Bi 
-G(w ) =the spline interpolate in the L-th interval for 
i, & 
the passive circuit conductance. 
2 OS ae fe 5 ; 1.7 pen 
G_ (wp) = the m-th coefficient of the spline in the L-th 


interval for the passive circuit conductance 


(i =. dig 2 pdy 4) v0 


Bo 2 ; ; : , 
Kp (AX rw )=the spline interpolate in the KL subrectangle 
pow a 
. He atti : Vawicma 
[A, , A, 3! x lw s Ora! for the active device 
conductance. 
me. 2° re : Ce te : ree 
GiAL , Ue) = the active device conductance determined at 
Kn 
2 : 2 ee 
the) poants.A, and» u.« 
k i 
7 «“ - 2 \ . . . > . . + v © . @e 
GM, (Ay , w,) = the ij-th coefficient of the expansion of 
bed J 
| , " Dc eae 
the device conductance as a function of A and 
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2 2 , 
G(Ay» Wy) = the active device conductance determined 


at the points AS and wee 
GM, 5 (Ag we) = the ij-th coefficient of the expansion of 
the device conductance as a function of noe and 
a2 for the MN subrectangle,. 
Gy. ht, 04) = the spline interpolate in the MN subrectan- 


gle [A ].for the active device 


rAyyy) * [Oy Oy yy 


conductance. 


G = total circuit complementary conductance of the 
Oscrllaeor. 

GA. = device electronic conductance at the n-th harmonic. 

ce = circuit conductance at the n-th harmonic. 

Gy = load conductance. 

Ge = device electronic conductance. 

G,. = shunt conductance “ofithe cavity. 

C = equivalent passive circuit conductance. 

IG, | = absolute value “Of G. at zero RF voltage amplitude. 

Ging (A, 0°) thie spline interpolate in the MN subrectangie 


for the active device complementary conductance. 


Gi (w) =circuit conductance, 


SG abo) 2 Spline interpolate for ce in the N-interval. 
Gp = diode conductance. 

i = i(t): current response of the active device. 
i. = i, (t): externally injected current. 


i, (t) = passive circuit current. 


i: (L) =.active device current. 
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I = peak amplitude of RF current for the constant CW 


case, 

I(t) = complex time varying current amplitude for the 
modulated CW case. 

i(t) = microwave frequency component of the current i(t). 


I(t) = bias component of the current i(t). 


T, (t) amplitude of the f-th component of the RF current. 
I .(t) =complex conjugate of Tp (t). 


T5641 = 00d harmonic component of I-(t). 


I5¢ = even harmonic component of I,-(t). 

Ty = complex current Sree at the fundamental 
frequency. 

Lyz, (t) = fundamental frequency component of the current 


response of the active device in the KL rectangle 


[A,r Ap yy] x [wp Hp]! 


Lic (t) = passive circuit current in the L-interval. 
8 = fundamental frequency component of the device 
currenG: 
T 4 = dc bias current. 
c 
Jo = dc current density. 
Jing 7 avalanche current density. 
ws = external current density. 
ind 


Jy J = electron and hole current densities (A/om*) . 


k = Ch.IV: Boltzman constant (1.38 ml 072 


Ch. V: constant = 9.91 x 1029 POD. ois. 


Joule/°K) 
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L = Ch. II: passive circuit inductance 
Ch.III: equivalent inductance of the highly 
| Simpl hied~osciilator circuit. 
eae linear differential commutative operator. 
L(jw,) = function operator in JW, for the CW case. 
L(jw+ oe) = differential function operator to operate on 


i(t) for the modulated CW case. 


Ly = equivalent inductance of the cavity. 

L, = equivalent inductance of the coupling iris. 

L, = equivalent inductance of the device. 

La = equivalent passive circuit inductance. 

L = plunger position with respect to the reference plane 
AA', 

Lon £4 = physical length of channels No. 1 and No. 2, 


respectively. 
L = Ch. VII: package inductance 


Chavillice Des sung va inpute i acne cance. 


LF = low frequency. 
Mn = mutual inductance for the n-th harmonic. 
N (2) = linear differential commutative operator. 


N(jw,) = function operator in JW, for the CW case. 


N(ju+ <p) = differential function operator to operate on 


v(t) for the modulated CW case. 


ie. -3 
Np? Ny = donor and acceptor densities (cm ~). 
eG) = PUpUr Ley econeentrar@zon.in.~(341L,3). interval, 
p,n = hole and electron densities (cm™ >). 
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incident power. 

output power. 

CW oscillator power. 

electronic charge (C). 

quality factor. 

response of the microwave passive circuit. 
microwave frequency. 


equivalent device resistance. 


equivalent resistance representing loss in the 
Cavacy . 

load resistance. 

resistance of unswept region of the Impatt diode 


at bias voltage WES ee 


R = real part of the diode terminal impedance. 

RG = return gain. 

Rj» Ro and R3 = resonator parameters - Ch. VII. 

Ry = bias current measuring resistor. 

R. = pulse generator source resistance. 

Rn = load termination. 

Sw (t)] = susceptance operator for the modulated CW case. 

ees yt Aap U 

S55 = hybrid tee scattering parameters. 

S55 = erbhetj-thcoerticrent of the bicubic Spline fitt- 
ing for the active’ device susceptance. 

See facnucn = spline interpolate in the KL subrectangle for 


the active device B/w. 
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SM, 5 (AQ we) = the ij-th coefficient of the expansion of 
21 


the device B/w as a function of A and aed Lor 


the KL subrectangle. 


Sumy (AoW) = the spline interpolate in the MN subrectangle 
ane actice device B/w. 

S(AC wi) = the active B/w determined at the points AC 
and we 

SM, 5 (Aye uy) = SM, 5 (AQ, we) ; at M= k and N = L. 

c = time. 

Ak = absolute temperature. 

toa scattering limited velocity of the charge carriers. 

Vv = v(t) = applied terminal voltage. 

V(t) = complex time varying voltage amplitude. 

v(t) = microwave frequency component of the voltage v(t). 

NOs ie bias component of the voltage v(t). 

v*(t) = complex conjugate of V(t). 

V = bias voltage. 

VRF = RF voltage amplitude of the n-th harmonic. 

si - CO eine sropaaaMmpLityloe at tC =.:0), 

Voias * Ve 

Vy = dc breakdown voltage. 

ee = the voltage drop across the avalanche region. 

Vie = dc voltage output from the network analyzer pro- 
portional to the phase angle of the reflection 
coefficient. 

Vy = voltage amplitude of the incident wave at plane AA, 

Yop = RF voltage across the diode. 
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Ww = width of the depletion region. 
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X (A) = equivalent device electronic reactance. 


x = width of the avalanche region. 

X-det= response of the detector in H-arm of the hybrid 
tee. 

BL, , AL = LEeSOnactor parameters. 

Xp (WV) = reactance of the diode at bias voltage V. 

XY, J) = passive circuit CW admittance. 


¥,[jw(t)] = admittance operator for the modulated CW cas 


' ae A abe i 5 

ly = admittance characterizing the general model of 
the actice device. 

Y5 = admittance of the simplified model of the active 


device. 

Bf = passive circuit admittance. 

2 = total circuit complementary admittance of the 
escillator. 

ry, = total circuit admittance of the oscillator. 

Yp = diode admittance. 

Y¥(short) = admittance of the short circuit. 

Y¥-det = response of the detector in the E-arm of the 


hybrid tee. 


Yun = the active device admittance determined at the 
o) ; ? 
points A, and wy. 
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= the active device admittance as a function of 


ne and ne 


device electronic admittance. 


transformed circuit admittance as seen by the 


packaged diode. 
device electronic impedance. 


circuit impedance. 


device impedance at zero RF voltage amplitude. 


characteristic impedance of the transmission 


line. 

normalized impedance at plane 
the diode. 

impedance of the active layer 
at, bias, voltage. V,. 

measured normalized microwave 


AA' looking towards the diode 


AA' looking towards 


of the Impatt diode 


impedance at plane 


< 
ats iV, Vp: 


Ch. II: positive constant characterizing the 


device. 


Gho lll: oscillation, growth: factor. 


Che Vo 2 sattenuatlOon coctricient: of the transmission 


line. 


ionization rate of the electrons. 


ionization rate of the holes. 


average ionization rate of charge carriers. 


derivative of the ionization rate with respect to 


the electric field, 


phase coefficient of the transmission line. 
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positive constant characterizing the device. 
diode voltage sensitivity. 

reflection coefficient. 

measured reflection coefficient. 

retlection coetficient at plane AA . 
incremental change of the diode resistance for 


eral = Ma’ 


a voltage change AV = 

change in the diode capacitance between voltages 

V5 and Veale 

change in the bias voltage (View - V5). 

error parameter controlling the iteration proce- 

dure for the evaluation of the steady state 

solution. 

dielectric) constant (F/cm). 

transit angle. 

diode wafer resistivity in the interval j, jtl 

(ecm) 

Ch.III: constant phase angle of v at t = 0. 

Ch. V : phase constant of the reflection coeffi- 
cient in degrees. 

Otherwise = ¢?(t);. time dependent phase of the 
applied RF voltage v(t). 

constant phase of the CW driving force. 

quasistationary frequency. 

Kurokawa: time and phase dependent angular 

frequency operator. 


Ch. Til: angular frequency of. oscillations. 


Otherwise: angular frequency. 
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CW angular frequency. 
instantaneous angular frequency. 

=angular frequency operator. 
angular frequency at the start of the L-interval. 
resonant angular frequency of an LC circuit. 
real part of the complex angular frequency w. 
avalanche resonant angular frequency. 

= hole and electron mobilities (cm*/v.s). 


guided wavelength. 
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CHAPTER I 
INTRODUCTION 


Microwave one port active semiconductor devices 
are finding wide applicability in the output stages of 
many communications systems requiring RF power levels 
Gt the erdersof. SotowloOu.watts.. crghsreiwability, long 
life and extreme compactness are inherent characteristics 
of semiconductor devices and therefore fully solid state 
microwave systems are a considerable improvement over 
existing tube or hybrid types[1,2]. In the communications 
context, the two most important solid state devices are 
theslmpatt and the.uGunn:diodes.,. .In,particular,,the 
Impatt diode, capable of delivering up to 12 watts at 
6 GHz [3], is the most powerful solid state CW microwave 
source presently available [4]. 

A significant area of growth in the communications 
industry is in computer based high speed data transmis- 
sion, where digital microwave systems play the major 
role [5,6]. Design trends for these systems indicate 
the strong possibility of the Impatt diode being used 
even more extensively in the near future. The modula- 
tion schemes most likely to be adopted for digital 
communications are phase and frequency shift keying in 
which the problem of amplitude nonlinearity inherent in 


solid state devices is not significant [7-9]. 
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It is unlikely that the Impatt diode, or any other 
solid state device will immediately replace tube devices 
in the near future. Rather, solid state technology will 
be applied selectively to perform specific system func- 
tions while freely capitalizing on the advantages of an 
all solid state approach. For this reason, and to ensure 
successful implementation in new system concepts, a 
precise knowledge of the characteristics and limitations 
of these solid state devices is essential. 

The electrical characteristics and properties of 
microwave circuits are governed by the active device. 

In system applications,either in oscillator or amplifier 
circuitry, the interaction between the active device and 
the circuit is very important for design considerations. 
Under transient conditions, the dynamic characteristics 
of the circuit account for this interaction and yield the 
transient response of the circuit. 

To optimize efficiency and economy, which are 
key figures of merit in power applications, the active 
devices are strongly driven into the nonlinear region, 
Consequently any detailed description of the dynamic 
hehaviour of active circuits, which is of primary 
importance to the designers, is very complex, and 
requires extensive numerical calculations. The numeri- 
cal methods are usually very different, depending upon 


the tvpe of active device. To simplify the calculations, 
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either the active device is very accurately described 
while only simple external circuitry is assumed, or in 
the opposite approach, external circuitry is accurately 
described and combined with an idealized negative 
conductance device possessing Van der Pol's cubic non- 
linearity [10]. The latter approach can only be applied 
to the simplest of these devices, such as the tunnel 
diode. On the other hand, high power devices have non- 
linear dependence of both conductance and susceptance 

on RF voltage and frequency. Consequently they are poorly 
approximated by a cubic nonlinearity. Hence in dynamic 
studies on these circuits a more accurate method based on 
detailed characteristics of the given active device must 
be employed. 

The presentation of this thesis starts with the 
development of a general model characterizing the dynamic 
and steady states of nonlinear microwave ieee ieee ue. 
The application to free-running oscillators is then 
described. Chapter III develops the basic theories of 
negative conductance solid state oscillators in the form 
of a simplified model. The basic principles of Impatt 
diodes, their structure, DC characteristics and operation 
under small and large signal conditions are discussed in 
Chapter IV. Chapter V describes the measurement tech- 
nique and method used to evaluate the large signal admit- 


tance of Impatt diodes together with their embeddin 
p 


, ie td . ae 
bg Pr etek x Los EAU S wae 7 \ ae 


ay 


yaxpode- St. Yideoz tS, ‘an SER 


an 1 
] 4 . _e ; : ¥  — : 
baweeean pt ytd Goa ee 
"a 1y ia 


* ¥ 7 
~ tT — 7 
‘i a a ' Pi a] pete , 
* . = E = 
- i = ye 
ea <! yee Ah 
a es 
~ 4 ’ 
SORE Shia - Sy LyY 
* 
= ¢ } 5 
_— 2 oF, 
: S,, a € * 4 
- a 
, #% a oe C.e a 
» A] 
¥ 7. iad ~ “ me 4 -' ~ <a 
_—e ‘ - (ot ee 
| yr 
Pn te 5 “Me > | * 
“ _ 
oJ 
- r “| 4 “ 
A OS op ht S27 AD, ee 
« > =a - - fe ba 
2 3 +“ = é V 4 ¢ 


MP iegt a 
bevel qa 
bas a ii 
» 
S212 a609 > 2S J 245-810 td rie eee mee -- 4 
= 7 » <= 7 ) Lo - 7 


- 


4 Sesio Leben Lsgetsy. & Sp aeSmaeer 

va - A 

7 

Mt  %es gh “ip cote nde 
ty “Tee OR at thee of nobswok Lam 
neds ofzed sid doékaveb IIT sezanne 
«ago steta bil22) eo rtoukaes oly 
i bad v 2 

, 7 
ms ee 
Sor. 2 QsaSw Sift enon Sea Sh odeperet 


' : a 
20 AHENOR STARS JS qRaaOUS 22 sag nah 


Wo2th wget 2notsihn08 ‘Langts” 4 ron Bie 
+ t4imeeuebor edd padlsoesh V -neaaaitg 
Ap aa 

$j pel aifa inca nit9: og tee ; 


Circuitry. A new frequency discriminator to measure the 
instantaneous power and frequency is discussed in 
Chapter VI. The model developed in Chapter II has 

been used for a theoretical study of the transient and 
steady state behaviour of the oscillator in Chapter VII. 
The results are applied to the Impatt diode osciliator 

in particular, using the large signal parameters deter- 
mined in Chapter V. The following chapter deals with 
experimental studies on Impatt diode oscillators, and 

the results are compared with the theoretical predictions 
from the previous chapter. This chapter also discusses 
the operation of a bias modulated oscillator and develops 
a digital modulation scheme for transmitting data at 
rates up to 40 Mb/s. The concluding Chapter IX summarizes 
the findings of this study, and presents suggestions for 


further work. 
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CHAPTER II 


METHODS FOR DYNAMIC STUDIES ON MICROWAVE CIRCUITS 


EMPLOYING ONE-PORT ACTIVE DEVICES 


Ziv tNCrOanctLLOn 


Microwave active solid state devices are usually 
very small compared to the guided wavelength so that the 
electric and magnetic fields do not vary significantly 
along their dimensions. Therefore, conventional circuit 
theory can be used in the analysis of these devices even 
up to the millimeter wavelength range. The device res- 
ponse to a single frequency excitation, which may be 
either RF voltage or current depending on the embedding 
Circuitry, contains an appreciable number of harmonics 
the amplitudes of which depend upon the signal strength. 
The higher harmonics, however, are greatly attenuated 
by tuning the high Q embedding circuitry to the funda- 
mental frequencv. This enables the active device to 
be uniquely characterized at the signal frequency by 
a nonlinear immittance dependent upon the signal 
amplitude and frequency. Consequently, equivalent 
circuit representation can be employed when investigating 
the active circuit response to unmodulated CW signals. The 


steady state amplitude and phase of the device response at 
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the signal frequency can be obtained by solving nonlinear 
algebraic circuit equations. However, in many applica- 
CL1ONs ,esome aqditronal information is vital, e.g. 
stability and noise properties of the steady state 
solution, response of the active circuit to modulated 
Signals (in particular nonlinear and phase distortion, 
AM to PM conversion) and intermodulation characteristics of 
the device. The desired method must be suitable for dynamic 
studies on active one-port circuits exposed to both CW 
and modulated signals in order to give some insight into 
these problems and therefore facilitate better circuit 
designs. 

The oldest method for dynamic studies on one-port 
devices is that developed by Van der Pol [10] for nega- 
tive conductance oscillations. He approximated the 


current response of the device, i, by 


ta=yoav + yee (251) 


where v is the applied voltage and a,y are positive cons- 
tants characterizing the given device. Representing the 
embedding circuitry by a single tuned parallel resonant 
circuit, whereby all package parasitic and device suscep- 
tances are included in the embedding circuitry, Van der 
Pol arrived at the following second order nonlinear 


differential equation: 
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where L, C and G are positive circuit parameters and i; is 
the externally injected current. To solve the above 
equation, Van der Pol proposed a method of slowly vary- 
ing parameters which was later refined by Bogoliubov [1l]. 
According to this method, the solution is sought in the 
form 


lO eas = a lh i gt IT a es RO Mile C23) 


Thus, the rapid time variation of v(t) is assumed to be 
harmonic, being characterized by a carrier or quasi- 

Stationary angular frequency @® , and the slow variations 
of v(t) are described by the unknown functions A(t) and 


o(t) which satisfy the relationships 


ieee ee doen 
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In other words, A(t) and ¢(t) are functions which are 
Slowly varying in time intervals comparable to 27/w . 
Substituting, Boe 2eo) intoekq. (2.2), two first order 
differential equations are obtained from which A(t) and 
o(t) can be conveniently evaluated using numerical 
methods. 

There are a few serious short comings when apply- 
ing Van der Pol's approach to active solid state devices. 
First, solid state devices are characterized by non- 
linear conductance as well as susceptance; Kuno [12] has 


recently extended the method by including nonlinear 
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inductance . Secondly, because of the single tuned 
circuit representation of the embedding circuitry, the 
active solid state microwave circuit, which usually 
Operates in a wide frequency band, cannot be well des- 
cribed. Finally, the active device is described by means 
of various constants, e.g. y, which are rather general 
and cannot be accurately determined for a given device. 
Many present dynamic studies are based on a method 
developed by Kurokawa [13] in 1969 for noise and stability 
investigations on CW oscillators. Kurokawa expresses the 


diode voltage and current in the form 
v(t) = A(t) Refexp jlot+ o(t)]}. 


He does not restrict the passive circuitry to a single 
tuned model. Instead, he uses a generalized immittance 
expression in which the w-factor is replaced by the 


operator 
tee eaee | 2 at: 


and the active device is represented by its large signal 
immittance. Expanding the immittance expression around w, 
Kurokawa derives two first order differential equations 
for o(t) and Alt). “Since the coefficients in these 
equations are related to the device immittance, this 
method can describe a particular device accurately. 

The great potential of this method lies in the studies 


on steady state CW operation and small noise disturbances. 
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Since the device nonlinearities are not fully taken into 
account, this -mevtnod' 1s hardly suitable for/the study of 
active devices exposed to strongly modulated or transient 
Signals. 

In 1972, another method was proposed by Gustafsson, 
Hansson and Lundstrom [14]. Representing the active cir- 
cuit by a feedback model containing a nonlinear element, 
they were able to apply a describing function technique. 
Since this technique has been intensively used in the 
study of nonlinear control systems, considerable experience 
is made available for the study of nonlinear problems at 
microwave frequencies. When dealing with strongly modu- 
lated Signals, this technique is only slightly better 
than Kurokawa's approach. 

Se€crions) 2. atone. OL this chapter deal with sine 
development of a theoretical technique for the charac- 
terization OL Nonianear mLcrowavVe active Circuits, nis 
technique is applied to free-running oscillators in 


section 2.6. 


2.2 ) General circuit, approach 


When studying the dynamic behaviour of a micro- 
wave circuit in which is embedded an active one-port 
device, it is very useful to divide the circuit into a 
passive one-port circuit and an active one-port device, 


as shown in Fig. 2.1. In this manner it is possible to 
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PASSIVE | ACTIVE 
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Figaz..) Representation Of microwave Circuit 
by two separate passive and active 
sections. 
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separate the characteristics of the nonlinear active 
element that vary relatively slowly with frequency 
from those of the passive circuit which vary rapidly 
with frequency. 

Assuming the terminal voltage, v(t), to be the 
driving force, the operation of the complete active 


circuit is described by 


igtv(t)] ef: iptv(t)] =2i0 k (2a) 


FE the active circuit is disturbed by an externally 
injected signal, is, the above relation, which is essen- 


Clavly —“Kirenhort'securrent Law, will change’ to 
oe AG a on ee (27.6) 


If the current response of the network to the 
applied terminal voltage can be determined, the voltage 
at the terminals of the active device, excited by any 
externally injected signal, can be evaluated using Eq. 
(2.6). Knowing the terminal voltage (at the terminals 
of the passive circuit), the current, flowing into the 
passive circuit, i,(t), can be evaluated. Using linear 
circuit techniques, the current and voltage anywhere in 


the passive circuit can thus be determined. 


2aeePassive iCLrCuLe characterization 


In analogy with low frequency networks, the 


driving force, f(t), and the corresponding response 
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of a passive microwave circuit, r(t), are related by 


[15] 


where Agreee ra Doreee rd, are characteristic parameters 


nN’ 
of the given circuit. Depending on the circuit excita- 
tion, the terminal voltage can be either the driving 
force itself, as is assumed in this study, or the circuit 


response. Eq. (2.7) can be written in operator notation 


as 


d Hh d 
Tis) tatty =8 Nae te tats (27.13)) 


where the linear differential commutative operators, L 
and N, operate on the functions enclosed in the { } 
brackets. 

The circuit response to an unmodulated CW driving 
Force can be evaluated from Eq, (2.8) by Fourier trans- 
forming it into an algebraic equation. The steady state 
amplitudes of the driving force and the circuit response 
are then related through a frequency dependent circuit 
immittance. 

A single frequency CW voltage excitation 

ek 

f(t) = v(t) =A e ° (2/9) 

produces a current response r(t), which will have the 


general form: 
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r(t) = i(t)= I,exp (jw,t) = exp (jut) ¥, (jw JA, (2.10) 


The actual physical signals are given by Re{v(t)}and 
Re{i(t)}. The term XY, (9%) 4 called the passive circuit 


CW admittance, will be given by 


° aS N(jw,) 
XY, (94) = AS = L(ja) C211) 


where N(jw)) and L(j4,) follow directly from the Fourier 
transform of Eq. oat 

For a modulated CW driving force, the evaluation 
of the circuit response is relatively complicated. Then 


the dirving force is usually given by 


f(t) = v(t) = A(t) exp] (wt + o(t)) = Vi(t)exp jut 


(Zeken 


where wW is the quasistationary angular frequency. The 
Signal V(t) is the complex time-varying voltage amplitude 
which has a magnitude A(t) and a phase angle ¢(t), where 
both A and ¢ are real quantities. 

| The general circuit response can be expressed in 


a foOrmysimilar to Eq. (2.12). apse. 


- 


The operator argument of N and lL changes depending on 
whether one is treating modulated or CW signals, or on 


whether one is operating in the time or frequency domain. 
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Tilt dares it) elt) expay at (2 13,) 


where I(t) is the complex, time-varying current amplitude. 


Dig rerentiating Eqs. (2.12) and(2.23), we obtain 
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eo ovAt) sexp jt [50 + $1 tv (t) | (2.14) 
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obllCs ies Wn ape fie nre mals 
Sony gee, TEL UT ee ; 


Thus the derivatives of v and i for the modulated frequency 
case and for the constant CW case can be symbolically 
expressed in the same way if one replaces the algebraic 
factor (j0)of the latter by the differential .factor 

(ju + es which then operates on the slowly varying 
complex amplitudes of the non-constant CW quantities. 


SubseuLcucing wos.) 2.4) wand) (2050) aintol) Haw) (2.7), and 


Utivazung =the operator notation’ of Eq: (278) , we’ obtain 
Lijo+ Sf (t)}= wijo+ Spit} (2.16) 
an de dt ; 


Analogous to the CW case, we can define a linear 


admittance operator, Y [ju(t)), for the non-CW case by 
. nN a. ’ aa aes Ow ak, 
¥[jw(t)] Y [jut ae] —— : (Zn) 
The complex amplitude of the current response is then 


given by 


I(t) = ¥ [jw (t)] {A (t) exp So Ch ly .). (2518) 
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The admittance operator Ytjw(t)] operates on the term in 
the { } brackets in this notation. We can also define the 


frequency operator 


[w(t)] = [a - 5 I (2.19) 


where w(t) is the instantaneous frequency, and w is the 
quasistationary frequency. The instantaneous frequency 


Ory (2) Ss 


w(t) = (ott o(t)ead SR... (2.20) 


In many microwave applications, the amplitude A(t) 
and phase angle ¢(t) usually change insignificantly within 
one RF period, i.e. 


GA ee Ga 


: dd 
ac ; at oe (Ly e (2-427) 


Accordingly, the admittance operator can be expanded into 


a Taylor series about w: 


n sae n 
Se Sa (2.22) 
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Substituting Amen g (=) Pron wt 2 ar Oo) win Cobos 4 (2 -522))) sand 


using the result for Yo img» (2.18), we have 
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Ete] (jw(t) - ja | A(t) exp j¢(t) } 


(2525) 


Equation (2.23) is known as the Carson-Fry series [15]. 
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The convergence of this series was discussed by Stumpers 
[15] for FM modulated signals. Stumpers pointed out 
that any terminated series approximating Eq. (2.23) 
must be discussed together with the remainder since the 
complete series may not be convergent for frequencies 
far removed from w. 

Inspection of Eq. (2.23) reveals that the operator 
XY [jw(t)] will be composed of terms jw, (ju4; (ju)>, 
GEC. Separating the real and imaginary parts of 


Y[jw(t)] results in 
Yo [jw(t)] = Glw*(t)] + jutt)Sfw* (t) ] (2.24) 


where GL (t) J and S fu aGEbel are conductance and capaci- 
tive (plus sign) or inductive (negative sign) susceptance 
operators, respectively. In Eq. (2.24), the Hien is 


an operator which, by using Eq. (2.19), is given by 


eS (a a em nA Daan (2.25) 
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Expanding G burnal and Slee lee by a Taylor series about 
uae we Obtain, using sha. (2.25), 
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Substitution of these results into (2.24) yields the 


desired passive circuit non-CW admittance operator: 


+ ju(t) — 
eae dia - 


a eee Dee 2 -2.n 
ted Gey l= J ete aS (w L| a Ce }e=-we,) 
=0 
(2228) 


Then the current amplitude is given by substituting 
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The convergence restriction applying to Eq. (2.23) also 
noldserpGohG. (2.29)... Since the frequency deviation 
(usuaily a few percent) at microwave frequencies is 
usually within the series convergence interval, a finite 
number of terms will approximate the current amplitude 


with the desired accuracy for practical purposes. 


2.4 Active. circuit model 


The relation between the terminal voltage and the 
current response of a nonlinear active circuit can be 
obtained from Eq. (2.7) by allowing the coefficients 


Bese a and Oe, De cOedependa on the driving force. 
fe) n O p 
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Expressing this dependence as a power series, Eq. (2.7) 
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C2r 
where Aah and ooh are constants depending both on the 
nature of the embedding circuitry and the active device. 
For H = 0, the above equation simplifies to the relation 
descrabing?passive icircuits.(Eq.—(@2.7)) and for M = 0, 

P = 0 and Aah = 0 for h>0, the above relation is 
reduced to a generalized Van der Pol's equation for 
active devices. For nonzero M, P and H the current 


response of the active circuit is generally a nonlinear 


function of the driving force, as well as its derivatives, 


ana “ts “integrals. 
Unlike the case for the passive circuits, the bias 
components here cannot be disregarded, and the terminal 


voltage and current will have the form 


vie) =) A(t) o+ v(t) (2 SA5) 
and 
aot) oat oe ah Z(t) (25732) 


where Vy (t) and T, ft) are the slowly varying bias com- 
ponents and v(t) and i(t) are the microwave frequency 
(RF) components. The time dependence is assumed implied 
in the following derivations (e.g., Ty (t) = Loreees etc.) 
for the purpose of simplicity of the expressions. Subs- 


Piet roe (oh) oldemds «(25520 anto, ba. (2.30) yields 
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dre 3: 
By expanding the above summation in (vote), Collecting 
like terms in wo and then separating the bias and LF, 


and the RF components, Eq. (2.33) yields 
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In case of low efficiency two terminal devices, 
e.g. Impatt-, Transferred electron= or Barrit diodes, 


the rectification components, expressed by the [ lr 


Lelmo ii bd. (2.54) are instaniricant compared Lo the 
bias terms. Then only the bias terms, as expressed by 
M any 


Pp 
(a Ae ee as 2 (286) 
seen ae p=0 PO atP 


describe sufficiently well the dynamics of the bias con- 
ditions, ,Only for very large signals, can the rectifi- 
cation and time-dependent thermal conditions affect the 
bias components of active solid state devices. Because 
of large time constants in the bias circuit and 

large thermal time constants of the mounted active device, 
the changes of the bias conditions in microwave amplifiers 
and oscillators are very slow as compared to the signal 
modulation. In oscillators specifically designed fow 
bias modulation, the almost step-like changes of the bias 
voltage are usually so short that most of the oscillation 
transients take place under nearly constant bias condi- 
tions. Moreover, since the signal modulation has hardly 
any effect on the bias conditions, the dynamics of RF 
signals in active circuits are usually well described 

by 


M H es; my Pp H ’ p~ 
[Ty ae SS} = [CS be 2] 
m=0 h=0 at) RF p=0 h= dtP Jrr 


(2737) 


where the coefficients aah and b are time independent. 
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The coefficients in the above differential equa- 
tion depend on the terminal voltage (driving force). 
In many applications (negative conductance amplifiers 
and oscillators), single carrier frequency operation 
(as to the power and energy flow) is the desired opera- 
tional mode. Hence, a single carrier voltage, either 
externally excited or self-generated, is the desired 
voltage waveform at the device terminals. In practice, 
however, the terminal voltage is composed of a number of 
harmonic frequency components. Their amplitudes depend 
upon the amplitude of the corresponding components of the 
nonlinear device current response and the impedance of 
the embedding circuitry at that particular harmonic 
frequency. Because of the very nature of microwave 
circuits, the impedance of the embedding microwave 
circuitry does not decrease monotonically to zero with 
increasing frequency. However, except, for example, in 
oscillators with tuned harmonics, this impedance is 
usually insignificant at harmonic frequencies. Very 
often the device capacitance acts as a low-pass filter, 
drastically reducing the higher harmonic voltage com- 
ponents at the device terminals. Then the single carrier 


driving force 
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very well approximates the operating conditions of many 
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active devices. The quantities A(t) and ¢(t) are slowly 
varying functions of time, obeying the restrictions given 
DYeEG. (2.21). No other restriction is imposed’ on these 
functions. Thus, a large variety of signals are des- 
cribed by Eq. (2.38a), e.g. frequency and/or amplitude 
modulated signals; single or multiple frequency modula- 
tion; transient signals and two or more signals with 
closely spaced carrier frequencies. In some cases it 

may not be possible to separate the slowly varying part 
of the signal into its amplitude and phase components. 


Then the driving force has to be expressed by 


x A is Ge oe -jot 
a eee V(t) e ; V*(t) Ee 
= Re[V(t)] cos wt+ Im[V(t)] sin wt (2436p) 


where V(t) is the complex time-varying voltage amplitude 
and the following derivation should be carried out for 
the above driving force. 

In theoretical studies of active microwave circuits, 
single frequency excitation, while considerably simplify- 
ing the mathematical treatment, is usually assumed. Only 
one exponential term (the one used in the passive 
circuits) is considered in almost all studies. In noise 
and stability investigations, this does not introduce 
any significant error. However to properly describe the 


nonlinear operation of the active device with strongly 
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modulated signals, both the exponential terms (which 
express the real driving function) have to be considered. 
No harmonic restrictions are imposed on the circuit 
response, which is expressed in the form 
F 
Tee ght Ct 


£40 
-F 


yeifot (2.39) 


Here, ) denotes the summation of all terms with the 
£40 
exception of the bias term, i.e. £40, and the complex 


amplitudes which satisfy 
I, (t) = Ti, (t) (2.40) 


carry information on both the amplitude and the phase of 
a particular current component. 

The summation in #" from Eq. (2.37) is expanded 
using Eq. (2.38a) “and if Iike-terms are collected, the 


; : ~h ; 
summation inv can be written as: 
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Similarly, differentvating Eq. (2.39) m times and Eq. 


(2.38a) p times yields 
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PaAnEst (2.41), the summations are carried out to the 
highest integral number H/2 or = and the coefficients 
oF ad tard with different g's are related through simple 
algebraic relations to the coefficients anh OfSEqaulans7 >. 
Substituting Eqs. (2.41)-(2.43) into Eq. (2.37) 


yields 
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The ane set of nonlinear differential equations 
for complex current amplitudes describes the active 
Gevice response to a sSingle-carrier driving force with 
non-constant amplitude. Even if all the parameters a 
and b were known, the complexity of the above description 
would make it unsuitable for any numerical studies. The 
description of the active device response becomes simpler 
if all passive circuit elements between the device ter- 
minals and the active device itself, namely all device 
package parasitics, very often involving energy storage 
elements, are included into the external pasSive circui- 
try. The transformation of the nonlinear device quanti- 
ties into the terminal quantities is then eliminated. 
Furthermore, since the external circuitry is assumed to 
present an infinite admittance at the device terminals 
Giealt trequenclessovher than ,thestundamental frequency, 
the current of this simplified active circuit does not 
contribute to higher harmonic generation. Thus 


eer i 0 LOL youl (Ue aren |e 0 (2.45) 


and the equation set (2.44) reduces to a set of uncou- 
pled differential equations with constant coefficients. 


Thus Eq. (2.44) reduces to 
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If in the above equation, one considers the current 
amplitudes of odd and even harmonics separately, then 


for the odd order one gets: 
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Equations (2.47) can be regarded as generalized 
admittance relations describing the current response of 
an active device to a Single carrier driving force with 
nonconstant amplitude. The amplitudes of the current 
response are nonlinear functions of the applied vol- 
tage amplitude, its derivatives and its integrals. 
Knowing both sets of the coefficients, i.e. a and b, 
the current response of the active device could be 
evaluated for any amplitude and frequency of the driving 
force. Disregarding the fact that it may be very 
difficult to find the appropriate sets of coefficients, 
the above description of the active device response is 
still extremely complex. It would lead to integro- 
differential equations, and is hardly useful for numeri- 


cal studies. The active device description can be 
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Significantly simplified by the following choice of the 


a coefficients 


a =]; auirm <0? forsi ame a0l ss (2.48) 


Then the square-bracketed [ ] differential operator in 
Eqs.(2.47) reduces to unity and the components of the 


current response are given by 
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In this simpler description, the current response of an 
active device is a nonlinear function of the voltage 
amplitude and its derivatives. Thus, a given set of b 
parameters will accurately describe the device response 
in a more or less limited frequency range at large 
amplitudes of the applied voltage. Using numerical 


methods based on bicubic spline approximation, both the 
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b parameter evaluation and the current response descrip- 
tion for any amplitude and frequency of the driving 


force can be accomplished in a piece-wise manner. 


2.5) Model. characterization 


The previously described techniques are designed 
for studies on active microwave circuits with non- 
constant CW signals. For a constant CW driving force | 


characterized by 


A(t) =A;  o(t) = o, = 0 (2.50) 


the techniques simplify to a conventional steady state 
approach involving large signal immittance relations at 
the basic frequency component. Thus by substituting the 
values of A(t) and CGE given Dyno. (2.00)) Into. Bq: 
(2.47a) and rearranging terms in powers of us and ae ; 


the admittance of the general model characterizing an 


active device is given by 
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where qT is the complex current amplitude at the 
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fundamental frequency and, according to Eq. (2.49a), the 
admittance of the simpler active circuit model is given 


by 


where the summation is carried out to the highest integral 
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case, the admittance of the passive circuit is given 
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In Eqs. (2.51)-(2.53) the d coefficients are related to 


the b coefficients by 
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Since they are described by linear differential operators, 
passive circuits are uniquely characterized by a single 
set of coefficients which can be valuated from circuit 


admittance measurements. Active circuits are described 
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by nonlinear differential operators, which are not commu- 
tative, and thus two sets of coefficients uniquely des- 
Ccribe active circuits. In analogy with passive circuits, 
the large signal admittance of active circuits provides 
one set of matching conditions while the other set of 
conditions should be provided by some other nonlinear 
characteristics of the circuit. Even if two proper sets 
of conditions are known, the evaluation of coefficients 
may be extremely difficult. Based on theoretical studies, 


one can sometimes guess one set of coefficients, e.g. 


in the case of Impatt diodes a5. i 0<a,< 1 and ay = 
Stat ee 0. A single set of unknown coefficients then 
remains. In amplifier and oscillator applications, the 


large signal considerations are of primary importance 
while the nonlinearities associated with higher powers 
of derivatives of the terminal voltage are of minor 
importance. In these applications, the simpler model of 
the active device, for which the single set of coeffi- 
cients can be determined using the device large signal 
immittance, is more than adequate for dynamic studies. 

If the passive or actice circuit admittances are 
analytically specified, the differentiation of the 
admittance expression can yield the model coefficients. 
Most often, however, the admittances are specified by a 
series of measured or calculated data points at different 
Signal frequencies and amplitudes. Even though numerical 


differentiation can be used, the natural cubic and bicubic 
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Sprune Curve fitting (into the: characteristics of a: given 
passive and active circuit) is most efficient. 

The natural cubic spline, the first two deriva- 
tives of which are continuous, was shown to best appro- 
Ximate the curve joining data points with minimum curva- 
ture. Given*a*set of M*data points,” e.g-, be Pieces bes a 


computer routine returns the spline coefficients in each 


of the (M-1) intervals. The spline interpolate in the 


L-th interval, (ae - men) can then be constructed by 
4 m-1 
Ze), 2 2 2 
a yee ee eb Com we) : (2555) 


m=1 


More than three data points are needed for the evaluation 
of the spline coefficients. The approximation accuracy 
increases with the number of given data points. The 
coefficients in a given interval, which are related to 


the first three derivatives of the approximation func- 
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tion at the interval boundaries, depend on all data points. 


Thus the effect of nonsystematic errors of individual 
data points on the numerical derivatives is inherently 
reduced resulting in "smooth derivatives" [16] (this 
property is partially attributable to the best approxi- 
mation characteristic, and partially to the minimum 
curvature property of the cubic spline). Unlike a Taylor 
series approximation, the question of convergence is 
irrelevant here. Since all spline coefficients are 


equally important, none of them can be neglected without 
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seriously reducing the approximation accuracy. 
The natural bicubic spline is fitted through Mx N 
functional values, e.g. G ge oraigd & , corresponding to 
Lage. M,N 


EWO, Sets Of Grid points, e.q.-A reeerAne Wy reser Wye The 
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computer routine returns the coefficients of the bicubic 
Spline interpolate on the KL subrectangle [A,r Any] 


[w ] in a 4 by 4 matrix. The spline interpolate 


ee ee 


is then given by 
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To evaluate the d coefficients, the device admit- 


tance data Yen should be first expressed as 
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Finding the spline interpolates for the G and S functions 
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and equating the coefficients from Eq. (2.58) and Eq, 
(2.52) at various powers of A and w, the da coefficients 


can be numerically evaluated. Since the b., Hives and 
Ee sety 
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b, h.2 coefficients are linked through simple algebraic 
e , 
relations, they can be evaluated from the known d coeffi- 


cients using the formulas 
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Se2nrl _/2h+1 
Oe eee oh os ee pe oh, © (2.59) 
and 
a -ly Eyceaaie 
Corl Ziweie@2oT lenge Sth 2pt leis O: 


The fundamental frequency component of the current 
response of the given device to the nonconstant CW driv- 
ing force which is given by Eq. (2.49a) can then be 


analytically evaluated from 


. es 
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5 eee t)el (-J jw tag) 15 ») | 
(2.60) 


where the b coefficients are related through the numeri- 
cal spline coefficients to all admittance data points 
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under nonconstant CW driving .OLrce,epy, using: BOs 4 (2.13) y 
(Ze.7) and (2.18), the passive circuit cttrrent may be 


analytically evaluated from 
, nv jwt ; d jJo(t) . 
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where the b©° coefficients are related through the numeri- 
Calesuline: coer sr CLeEnts OL) thescircuit for the interval 
Ie 

The above technique will be applied for the analy- 
sis of free-running oscillators as these have wide 
applications in communication systems. The resulting 
method can be extended to analyze injection locked 


oscillators and amplifiers. 


2.6 Theoretical description of a free-running negative 


2 a 


Operation of a free-running oscillator is described 
by Eq. (2.5). This nonlinear differential equation 
together with proper initial conditions uniquely describes 


the time dependent terminal voltage of an active device. 
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If the transients in the bias circuit can be disregarded 
ang the vembedding circuitry js designed to ensure ‘single 
frequency oscillation, then Eq. (2.5) will describe the 
transients ofzthe RF oscillations.’ “in this’ case, only 
the fundamental frequency components of both the active 
device and the passive circuitry may be considered. 
Substituting for the device current as given by Eq. (2.60), 
Sno wore Chew pLOper circuit Curvenc) (ha. (2.61)) anto 

Eq. (2.5), the following differential equation for the 
unknown, Slowly varying amplitude A(t) and phase (tt) 


of the RF oscillations can be obtained 
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c . dP JO) eae 
+ op hr celaueae! Pt) } =a Oe: (2262) 


In Eq. (2.62) the b coefficients are related to 
a few discrete admittance values of both the embedding 
Ciorcuiueana, tie matllve Uevice. COMInNING tie cCeéxiis On 


the Jert Nand side oft Eq. (2.62), the differential 
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equation characterizing the slowly varying components 


of the terminal voltage can be obtained 


7 2S p 4 
a iy +2 jo (t) 
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C2505) 
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b = bo = he b =70 2.64 
Kom 7O, oe amie Oso re NlSD O70 Jenlai yO; es. (2.64a) 
and 
b Za fed 5 gees (2.64b) 
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If a proper KL~-subrectangle is used, i.e. the 
instantaneous voltage amplitude and dd/dt are within 


the limits. Ay, < A,(t). < and 0 < d¢/dt 


k <4 S Apay = Sore ala C 


L’ 
then the differential equation (2.63) describes accurately 
the dynamic conditions of a given active microwave cir- 
cuit. Performing the differentiation, separating the 

real and imaginary parts and neglecting the higher 

order derivatives of the products of the amplitude and 
phase, the following set of differential equations 


can be derived in the slowly varying quantities A(t) 


and. ¢(t): 


Sh 


Vie oa a ; ba a 
iy ve Y Ny | bats to 
nan ny ] J) at 
' . A " ri -_ 7 : 
an A 
i i i 
; av . 
oe 
‘ i] ‘ rin 
~ Pe ahs eae ki 
fe of ia ve, he ) wal wa ‘ 
7 ; i iy Piss ft Wyola Hp Eh 
. rr U ‘ ‘ 


' <r, i> pe é 5 mT er ate an ae 
> atredogensy yabyisy yowolm aie: 


— 7 , Aaa 
bbintatde: oc f - 
z) 


ae! wy 


' 
ef 


af » “i 
’ SUpLy nnd Aer a. i 
L : ; y i 7 : io ‘a 7 
Th esha 
: me a “ vi 
Mm, be Vie C+) ge Ye ae ; ¢ 
: A ae 4 , ‘a i 


a4 
| iy iS) aa 4 %?* ; >y% 
ty . : 
2 : ‘ ‘ o La 7 i 
. : ates cooky & 
ay. 
iJ i bea at Obpraspetds TENE negara oat 
aa beim storey tLe ns Oped sen | ‘auton 
is d mea 
7 Y a a 
.! ha Die (ra? (a)4 e a2 b. 
’ a OL: IBUD? ta ignedessib = 
i} ) ral Vos ey ai 
j a < 
‘wvisee ov ty 5 ke anoisibioe ‘obmsingh ei 
| ji “ ‘ i‘ GM att 
saad, (ES phaAmeets tb aid | eeeriahal 
Ny rts 
| ta ‘a iy 4 a Tr ad 7 7s f Py mi bra 2 


Cree aetakiehtes nti, Meith es paoubon esi to. ies is 


oy 


itevips Lol Ldnetat i LB :.20 


$s) hom revere ae slat ine aaa 


38 


6 n 6 n O 6 n-l 
aco y x oe. en xp" +34 } TALI Dn = 0 
n=0 ° n=0 n=l : 
(2.65a) 
and 
n+1 6 n fe) 6 2 n=L 
‘e) ss n x n xOA ntl= 
199 —DG+A) =DBt=— ) =~, D  Ge=0 
n=l 2! fe at 2 2, 1) 
(26.5) 
where in Eq. (2.65) D is an operator representing partial 


differentiation’ with respect to w, x = dd/dt, AC = dA/vaGc, 
G and B are the total circuit conductance and suscep- 
tance Of athe escallator, G and B are the total circuit 
complementary conductance and susceptance of the oscil- 
lator. The two nonlinear equations of (2.65) may be 
written in a form suitable for steady and transient 


(numerical) solutions in oscillators as 
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InfEqs. (2°65) and (2.66), the Y and Y' “are related to 
the passive and active circuit parameters through the 


following equations: 
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In these equations, the Lo parameters characterizing the 
passive circuit in the L-interval are related to the oh 
coefficients of Eq. (2.64a) through the equations given 
by TEq. 2(2054/"'for*h.='0.9 They-are-directly computed 
from coefficients evaluated by fitting the discrete 


admittance values of the pasSive circuit to the one- 
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dimensional cubic spline. The c's are nonlinear func- 
ti0ns of the voltage amplitude A up to the “sixth order 


and can be written as 


3 
_7,\m 2n 5 
Rives oe bis ae nat aap Oa ee ee 
= m 2n 
DQ = ia v 
KL ome, 2 pie 1) ebeames wa fads Oe zie 
3 m 
an (+1) 2n A 
KL" om, 2 deer (antl) KL¢2m, 2n” at at oa 
| ; m 2n 
Kno amtr = 2 (2) (2041) ep Goma anh 1 M= 071,273 


The active device parameters (c§'s) can be deter- 
mined for any subrectangle KL depending on the value of 
the amplitude A and that of the instantaneous frequency 
wtx. The d's, here, are the active device parameters 
and are given in terms of the b coefficients through 
Eq. 7 (2.54) 2" They are calculated from the coefficients 
resulting from the bicubic spline which fits the device 
large signal admittance data. 

The transients of a given oscillator, with known 
d coefficients (passive and active circuit parameters), 
can be studied by solving the nonlinear equations (2.66), 
Simultaneously in A and x. This may only be achieved 
fora given set of initial COonadterons » (for A "andi sx). 


The physically acceptable initial conditions, in free- 
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running oscillators, are those which correspond (i) to 
the start of oscillations or (ii) those corresponding 
to the almost steady state. The solutions give the 
amount of time taken by the oscillations to build up 
to the steady state amplitude and frequency, which is 
an important factor in the design of bias modulated 


solid state oscillators. Because of the complexity of 


calculations, Eqs. (2.66) can only be solved numerically 


on a digital computer. 
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CHAPTER III 


BASIC THEORIES OF NEGATIVE CONDUCTANCE 


SOLID STATE OSCILLATORS 


So) (introduction 


In practice, negative conductance solid state 
oscillators comprise an active device such as Impatt, 
Transferred electron or Tunnel diodes, with a suitable 
BassivVe@ourcul ce.  Inethesschemetre  ciagram, Fig. el, 
the oscillator is represented by a waveguide cavity 
containing a negative conductance device. An adjustable 
short circuit plunger and a slide-screw tuner or other 
Suitable device optimizes the output power. In prac- 
tice, the cavity may take varrous forms), cylindmical or 
coaxial, for example. In the case where the output port 
is of the coaxial form, the slide screw tuner is most 
likely replaced by slug tuners. Coupling to the load 
may be accomplished through either a window or a probe. 

Because the dimensions of the semiconductor chip 
are generally very small compared to the free space wave 
length, at the fundamental frequency of oscillation, the 
conventional voltage and current values can be deter- 
mined at the terminals of the active device. Therefore, 
even though equivalent circuit components for various 


devices may have differing voltage and frequency 
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Fig. 
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Schematic diagram of a microwave 
GSscillator. 


Equivalent circuit representation 
of a microwave oscillator with 
separate resonant circuits 
representing each mode. 
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dependence, the general theory for negative conductance 
oscillators will be the same for any active device. 

If the harmonic content in the waveguide at the 
output port is small, at each harmonic frequency the 
device sees an impedance un-affected by the load con- 


dition. Under these circumstances, the active device 


admittance at the fundamental frequency exhibits a certain 


nonlinear behaviour necessary for it to achieve a stable 
operating condition. 

In general, the current and voltage wave forms 
generated by an active device are nonsinusoidal. Depend- 
ing on the embedding circuitry, it can be assumed that 
either the voltage or the current contains the basic 
frequency component only (voltage or current excitation 
of the active element) while the device response, i.e. 
current or voltage, contains a number of harmonics. By 
Fourier analysis of the oscillator wave forms or by 
measuring the large signal admittance, it is possible 
to obtain the large signal admittance of the active 


devices as defined by: 
or YA = Selb jJBL)A eos oy 


where A and I are the fundamental frequency components of 
the device voltage and current; - is the device's elec- 
tronic admittance. The electronic parameters of the 


device, i.e., Gy and Boe are nonlinear functions of the 
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oscillation frequency, Wor and the voltage amplitude, A, 
and abso depend On tie device bias and thermal condi— 
tions. On the other hand, the passive circuit admittance 
isa nonlinear Lunction OL trequency OnLy. 

Circuit models must take into account the inherent 
nonlinearity of the active device. To aid the under- 
standing of the operation of microwave oscillators, a 
very Simple oscillator model is first introduced and the 
behaviour of these oscillators under transient and steady 
state conditions is examined. Later in this chapter 
Kurokawa's geometric model of [17] microwave oscillators 
based on device line and impedance locus will be briefly 
presented to discuss the behaviour of negative conductance 


solid state oscillators. 


Baza toUuvalrentl Crreure OL a negative conductance micro 


wave oscillator 


Unlike their counterparts at lower frequencies, 
microwave resonant cavity parameters are defined based 
upon the choice of relevant field quantities. The stored 
Magnetic and electric energies in the electromagnetic 
field inside the cavity determine the equivalent induc- 
tance and capacitance, whereas ohmic losses in the 
resonator walls and dielectric losses determine the 
equivalent conductance. A cavity can resonate at an 


infinite number of frequencies. 
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A microwave cavity without excessive losses which 
is coupled by some arbitrary transducer to a waveguide 
is accurately represented by the equivalent circuit 
shown in Fig.: 3.2., The parallel reSonance) circuits, of 


resonant frequencies Wor Ypres eons represent the 


on 
resonant modes of the microwave cavity. The load admit- 
tance (Go + jB,) is coupled to each of the resonator modes 
by the mutual inductance M, of a coupling element having 
a self inductance L,- An active element within the cavity 


can excite these modes and is characterized by admittances 


(Capt jbo, (Cea. Jb) mene: electronic admittances 


voltage at the appropriate terminal pair; that is, 
4 = GredadB = (3.3.2) 


wheresmraH (O71 )2 , ste pls 

Resonators for oscillators are normally designed to 
possess fairly well-separated modes, permitting oscillation 
only in the vicinity of the desired frequency. According 
to the Foster Reactance [18,19] Theorem the behaviour of 
a resonator within a narrow range near the selected fre- 
quency may be represented by retaining the parts of the 
circuit representing the corresponding mode, Fig. 3.3. 


In this diagram: 
Se =“shunteconductance of the cavity 


L,C) = equivalent inductance and capacitance 
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L, = equivalent inductance of the coupling iris 
(window) 


M = mutual inductance between L. and Le 


1 
The diode package and bias circuit losses and reactances 
are lumped with those of the cavity since the diode is 


mounted inside the cavity. 


Pep OcClleatOle: UNucr sledGvestaue aia transtene state 


conditions 


3.96... Analysis offa highly simplified; model of a. two 


terminal oscibliatorveircuit 


In this section a simplified analysis for micro- 


wave oscillators is introduced. Although such a simplified 


model cannot represent the complete behaviour of micro- 
wave oscillators, the analysis helps to explain some of 
their features. 

The equivalent circuit of a two terminal microwave 
oscillator developed in the last section (Fig. 3.3) can 
be transformed into the simplified form shown in Fig.3.4a. 
The electronic susceptance, Bar may be considered to 
result from an equivalent inductance L Or capacitance 
Co: Therefore, Fig. 3.4a may be reduced to the simplified 


form in’ Fig. 3.4b94im Wheen 


(3.3a) 
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Fig. 3.0  EQuivalent.circust Of “a (microwave 
oscillator designed to operate ina 
Single frequency mode. 
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Pro. 3.4a° Equivalent circurvt of a. solid state 
microwave oscillator tuned to one 
mode Of Oscillation. 


Fig. 3.4b Highly simalified equivalent circuit 
of a microwave oscillator. 
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Since Be can be either inductive or capacitive, L, and 
Cc. are not present simultaneously in the equivalent 
CLreult. 

The electronic admittance, Yor is not constant 
with amplitude and frequency; hence the relation between 
the voltage across the active device and the current 
through it is not linear. Since one consequence of non- 
linearity of circuit elements is the generation of 
harmonics in current and voltage, 1t follows that. varia- 
tions of Ye with amplitude and frequency are always 
associated with the presence of harmonics. A strictly 
rigorous analysis would have to take harmonics into 
account or else show that the presence cf harmonics has 
negligible effect upon the steady state frequency and 
amplitude. Because the following treatment is intended 
to be of value principally in explaining qualitatively 
some of the observed characteristics of oscillator 
behaviour, harmonics of the driving force (voltage or 


current) will be neglected. It may also be noted that 


the dependence of Be upon the amplitude causes some change 


in the frequency as the amplitude of oscillation builds 
up. 

Under the rather severe assumption that all para- 
meters in the circuit of Fig. 3.4b are constant and 


linear, summation of the currents at one of the nodal 
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points yields the following differential equation for the 


voltage v across the elements: 


2 


inv G dv Vv 

—~tawtwa=0. (4) 

at Gidt CL 
TRERSOLUCTON- Ofahcr. (3.4) 1s 

2 Ce 

v= Vee Sin(wt + ¢) (35) 
in which 

@ Sena ear-iG,d GackG. 2G 

2c ve a e 
ss mice Ge) 120 C30) 
i) Sree 72M Tie Bhe (3.7) 


and the constants Ne and ~ are the amplitude and phase 
angle of the voltage at t = 0, respectively. 

PQ. 3.7) snows, thatthe circultiucan) be oscillatory only 
TEC os oe Examiiatdlonsor Eqs.) &(3.5)) and, (3:6) 
discloses that since the resonator and load conductances 
oar and G, are positive, sustained Sinusoidal oscillation 
is possible only if the electronic conductance Gs is 
negative and equal to or greater in magnitude than Gu: 
Le Ge is negative and equal in magnitude to Of, Cs 
zero and the oscillation is of constant amplitude. 
Theoretically, 2 Ge is negative and greater in magnitude 


thaneG, the anplatudesot cscitilation will inerease con- 
— 


tinuously. However, the magnitude of Ge is dependent on 
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the amplitude of oscillation. G, may remain constant or 
increase in magnitude for some amplitude range, but there 
is always an amplitude above which it decreases. Even- 
tually, an equilibrium amplitude is reached at which the 
dynamic value of Ge is equal in magnitude to Gu: What 
the value of the equilibrium amplitude may be depends 
upon the manner in which Ga varies with amplitude, the 
resonator losses, and the load. The time which oscilla- 
tions require to reach the steady state amplitude depends 
upon the behaviour of e with frequency as well as the 
effect of voltage amplitude and frequency on G.: 

Eqevs (3.5) indicates. that, since v is zero when We 
is zero, oscillation cannot take place unless initiated 
in some manner. In an actual case, the thermal motion 
of electrons sets up random pulses which will initiate 
Oscillation i£ the static valuewof Gx. at zero amplitude 


is medative.and |G ea Cie: 
e c 


a eoeee DULLaYD and limitation of amplitude of oscillation 


Eqs. (3.4)-(3.7) indicate whether or not sustained 
oscillation is possible; whether the amplitude at any 
instant is constant, increasing or decreasing; and the 
frequency at which oscillation takes place. They give 
no indication of the equilibrium amplitude of oscillation, 
the manner in which the amplitude builds up, or how the 


equilibrium amplitude varies with circuit parameters. 
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The equilibrium amplitude of oscillation can be 
determined only if the manner in which the electronic 
conductance G. varies with amplitude of oscillation is 
known. However, if a curve of Ge versus amplitude can 
be constructed, a simple graphical method may be used. 
The trend of such curves is shown in Fig. 3.5 along with 
a curve of Gy versus amplitude. Useful qualitative 
information concerning the behaviour of oscillators may 
be derived from the analysis of the curves in the figure 
even though it does not necessarily represent any parti- 


> G_, the exponent 


Cularvoscillater.@ If G.< 0 andi7G |e > 
e elo Cc 


Oeln Baw. (3.0) 1S 'zero! or positive.» Any small initial 
voltage WA resulting from random noise or other circuit 
disturbance initiates oscillation. When rent exceeds G_ 
the exponent a is initially positive and the amplitude 

of oscillation increases. As the amplitude increases, 
IG. | decreases, until the amplitude reaches the value 
corresponding to the intersection of the |G. and G, 
curves. Further increase of amplitude would cause |G. 
to become smaller than Gus Consequently, a would become 
negative. The amplitude would therefore decrease to the 
value corresponding to the intersection. This analysis 
indicates that the intersection of the curves corresponds 
LOpaw cond ton Of stable. cqultibrium in the circuit ‘and 
that sustained oscillation takes place at the amplitude 


corresponding to this point. 
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The curves in Fig. 3.5 show that the amplitude of 
oscillation can be increased by two methods. The value 
OF Go may be reduced by changing the load or the resonator 
16ss: or: the IG. | curve may be raised, which is ordinarily 
accomplished by changing one or more of the operating 
parameters of the electronic device which affect the 
negative electronic conductance. In some oscillators, 
however, the magnitude of the electron te conductance 
increases at low amplitudes, as indicated, for example, 
by @ne cumvesPotfrigs 3.5. “Again, oscillation cannot 
start unless ioe is equal to or greater than Go: FOr 
cutve (a) of Fig. 3.5, ous just exceeds G, and a is very 
small, but positive. As the amplitude increases, |G, | 
also increases, a grows, and the amplitude increases more 
rapidly. The action is cumulative and the amplitude 
increases rapidly until it reaches a value corresponding 
to the intersection of the curves at. pointoP. <tf the 
curve of the electronic conductance is then raised by 
changing some parameter of the source of electronic 
conductance, the amplitude increases continuously with 
change in the parameter. When the parameter is changed 
in the opposite direction, so as to lower the curve of 
electronic conductance, the amplitude decreases con- 
tinuously with change of parameter until the curve becomes 
tangent to the ee line, as shown by curve (b) of Fig. 


3.9. <Any further loweming of the electronic conductance 
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curves causes IG. | to be smaller than G, and a to be 
negative, whereupon oscillation ceases. It can then 
be started again only by raising the curve to a point 
at which Peele exceeds Gx and the amplitude again jumps 
te*that correspGuding*+tospoint P. 

If the electronic conductance curve has a maximum, 
as shown in all Gj-curves OG¢ EI Gte33 5 yjwand LGaihe is 
adjusted to a value that is somewhat smaller than Gur as in 
(GS) 10h .higegse5y, ERG.oscillationsis-not self-starting. 
However, if the circuit is driven by external means to 
an amplitude equal to or greater than that corresponding 
to intersection 1 of the conductance curves, the ampli- 
tude will increase abruptly to the value corresponding to 
intersection 2 and oscillation will continue when the 


driving force is removed. 


3.4 Kurokawa's simple theory of oscillators: free-running 


oscillators 


In practical oscillators, the device impedance is 
a function of both amplitude and frequency of the RF 
current through the device. A graphical approach may 
clarify the behaviour of oscillators even when a simpli- 
fied version of the device impedance has been used. 

The characteristics of negative conductance 
oscillators will be presented along the lines proposed 


by Kurokawa's simple theory of oscillators [13,17,20]. 
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Since the model proposed considers the equivalent circuit 
of the oscillator in series representation, the impedance 


concepts will be the most convenient to use. 


ay cae Ul Equivalent circuit 


In this theory, the active device is treated as a 
general impedance, -Z (A), with real and imaginary parts, 
R, (A) and X (A), which are arbitrary functions of the 


amplitude of the RF current, A, through the device, 


driving a general passive impedance Z,(w). The dependence 


of the device impedance on the frequency is neglected. 
This is justifiable since the device impedance is usually 
a slowly! varying function of w as compared to the circuit 
impedance. For the purposes of qualitative analysis this 
dependence of ~Z (A) on w will not affect the results of 
the discussion, but simplify the analysis considerably. 
The equivalent circuit for the proposed model ae shown 

in Pig. 8.6... The equation for free-running oscillation 


is given by 
= = ) 
[Z (w) Z (A) ] a Oi. (37330) 


This indicates that the total impedance times the current 
is equal to zero; namely, there is no applied voltage. 
Because of its simplicity, the oscillator model represen- 
ted by Eq. (3.8) reveals a number of important aspects 


of microwave solid state oscillators. Some of those 
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Equiavalent-circnit of. a 
free-running microwave 
oscillator. 
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Impedance locus and device line. 


Day, 


oY 


i 


Ss | 


related to free running-oscillators will be considered. 


3.4.2 Impedance locus, device line, and operating point 


To demonstrate the real meaning of Eq. (3.8), the 
locus of the circuit impedance and that of the negative 
of the device impedance are drawn on the complex plane by 
Varyand) © and A, as shown in Fig s,3.7..° The locus of 


Z() is called the impedance locus and that of Z (A) is 


called the device line. The arrows indicate the directions 


of increasing w and A. The intersection of these two loci 


correspond to 


Z,(w) = Z (A) A (3.3) 


When the RF current amplitude flowing through the device 
is A, the pornt. corresponding to Z (A) is called the 
Dperating pDOInNt.  ¢ eletsetinice, and Eq., (3.8)" 1s Satis 
fied, Z , (w) must be equal to Z (A). Consequently, in 
steady state, the operating point is located at the 
intersection between the impedance locus and the device 
line. The steady state oscillation frequency and ampli- 
tude can be determined at the intersection point of the 


EWO -LOCL. 


3.445 Prtects of cdistank retlection 


With this much understanding of the device line 
and impedance locus, typical behaviour of microwave 


oscillators like the one illustrated in Fig. 3.1 can be 
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discussed. First, the probable shape of the impedance 
locus must be deduced. This can be done conveniently 

on the Smith chart. The vector’) drawn from the center to 
a given impedance point on the Smith chart indicates the 
corresponding reflection at the reference plane where 
the impedance is measured, 

Reflections in an oscidlator Circuit,éexhibit a 
phase rotation with a change in frequency at a rate 
determined by the electrical distance of the reflector 
from the active device. The longer the distance the 
quicker the rotation. Physical structures near the 
device, such as the short circuit plunger or the wave- 
guide post used to bias the device, create large reflec- 
tions which rotate slowly. On the other hand, the slide 
screw tuner which is located at a distance from the 
device gives a reflection which rotates quickly. Because 
Of thiseaquick sotatien, the tip of the resultant vector 
most probably draws a loop or loops on the Smith. chart 
aseadmustrated inshig.93.9. Consequently, itt is most 
likely that the impedance locus has loops. It may be 
worth mentioning that the distant reflection has less 
tendency to form loops when the main reflection vector 
rotates quickly. One way to avoid loops in the impedance 
locus is to make the Q of the close-in circuit as high as 
other conditions permit and the distant reflection as 


small as possible. 
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If adjustments are made in some suitable circuit 
component such as the short circuit plunger, the posi- 
tion of the impedance locus changes with a loop down- 
Wards as shown,in Fig. 3.9a. ~The operating point Py 
first moves to the right as indicated by P. but as soon 
as the upper edge of the loop separates from the device 
line, the operating point jumps to the left. Eigen Seis 
will help understand the detailed behaviour near the 
upper edge of the loop. As the operating point Ps , 
approaches the upper edge of the loop, the intersecting 
angle becomes small. This indicates that the oscillation 
becomes noisy. If the locus is moved to the position 
indicated by the dotted line, the intersection disappears. 
No steady state oscillation is possible near the upper 
edge of the loop. Consequently, the operating point jumps 
to the other intersection Poe If the maximum power point 
is located at the position indicated by X, the power 
initially increases as the locus moves downward, since 
the operating point approaches the maximum power point. 


But when the operating point jumps to P the power 


bu! 
suddenly decreases and the frequency increases as illus- 
trated in Fig. 3.10@ee2£tthe adjustment is reversed’, the 
locus moves upward, the operating point moves to the 
right, and the frequency as well as the power follow a 


different path; and when the lower edge of the loop 


separates from the device line, the frequency and power 
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operating point. 
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suddenly jump again. In this particular case, both fre- 
quency.and power decrease... This behaviour is typical of 
microwave solid state oscillators having distant reflec- 
tions and in accord with experience. In practice, 
however, the impedance locus may contain many loops and 
the loops move in various directions while changing their 
Sizes. The result is a more complicated and frustrating 
oscillator behaviour. 

Similar complications will arise when the bias 
supply is adjusted while keeping the rest of the circuit 
unchanged. In this case, the device line moves with 
respect to a stationary impedance locus resulting in an 


almost identical effect. 


So424 VOnset For oscillation: transient and steady state 


conditions 


In this section, the steady and Serene state 
conditions of a microwave oscillator will be discussed 
along the lines proposed by Kurokawa [17]. 

Because of the white-band noise in the lossy 
components, there will be present in the embedding cir- 
cuitry small RF currents of proper frequency after the 
bias supply is turned on. Under these conditions, the 
equation of the free-running oscillation is approximated 
by 


[4 (w) = Z,(O0) 1 = 0 (35 .0)) 
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where =2.-(0,) is the small-signal device impedance. Since 
I is small but nonzero because of the switch-on distur- 
bances or noise, and since Z(w) is not equal to Z (0) Lo 
any real w, w must be complex. The real part will be 
expressed by w' and the imaginary part by -a. Since the 


Was te avs 
aie aaks positive 


time factor corresponding to w'-ja is e 
a indicates the amplitude is increasing with time and 


negative a indicates a decreasing amplitude. A concrete 


example in which 2 (w) represents a series resonant circuit 


will help make the point clear. 4,(w) is given by 


G3 s4-19) 


ee _ 1 
Zt) = j(wL De + Ry ton Ry 


where Ro represents the loss in the cavity and Ry the 


LOa@ erin the vicinrcy oL resonance, 4.,(w) is approximated 


by 
Z(w) os 2j(w-wi)L + R, + Ry, Coeck2) 
=) 24(ml — ) 
23 (w wi L + 20L + Ri+ RB, Gaede) 
where 
Qj. = ake : (3.5145 
Sn rane 


The impedance locus, the locus of Zw) with real w, is 
a straight line parallel to the imaginary axis. The a- 
constant loci are all parallel to the impedance locus 
and} the w'-constant loci are perpendicular to it. The 


right-hand side of the impedance locus corresponds to 
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Typical frequency and power 
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multiple-tuned oscillator. 
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positive a's and the left-hand side to negative a's as 
inarocated “Ln ir ld.s osu. “Ime thverciagram Z (0) corresponds 
cova positive ‘a. Consequently sat Eq: “S00 )*hdotds ,—d 
must be positive whenever I # 0. In other words, the 
circuit device interaction forces the current amplitude 
to increase. As the current amplitude increases, the 
operating point moves left along the device line toward 
the intersection with the impedance locus. The figure 
indicates that the frequency as well as the amplitude, 

in general, changes as the operating point moves since 
the value of w' changes along the device line. Since a 
is positive and the amplitude increases as long as the 
operating point is located on the right-hand side of the 
impedance locus, amplitude growth continues until the 
Operating point reaches the intersection satisfying 

Z,(w,) - Z(AQ) = 0, or equivalently 
2L(w - w.) +X (AQ) = 0, eases ~RSA mes aa) C3 oe) 


where WS and AG are the frequency and amplitude corres- 
ponding to the intersection and X(A) is the imaginary 


DALLLOL ~Z (A) . 


S740 Nvolobr lit veOteerecentunning oscillators 


If the amplitude increases beyond the intersection, 


the circuit-device interaction forces the current to 
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decrease since the left-hand side of the impedance locus 
corresponds to negative a's. Consequently, the operating 
point returns to the intersection. In other words, the 
snLersectron. Ine hig..65 ub) Peeuite in a stable steady 
state operating point. 

In case if the impedance locus has a loop, there 
could be more than one intersection between the impedance 
locus and the device line. The same arguments developed 


above can be used to differentiate between stable and 


unstable operating points. If two stable operating points 


happen bo occur, the starting history of oselllation will 


determine which operating point is functional. 
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CHAPTER IV 


IMPATT DIODES 


lh Introduction 


Impatt (IMPact ionization Avalanche Transit Time) 
diodes are two-terminal devices utilizing the combined 
effects of avalanche breakdown and transit time to 
Cransiorm parc or ithe de energy into RF energy ina 
certain frequency range (from the circuit point of view, 
the real part of the diode impedance is negative in this 
frequency range). A p=-n junction can be operated in its 
Impatt mode when it is biased into reverse avalanche 
breakdown and mounted in a proper microwave cavity. At 
the present time, the Impatt diode is one of the most 
powerful solid state sources of microwave power. A CW 
efficiency of as high as 35.6% (with 2.9 W output power) 
at 10.4 GHz has been reported [21]. A hi-lo Schottky-barrier 
GaAs Impatt diode has given a CW power of 12 W at 6 GHz 
oe 

The possibility of producing high frequency oscil- 
lations from a reverse-biased p-n junction device operat- 
ing under avalanche breakdown conditions was first 
discussed in 1958 by Read [22], who first proposed the 
Read structure (in p-n-i-n configuration) which has a 


peculiar doping profile. Later, Missawa showed that a 
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negative conductance of Impatt nature can be presented by 

any p-n junction diode with almost any doping profile [23]. 
In thisrehapter /-the-devicesstructure: and its 

dc characteristics are considered first. Section 4.3 

describes the basic principles of operation of the Read 

diode. Small-and large-signal analyses are discussed in 


section 4.4, 


422°" Device structure: ana dC Characteristics 


A p-n junction is the simplest Impatt structure. 
When reverse biased, mobile charges from the regions 
close to the metallurgical junction are swept away by the 
field, leaving a positively charged (n-region) and a 
negatively charged (p-region) depletion region. At the 
same time, a small reverse saturation current flows. As 
the bias voltage increases, the depletion layer becomes 
thicker and the unswept portion decreases. These 


i " 


unswept portions of the "p" and 


" " 


n" sides introduce a 
resistance ‘which degrades diode efficiency [24]. When 
ere as Vy, (V) is the breakdown voltage), the introduced 
depletion zone acts as a nonlinear capacitor. Fig. 4.1 
shows the equivalent circuit for a p-n junction diode 
at a given bias voltage V (Vi<V)). tn—this- figure Rp (V) 
and C(V) depend on the semiconductor characteristics. 
If the reverse bias voltage is increased to the 


value of V the junction breaks down so that a small 
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Fag. 4.45 Impuc diode equivalence circuit at VC V,, where 
V = bias voltage and V, = breakdown voltage. 
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Fig. 4.2 Read diode 
{a} doping profile. (ptnant), 
(hy _ welect ramet 1610 sdistribution,.and 
(c) ionization integrand at breakdown. 
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increase of the bias voltage beyond this critical value 
results ina large current flow. This large current is 
the result of an impact ionization producing multiplica- 
tion of electrons and holes in the depletion region of 
the junction. This phenomenon is known as avalanche 
breakdown. 

Impatt diodes are designed in various structures. 
The Read diode, the one-sided abrupt p-n junction, the 
two-sided abrupt p-n junction, the p-i~-n diode and the 
lo-hi-lo structure are the most commonly used structures 
fest. The Read diode iS a convenient structure for 
analytical purposes and for obtaining an understanding 
of the physics of the Impatt diode operation, but it is 
difficult to manufacture. One and two sided abrupt p-n 
JunGtions are easier to fabricate but very difficult to 
analyze. Impatt mode oscillations are obtained from 
these more practical diode structures in much the same 
manner as from the Read diode. The most widely used 
Impatt diode structure has been the one-sided abrupt 
pen-n’ junction (p-mn is used instead of pn in order 
to reduce the diode parasitic resistance) which will be 
used in the experimental part of this study. 

Fig. 4.2 shows the doping profile, the electric 


field distribution, and the ionization integrand at 


, a ; ; : 2 + oo 
breakdown condition for an idealized Read diode (p -n-i-n ) 


: Y Pha porate : : At ‘ 
or its dual (n~p-i-p ). The ionization integrand is 
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where a and om are the ionization rates of electrons and 


a 


holes, respectively, and W is the depletion region width. 


The avalanche breakdown condition is given by 


WwW 
| <qedx = 1 . (4°°2)) 
0 


The avalanche region is highly localized. Therefore, 
because of the strong dependence of a's on the electric 
fmveld, most OLgthe multiplication process occurs: inla 
narrow region near the highest field, between 0 < x < Xn 
where Xn is defined as the width of the avalanche region. 

The voltage drop across the avalanche region is 
called Vas Both Xn and eS have profound effects on the 
optimum current density and the maximum efficiency of an 
Impatt diode [25]. 

The corresponding results (the doping profile, 
the mee hie Oss distribution and the ionization 
integrand) for a typical one-sided abrupt p-n junction 
are shown in Fig. 4.3. Again, there is a noticeable, 
highly localized avalanche region. Fig. 4.4 shows the 
results for a p-i-n diode which has a uniform field across 


the intrinsic layer. The avalanche region in this case 


corresponds to the full intrinsic layer width. The 
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4,3 One-sided abrupt p-n junction 


(a) 
(b) 
(ec) 


doping. profile (pin), 
electric field, and 
ionization integrand at 
breakdown. 
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Fig. 4-4 p-i-n diode 
Cay doping prorrie”(p=1-n); 
(b) electric field, and 
(c) ionization integrand at break- 
down. 
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breakdown voltage, the maximum electric field and the 
depletion layer width can be calculated from the break- 
down condition and the field dependence of the ioniza- 
tion rates. Schroeder and Haddad [26] carried out some 
computations on several diode structures. According to 
their calculations an n-p structure has a higher break- 
down voltage with respect to the on case of the same 
doping concentrations. The breakdown voltage for the 
asymmetrical two-sided Si and GaAs diodes is higher than 
that of the one-sided diodes of the same doping level. 
In general terms, the breakdown voltage depends mainly 
on the doping profile, the semiconductor material and 
the width of the avalanche and depletion regions. 

For a Simple structure as n=p or p'-n, the doping 
profile can be accurately deduced from the measured 
voltage dependence of the junction capacitance. For more 
complex structures, there is no simple procedure for 
evaluating the doping profile. Therefore, only in the 
case where the internal structure of the diode, its 
junction geometry, and its doping profile are known, can 
the Impatt -govameunel characteristics be evaluated. Unfor- 
tunately, manufacturers do not supply any information 


concerning the detailed structure of the diode. 
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#1) basic principles of impatt, operation 


A good understanding of the ac characteristics 
of Impatt diodes can best be obtained by considering 
the operation of a Read diode [22,27,28]. The 
structure, field distribution, and current waveforms 
for a reverse-biased Read diode are shown in Fig. 4.5. 
The Read diode consists of two regions: a narrow 
avalanche region (p-region) in which carrier multipli- 
Gatleneby, impact j10n1zation.occuss, sand .a.~drift region 
(intrinsic region) sin ~which.the+carniers -drift.at 
Saturated or field-independent velocities. No impact 
1.Gnt790ELOR,OCGUES wlan thesdmitt @eregion,. 

The negative conductance or resistance of a Read 
diode is attributed to a phase shift between the current 
through the diode and the voltage across it. This phase 
shift consists of two components. There is a phase 
delay in the external circuit current caused by the 
avalanche multiplication process and a delay caused by 
the finite transit time of the charge carriers drifting 
through the dEift.region 

The phase delay caused by the avalanche process 
results from the fact that the rate of generation of 
electron-hole pairs in the avalanche region is propor- 
tional both to the electric field and to the density of 
electrons and holes already present in the region. If 


a small RF voltage of sufficiently high frequency is 
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assumed to be applied along with a dc voltage near 
breakdown across the diode, the rate of generation of 
electron-hole pairs will exceed the rate at which the 
pairs leave the avalanche region. Thus, the carrier 
density in the avalanche region will grow exponentially 
with time as long as the RF and dc voltages add to give 
a total field above breakdown. When the RF voltage 
changes sign and subtracts from the de voltage, causing 
the field to fall below breakdown, the avalanche multi- 
plication process drops almost to zero, causing the 
current to decrease. Thus the avalanche current will 
have its maximum when the RF field goes through zero. 
Thus, as shown in Fig. 4.5, the avalanche process con- 
tributes a..90? inductive phase lag to the current 
generated in the avalanche region. This current is then 
injected into the drift region of the diode. The current 
anauced in) thelexternal circuit by 4enis charge 7 drifting 
threugh the (drift: reqron at -avsaturated, drift velocity, 
Ho@ehownaat the DOtEOMuot Fig.) 4.5. ),.Lt Ls obvious that 
the fundamental component of the external current is 
more than 90° out of phase with the RF voltage, causing 
the resistance or the conductance of the diode to be 
negative. 

Thus a net transfer of energy into the ac field 
from the dc field has occurred, and the ac voltage can 
grow. In an appropriate circuit, self-sustained oscil- 


lation Can be obtained: 
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Fig. 4.5 Read diode structure, field 
distribution, and voltage and 
current waveforms. 
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4.4 


Small-and large-signal operation 


The equations which govern the operation of an 


Impatt diode are Poisson's equation and the continuity 


equations for holes and electrons. 


these equations are 


and 


where 


where 


an 
ot 


M |4Q 


In one dimension, 


(N) - Natp-n) (43) 
og 
D 
i gq ox 
(ama) 
' og 
bhai baat ia 


electric field (V/cm) 


electronic charge 


(C) 


dielectric constant (F/cm) 


donor and acceptor densities 
hole and 
hole and 


hole and 


rates 


hole and 


(cm™ >) 


electron densities (cm 


2 


electron current densities (A/cm) 


electron generation recombination 


electron mobilities (om*/v.s) 
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k Boltzmann constant (1.38 yes ea Joule/°K) 


i absolute temperature (°K) 


Usually, several simplifying assumptions are made 
En Order (to obtain a. closed form solution for the diode 
small-signal admittances. In many small-signal studies, 
ionization coefficients of electrons and holes are 
assumed to be equal, and recombination effects together 
with the diffusion current are neglected. Even with 
these assumptions, closed form solutions are obtainable 
only for simple structures such as the Read and the p-i-n 
structures. For more complicated structures, solutions 
are usually obtained by the use of the digital computers. 

Small-signal analysis of a Read diode gives the 
following expression for the real part of the diode 


terminal impedance [29] 


R = W r 2 - cos 6 | [a7 6) 


where W is the depletion region width, A the diode area, 

@ the transit angle given by uW/v,, (vi, is the scattering 
limited velocity of charge carriers) and Wy the avalanche 
resonant frequency given by 


Bag 
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DO 


in which the quantity a' is the derivative of the loniza- 
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tion rate with respect to the electric field and Je is the 
dc current density. From Eq. (4.6), the diode negative 
resistance can only be negative if w > Ws 

The variation of the negative resistance with 
the transit angle is plotted in Fig. 4.6. -The peak value 
occurs near 0=7. For transit angles larger than tm and 
approaching 37/2, the negative resistance of the device 
decreases rapidly. 

The negative conductance at a fixed operation 
frequency as a function of the biasing dc current density 
is shown in Fig. 4.7(a) for the Read diode (p'-n-i-n), the 
abrupt pen junction ,\and the p-i-nvdiode. This figure 
shows that the negative conductance for a Read diode 
rises to a maximum value at a relatively low current 
density, then decreases, and eventually disappears at 
current density at which the avalanche resonant frequency 
(which is proportional to Was Ref. [30]) equals*the 
operating frequency. The other structures show very 
Similar characteristics. However, the conductance at 
low current density is smaller and the peak negative 
Ponductae ae «docure at larger current for the abrupt 
junction. Even larger current density is needed for the 
Misawa (p-i-n) diode. 

At a fixed current density, the negative conduc- 
tance versus operating frequency for the above structures 


is shown [30] in Fig. 4.7). The Read diode has a relatively 
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Read diode: 
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Effect of doping profile on conductance for 

Ly Read, \2)) abrupt, 3) sMissawa (uniform 

aie profiles. (Reproduced from Ref. 

30) 1). 

(a) Negative conductance vs. dc current at 
one frequency. 

(b) Negative conductance vs. frequency at 
constant dc current density. 


Idealized negative resistance vs. 
transit angle. 
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narrow peak around the transit time frequency. A wider 
frequency range is obtainable from the abrupt p-n junc- 
tion with some reduction of peak conductance value due 
to nonoptimum phase relation. The reason that the 
negative conductance disappears at low frequencies (smail 
6) is that space charge in the drift space can exhibit a 
large positive resistance at low frequencies [31]. 
The Misawa diode, however, exhibits a small negative 
conductance, nearly constant at frequencies from zero 
to. 8 = 27. 

The foregoing discussion has considered the small- 
Signal conductance of several diode structures to be a 
function of the dc bias current and frequency. However, 
the dependence of the diode susceptance on the bias 
current and operation frequency is also important. The 
sensitivity of the small-signal diode admittance to the 
bias current and frequency, as given by C.A. Bracket 
[32] (for n’-p-p” structure), is shown in Fig. 4.8. 
Ler kee 4.8, the admittance characteristic is seen to 
be very similar to that of the Read dicde [22]. Gen- 
erally over an appreciable frequency range, the diode 
susceptance decreases and the diode conductance increases 
with the increase in the de bias voltage, 


At large current densities the simplifying assump- 


tions used in the small signal theories are not justified. 


Small-signal results are useful only at the initial 
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Fig. 4.8 Small signal admittance plane plot for Ge 
ntp-pt+ diode from Gummel-Scharfetter program. 
(Reproduced from Ref. [32]). 
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stages of oscillator design in order to indicate the 
oscillation frequency and the maximum negative conduc- 
tance obtained under certain bias conditions. Under 
normal operating conditions ,/e2q.84n an oscillator 
Circuit, the RF voltage is comparable to the bias 
voltage and the diode RF characteristics become ampli- 
tude dependent. 

Under large-signal conditions the ionization rates 
for holies and electrons depend strongly on the electric 
field. This causes the multiplication current to respond 
nonlinearly to the large voltage amplitude across the 
avalanche zone [83] ." Ini this™casey the’ iconduction current 
will have a high harmonic content. 

Several large-signal analyses have been undertaken 
which provide insight into the nonlinear properties of 
the Impatt diode [22,28,34, 8]. The most accurate and 
complete of these, a numerical solution, was given 2m L369 
by Scharfetter and Gummel [28]. Basically, the program 
they used performs a simultaneous time-dependent solution 
of Poisson's equation, the current equations, and the 
continuity equations for a one dimensional structure. 
Included in the program are accurate expressions for 
electron and hole mobilities, ionization rates, and 
thermal generation and recombination. Although these 
studies probably give a most accurate and detailed theore- 


tical description of the large-signal operating conditions, 


teins 248 | 
sede thine ne aya is reds 
gekd oat ios otdasato 42 i 6B 


~ i as ancoed aoisend ORS aia a 
. ’ oy ae j 

1 ita a a ee nae ¥ vs 4 Ft 
hs alll Sak Sent, ie ; y 


gota agtieriadk: ‘ats anott thee: Sue 
Lier ee; J 
jason lS wer A gipnor3e" 


ue i a Ae 7 
divoyaers OF InoaRHe voi nea sagtietin oe or 
aay aaoNor atid L.Cqpte spgstowiopnal 


hes Me oti teh. 100 dg (wees gna a phen 


+0 


it ae 
BGT ot ame ao Aguas dete uP 
+ ? | fi he - 
aotarrebay Aged! wpa sacl 08 Tiegh : 
v4 te. a | 
16 ‘aekstaqoxe ‘santh aoa “bt bine di oF, 


bia actaitnrcn, phon eit” Ube Says \, 

eae arp | 

goer wt avy ngewe tawihet Poolhronits Saaadd ie setan tf 
| Leas ae 

amapond wit eae bon tee) se el “aden aad ss 


+ 
7 


ame es aw le hs Be 


10i4u fou sophisti epaaared Ligeta” B npriotz9g boaw edt} 


ae perme? gow se i 


eh 


au 
a iy ” 


v, 
7 


Js 200 A/cm 


120 
100 
80 
UO 
= 
KG 
=a 
ay 
ea} 
60 oO 
op 
= 
ap) 
40 
© VOLTS AC 20 
~. 5.4 
L Jer, oe ; 
=2% ~20 ~15 ~10 -5 Q 
CONDUCTANCE 


Fig. 4.9 Diode admittance (susceptance vs. 
conductance) as a function of frequency 
and RF voltage amplitude, with the 
resultant efficigney indicated. Current 

. density 200 A/cm”. (After Gummel and 
Scharfetter |28)). 
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the large computational time necessary generally pre- 
cludes detailed studies of the operating characteristics 
under a wide range of parameter variations. 

Fig. 4.9 is a compact presentation of the impor- 
tant results of their calculations [28] for a Read 
structure. In most cases, the diode susceptance increases 
with increasing RF voltage amplitude, whereas the diode 
negative conductance decreases. The results show diode 
efficiencies as high as 18%, although neither oscillator 
power Mor efficiency is saturated at that point., It is 
difficult, with that program, to determine the maximum 
efficiency or power output because these quantities depend 
on the proper termination of the harmonics of the funda- 
mental frequency; these constraints are not known a priori. 

Since 1969, most of the effort in Impatt large 
signal studies have been directed toward reducing the 
computational time required for obtaining a solution 
[8,34], without causing a significant reduction in 


accuracy. 
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CHAPTER V 


MEASUREMENT OF LARGE SIGNAL ADMITTANCE OF IMPATT 


DIODES AND ASSOCIATED PASSIVE CIRCUITS FOR OSCILLATORS 


SSL *Iineroduck ion 


The theoretical models of the active and passive 
Circuits developed earlier can be implemented in 
Dractical applications, such as,;oscililators, only if 
the active device, e.g., Impatt diode, is accurately 
characterized. Accordingly a technique for the measure- 
ment of large signal admittance and the associated 
passive circuitry is described. This technique, together 
with the theoretical models, will provide an accurate 
method for dynamic studies on Impatt diode active circuits. 
Several authors [9,28,35-38] have published methods 
and/or results of microwave impedance measurements on 
packaged diodes. However, in order to accurately evaluate 
the active wafer parameters, it is necessary to measure 
precisely the parasitics and impedance transformation 
properties of the diode package, the leads, and the mea- 
suring circuit. So far, this has been achieved by the 
use of two techniques. In one technique these elements 
are represented by a 2-port network connecting the mea- 
suring system to the diode wafer, then the parameters of 
the 2-port network can be determined by standard measure- 


ments [39,40]. This technique, however, requires 
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replacing the diode wafer by three pei impedances 

for calibration purposes. This necessitates the use 

of identical dummy packages for the three reference 
impedances. To accurately substitute for the diode under 
consideration, the diode packages must be fabricated in 
the same manufacturing batch, which is commercially 
unpracticable. Furthermore, replacing the diode by these 
packages for calibration, one may change the actual RF 
conditions as seen by the diode, since the reference 
Planeemay; sitre Sitighntiy. “Accordingly, it 1s dirtficuit 
to attain a reasonable measurement accuracy. The second 
technique used relies on the capacitance vs voltage (C-V) 
characteristic below breakdown of the Impatt diode, 
which can be measured accurately at 1 MHz. Gewartowski 
Strat (oo) sand 2cCLO eb al. 132) Used cits cechnigue to 
obtain the parameters of the 2-port network connecting 
the measuring system to the diode wafer. However, the 
parameters characterizing their 2-port network are not 
uniquely determined. 

Even though a considerable amount of work has been 
carried out-in improving the microwave measurement 
techniques for Impatt diodes, the measurement accuracy 
is still very limited [41]. The magnitude of the packaged 
diode's negative resistance is of the order of 10 ohms 
Orvtece, ror vraclical Giodes. Under Conditions of low 


dc currents, high RF voltages, or maximum power outputs, 
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this magnitude can be as low as 3 ohms. Van Iperen 
and Tjassens [41] claim that it is questionable whether - 
the standard VSWR (slot tsa dane technique) or reflection 
coefficient measurement techniques are sensitive enough 
to produce highly accurate results at these low impe- 
dances. They advocate the use of a microwave impedance 
bridge method that requires a very stable frequency source 
and critically machined components. However, their method 
is very time consuming, especially when the measurements 
arecmadenatbyseveralhfirequenceresc 

Underalargessignalinconditions pmbhe faccuracymot 
the measurements is further restricted because the Impatt 
diode is a nonlinear device and will produce harmonics of 
the fundamental frequency during the measurements. How- 
ever, if the diode is terminated in a 502 line, the 
impedance at various harmonics is also 502 and the 
harmonic content is well defined. The disadvantage of 
this arrangement is that a large amount of input RF power 
is needed to drive the diode into large-signal conditions. 
Diodes in matched resonator mounts of reflection type 
amplifiers can ie carefully designed to reduce the con- 
tribution of the second harmonic, while reasonable power 
levels can be used in the measurements. 

The technique developed here for the measurement 
of large-signal admittance of Impatt diodes can be sum- 


marized as follows. The measurement system, the diode 
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parasitics and diode package up to a defined imaginary 
plane intersecting the active wafer are characterized 

by two 2-port networks. By knowing the parameters of these 
networks, the admittance of the active wafer can be 
directly determined at various frequencies and for 
different RF power levels. The error parameters asso- 
Ciated with the first network are evaluated by measuring 
the response of the measuring system to different known 
impedances corresponding to various positions of a 
precision movable short. Replacing the movable short 
with the Impatt diode mounted in an amplifier assembly, 
the parameters of the second network can be evaluated by 
measuring the impedance of the Impatt diode at various 
dc voltages below breakdown. In the present system, the 
Signal frequency as well as the RF power level can be 
easily adjusted without physically disturbing the diode. 
The measurement system has a 50 dB dynamic range and is 
easy to use for measurements at different frequencies. 
The active diode and passive circuit admittances are 


measured using the HP 8410A network analyzer. 


Dee Experimental ayrangement for admittance measurements 


Fig. 581 ista blocki diagram which illustratesmthe 
set up for the admittance measurements on Impatt diodes 
under both small-and large-signal conditions. One of the 


two isolators is used to isolate the TWT amplifier from 
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the klystron source and the other is used to isolate the 
[Wl amplitier from the rest of the system. The TWTr 
amplifier has a maximum gain of 30 dB. 

The diode functions as the active element of a 
reflection type amplifier. The schematic diagram of this 
amplifier is shown in Fig.5.2. It consists of a waveguide 
circulator, an impedance transformer with a movable slug 
for tuning purposes and an Impatt diode. A coaxial 
transformer matches the 50 % characteristic impedance of 
the transmission line to that of the diode chip. This 
transformer has a movable slug which allows measurements 
COmbe carriedwout anpa frequency; range iuom 5.7 «to-6.6 GHZ 
without disturbing the measuring system. The waveguide 
circulator is used to separate the waves incident on and 
reflected from the diode. 

The principle of operation of the measuring system 
is simple and can be summarized as follows. The output 
power of the TWT amplifier is split into two signals by 
Means of the directional coupler #1. The signal in the 
secondary arm is attenuated by the 13 dB attenuator to 
the reference power level of the network analyzer. In 
the main arm of directional coupler #1, the second signal 
is attenuated, by attenuator #1, directed through the 
cinculator fang\yGinallyaganected into the diode. "| At a 
specific power level and frequency, a reading is taken 


first switching the RF energy into the diode, and then 
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switching it into an RF short which serves as an accurate 
reference impedance. This procedure is repeated at each 
frequency and power level. Therefore, the data is 
referenced to the accurately known reflection coefficient 
of the RF short to minimize the effects of the slow drift 
of the network analyzer. The reflected signal from the 
diode is compared with the reference signal within the 
network analyzer. 

The characteristics of the microwave circuitry are 
governed by the circulator, its bandwidth and isolation 
being the most important parameters. As the diode is 
mounted on a precision 50 ohm connector, it has excellent 
matching characteristics (VSWR< 1.05) to frequencies 
beyond 18 GHz. However, the four=port circulator has’a 
narrow bandwidth (as Y-band components are used) and 
outside this frequency range it presents a highly reactive 
impedance. Thus, any second harmonic components generated 
by the Impatt diode would be largely reflected from the 
circulator and confined within the diode mount. In this 
Situation, the RF voltage incident upon the Impatt diode 
will have some second harmonic content present. Another 
TMpOYrtant aspect OL the Characteristics” of the four-port 
circulator is the isolation between ports 1 and 3 (see 
EVO. Del). Lnis “1SGlatcion describpesythe degree of “attenu- 
ation experienced by the incident RF in reaching the path 


of the reflected signal.” ~The isolation of the four-port 


35 


ee | , st we iM ey a o Ew, 


Bie Ne 1H, ee came 
Howe te bostsager ek 
ek a¢sb edt sone 


i 7" 


\< 
319 POdRpeoR fetaoel® sz aweonex 


ARivh wols wns « aves Te. odd ost 
ae mead tetp te ‘been teen eat’ 


ote neetetw Laaeee sonexaney ot 


- ii 


wisiverta evawor> ta ats 20 sshawbapingsato ot 


a 


3 P- ; en = gts ixbiwhned avs ctemttropt a9 vst & 


‘ie 


: r, 
ei eboll efit @A yeaadomeraq snstuogqns geon one eat 


eo 


swelionks aod ai .sotowmos amo Ce noLarossg 6 ao per 
. 7 _ 


sebpasdpes? of (2 Spat,  gHBY” opida Medsesetts, 


= 6° get godelus1i3 Babee? eee yas dehiedaad 2 Sr. 


’ 


aca tebaw wre ates roe) tavad~\y en) “tho Leteted wciigad «| 
ipen videiei s as tian sang Sk spies < ondlipent aldy ‘obtedao 7 


betd ra (vei TY pig oo 2 rescmngee bienee yous ets .omshoqmk iv 
(Y mors’ bola l bet hove? 62 biuaw oe) ‘stoqnt ott yal. : at 


oe ee, 


gies at. ahem SBD iM sittsh mbisie Iona nat bal soteiuakte yy 
abo kit aes eye ot sitledionst opatlew i oa pen . 
relizoi  Saesag siege aipomre brompe ate 
4 yen apie ate ho sip advo pooenacts eit bayh uini si r 
yaa) heh aoadg reused, Lenaend re ° Saas 
staat Te cagcal bh aoe eB: pam ek ote a Bs a pat | 
fees eis sei avast nl Wi, saw btontt ita . Lf 
epost Ge Sota ath” anny lene 


ai 4 
hee 
“ Cw he he 


eee ‘ 
at Te 


96 


circulator used here was greater than 30 dB between ports 
E©and 3. 

Measured data consist essentially of two sets of 
dc voltages. One is proportional to the return gain, 
RG, and the other to the phase angle of the reflection 
coefficient, 6, such that ¢(degrees) = Vac (mv) /10. The 
return gain is defined as negative return loss in 
decibels and is converted to the reflection coefficient 
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All measurement systems introduce systematic errors 
into the measured results. In precision measurements, 
these errors must be greatly reduced by either designing 
a system very close to the ideal one or else by using a 
system which can be calibrated in situ by substituting 
reference standards for the device under test. At micro- 
wave frequencies, the former approach has usually proven 
to be impractical. This is because high frequency signals 
must be measured using such devices as directional 
couplers, circulators, switches, and other components 
which cannot be designed with ideal properties. The end 
result is a large measurement error. Furthermore, in 


using network analyzers for impedance measurements, the 
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slow drift in the low frequency circuitry creates a 
serious problem, especially when making point-by-point 
measurements of small impedance changes. 

Fortunately, at microwave frequencies, one can 
use impedance standards such as movable shorts. In 
this case, the response of the system to different 
positions of a movable short may be used to evaluate 
a set of error parameters characterizing the error 
network shown in Fig. 5.3. These error parameters are 
utilized to correct the measured results. 

The error network representing the measurement 
system between planes CC' and AA' can be characterized 
by replacing the Impatt amplifier and its bias components 
at plane AA' (see Fig. 5.1) by a precision movable short 
as mentioned earlier. The error network of the measure- 
ment system has to be measured at each frequency and 
power level by recording the return gain and phase for 
different short positions. The reflection coefficient 


if due-torthe RPeshorus(Pigy 5.4), 1s transformed into 


AA! 


an input reflection,.coefficient et through-the error 
network. For a known error network, the functional 


dependence between Ne and Daal has the form of a bilinear 


transformation [43] 


ANT a1 + BN 
rs + a, 
m CNI,art i 


where Pant is the reflection coefficient at plane AA' due 
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Fig. 5.3 Measurement system representation for 
evaluating the Impatt diode admittance. 


Fig. 5.4 Measurement system representation for 
evaluating error network parameters. 
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COisthe yliding shortrcincusit tat jayidiistance 2.5 ib- iis 


given by 


aie tan eerep lume slse, 


AA! (583)) 


where 2 is the plunger position with respect to the 
reference plane AA', 
a is the attenuation coefficient of the transmission 
rire, 
and 8B is the phase coefficient of the transmission line. 
However, available precision shorts behave almost ideally 
and therefore Eq. (5.3) may be simplified to 


-j472/2r 
tae 


Ny =/—e 


AA! (a = 0) (5.54) 


where eG is the guided wavelength. 

There will be m equations of the type (5.2) 
corresponding to m positions of the short. Since each 
of the quantities AN, BN and CN are complex numbers, 
there are six unknowns to be determined. This could be 
accomplished by solving six independent equations derived 
from three separate measurements. Therefore, if a large 
number of measurements are taken to determine AN, BN and 
CN, systematic errors could be minimized by utilizing 
numerical optimization techniques. This can be achieved 
by considering three measurements and thus obtaining 
three equations from Eq.(5.2) and then using a numerical 


routine for arriving at the values of AN, BN and CN. 
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These values are then taken to initialize the optimiza- 
tion routine. It has been found that the optimized 
values of AN, BN and CN do not depend upon the choice 
of the initializing values. 

The measurement system has been calibrated using 
the above procedure at various positions of the movable 
shert for the-frequencies 5.7, 5. 9%w6;2 and 625° GHZ. 
Figs. 5.6 and 5.7 show the measured reflection coeffi- 
Cient, ue ,and the corresponding adjusted and corrected 


values, [. for frequencies 5.7 and 6.5 GHz, respec- 


AAI! 
tively. It can be noted from the figures that the 
measured reflection coefficient values lie on the unity 
circle after correction. It may also be noted that the 
calibration procedure is frequency dependent as indicated 
by the value of the phase difference between the corrected 
and uncorrected results at the two extreme frequencies 
used, this phase difference being greater at the lowest 
frequency. 

In order to evaluate the active Impatt diode 
admittance, the impedance determined at AA' is trans- 
formed to an imaginary plane at the semiconductor chip. 
This is done by the use of the a, b and c parameters 


which will be defined in the following section. 
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5,4. 2aLnsertion.network characterization 


The insertion network may be fully characterized 


by three complex parameters a, b and c. Three termina- 
tions can be used to fully determine these parameters; 
also systematic errors in the measured results can be 


minimized by using a large number of terminations. 


However, replacing the diode by the other three termina- 


LLOnS,Vi2., an Open, alshortvor.a known Load, may 


change the diode RF conditions. In the case of the Impatt 


diode, impedance below breakdown voltage can be substi- 
tuted for the above terminations as it has a known 
impedance for a given bias voltage. 

When the bias voltage of the diode is smaller 
than V, ,the impedance Z5 of the active layer of an 


b 


Impatt diode is given by 


Z(V) = Ry(V) - j/wC(V) 


Ry (V) + 5X, (V) (5.5) 


where Rp (V) and C(V) are the resistance of the unswept 
region and the bias voltage dependent capacitance, 
respectively. Both Rp (VY) and C(V) decrease with the 
increase of the reverse bias voltage. 

Generally, any two-port insertion network is 


characterized by the fundamental parameters A, B, C and 


D (Pig. 5.3) [44]. The transmission line characteristic 
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impedance is denoted as Zoe The necessary circuit para- 
meters are determined by measuring the normalized 
impedance Zn , looking from the AA' plane towards the 


diode. Zn is given by 


AZ, + B 


oak ZA CEE GD ; 8! 
re) D 
Four such equations are required to evaluate the A, B, 
C and D parameters. Since only reciprocal elements form 


the insertion network, Eq. (5.6) can be simplified to [37] 


aZy+b 
where 
a= A/ (ZOD) 
b= B/(ZD) 
CoaG) 
c = C/D 


| = ADT-"CB = Z D* (a- beyrs 


a, b and c can be determined by evaluating Zn for 
three bias voltages (from the C-V characteristics below 
breakdown) and microwave impedance measurement of Zo at 
the same three dc bias points. This procedure can be 
summarized as follows: 
in At various bias voltages below breakdown the 


Capacitance of the diode wafer may be measured accurately 


with a Boonton /1. Ay UMizyacapacitance bridge. /fheiec-v 
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Characteristics are thus determined for a range of bias 
voltages (V) between zero and just below breakdown. 
Fortunately, for Impatt diodes based on a one-sided 
abrupt p-n junction it is possible to evaluate the 
incremental change in the diode resistance [36] AR, (J) 
(over a certain bias range below breakdown voltage) as 
a function of bias voltage with the help of the C-V 


characteristics. The incremental change (AR)) of the 


diode resistance which corresponds to the change of voltage 


AV (in the interval designated j, oe ea Ue. ee Vee has been 


expressed by [36] 


1 1 
Givpwaes @ Glo 


NEO apie 0. C10) 
D 5 V541 


where 


o(j) is the diode wafer resistivity in the inter- 
Vals jt, 
€ is the dielectric constant (ce =.1.036 
DF/cm tor S18 [36)"), 


SN) and. ety ) are the diode wafer capacitances 


af edk 
at bias voltage of V5 and pee respectively. 


Te Tesistivaty enti tor n-=type, Sl i(in units OQ. cm) 


is approximated using the expression [45] 
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and 1 
ae x(C, + Con (pF) 
VC = c, - Cray (pF) 
AV = ea ~ ve (volts) 


oO 
I 


junction diameter (cm) 


30 


ke=—sa9 constant =—9;.91*«410 Lore Se 


N, (3) = impurity Concentration in, (741,53) imterval. 


Theyapproximationssare valaderor a any PNG 2x mos 


This range is appropriate for diodes intended for normal 
Impatt mode operation at frequencies in the 3 to 8 GHz 
range. At each bias interval j, there will be a specific 
N. (3) and o(4), for which AR, (3) is determined. For the 
penn’ abrupt junction Impatt, the diode resistance Rp 
can be assumed approximately zero ohms at a bias suffi- 
ciently close to break down. Taking this resistance as 
a reference and extrapolating in terms of ARys« LE ies 


possible to arrive at the value R, corresponding to any 


D 


bias’ voltage. Therefore, the diode impedance Z 


D 
(Z, aeons: jXp) can be evaluated at different bias 
voltages. Fig; 5.8 shows xX fand R*sasifunctions of the 


D D 


bias voltage for one of the diodes used in the measure- 
ments. 
eae Reflection measurements are then made at the same 


bias voltages at the desired microwave frequency. From 
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these measurements, Zo as a function of the dc bias may 
be determined. 

Once the a, b and c parameters (of the insertion 
network) are known, the Impatt electronic admittance 
and the RF voltage across the diode can be obtained as 


outlined in the following section. 


5.5 Determination of the diode admittance and RF 


voltage [37] 


When the bias voltage becomes larger than V the 


Db’ 


diode is characterized by an admittance Y. with a nega- 


D 
Live: read part,,andytheicireustsofaFigeay57u5° canabe 
operated as a reflection type amplifier. The nor- 
malized impedance Z 4 of the diode (normalized to the 
characteristic impedance ZG of the transmission line) 


can be measured as a function of input RF signal levels. 


The electronic admittance Y. is then given by 


D 
eae! a)/(b- Z2.)=-G,+ 5B, . (52) 
The RF voltage Ver applied across Ypr the diode 
admittance, can be calculated as follows. Using the 


Doulraary CONGLcion at fA plane, Fig. O.o, ala tie 
definition of the A, B, C and D parameters, the follow- 
ing equations are derived [37] 
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where a is the voltage amplitude of the incident wave 
oEtlegnn we Diane. Subst tutingeraqt (5.8) tintorkg. (5.513) 5 


Vv is obtained 
RF 


at 
Lt 4 e 
Vee a aC Be) eater Oes ay ss slain (5.14) 


Lica Cm Csr Si 2) aCe Hes. (Sle) Yp as well as Me are 


calculated. 


BO Evaluation of «the spass ive CLEGUItE admittance 


To perform dynamic studies on Impatt diode oscil- 
lators, both the active device admittance and the passive 
circuit admittance seen by the diode are to be evaluated. 
In this section, the measurement of the admittance of 
the embedding circuitry, to be treated in Chapter VII, 
will be discussed. The methods developed to obtain the 
parameters of the matrices characterizing the two networks 
are repeated here. But in evaluating the latter network, 
the oscillator under consideration is to be placed at 
plane AA' (Fig. 5.1). Once these parameters are derived, 
then, by transformation of the 50 ohm characteristic 
impedance of the line, Zoe to the plane AA', the passive 
circuit admittance can be determined at the frequencies 


of interest. Thus the’*circuit admittance is given by 
ca, ua 
x 


a a (pus) 
Cc Z2,+b 
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5.7 Experimental procedure 


These measurements involve the determination of 
large - signal Impatt diode admittances together with 
the admittance parameters of its associated passive 
circuitrymin a practical oscillator (to be treated in 
Chapter VIII). The active measurements are performed 
over a range of power levels between -20 dBm and 23 dBm 
and at a number of frequencies lying between 5.7 GHz to 
6a) GHZ.2%e diode RF voltage range of 0.5 to 25 volts 
corresponds to the power levels mentioned above. The 
two diodes, Hp 5082-0437 and Hp 5082-0431, used in the 
measurements have a junction area of eres em* as obtained 
with microscopic measurements of the diameters of several 
damaged diodes taken from the same manufactured batch. 

The capacitance measurements at 1 MHz (C-V charac- 
teristics) are made using a Boonton meter, a diode mount 
with an APC-7 precision connector, a digital meter and 
a high resolution dc-power supply. The measured capaci- 
tance is that of the wafer plus the diode package. The 
package capacitance is accounted for by initially adjust- 
ing the Boonton meter to zero with an empty package 
placed in a specially désigned diode mount. Measurements 
are typically made at thirteen voltages in the range from 
zero to 100 V (100 V is slightly less than geebe Bag. 
5.9 shows the behaviour of the junction capacitance with 


bias voltage for the two diodes. The junction 
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Fig. 3.9 C-V characteristics of Impatt diodes #1 & #2. 
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Capacitances at 100 volts are .31 pF for diode #1 and 
2289 pF™for°diode #24 

Reflection measurements with the network analyzer 
are carried out at a specific frequency and at several 
power levels between -20 dBm and 23 dBm, as follows. 
sa Characterization of the error network, at the 
desired frequency is done at 0 dBm only, because it 
is experimentally found to be independent of the power 
level used. With the precision movable short placed at 
a plane AA' (Fig. 5.1), reflection coefficients are 
recorded for thirteen positions.of the short. This data 
is used to determine the AN, BN atta CN parameters of the 
network. 
23 For microwave measurements below and above break- 
down, the Impatt amplifier is tuned to the proper fre- 
quency. This is done by replacing the movable short 
at plane AA' with the Impatt amplifier. The klystron 
and the TWT amplifier are replaced by a sweep oscillator 
(4-8 GHz) for swept frequency measurements. The diode 
parasitics are tuned to the desired frequency by chang- 
ing’ the’ position of’ the’ slug’until a gain’ of about 10° dB 
is observed on the phase gain display of the network 
analyzer. Then the measurement system is returned to 
its original arrangement with the klystron source and 
TWI amplifier replacing the sweep oscillator. The system 
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the desired power level and frequency. 

Ss Reflection measurements are made with the values. 
of voltages See being the same as those utilized 

in the determination of the diode C-V characteristics. 
The data obtained are combined with the diode impedances 
as calculated from the C-V measurements in order to 
determine the a, b and c parameters of the insertion 
network. It should be mentioned that for each measure- 
ment, the return gain and phase are recorded, once for 
the RF signal directed into the diode and then again 

for the signal reflected back from the RF short (by 
using the RF switch). This procedure enables the correc- 
tion of each measurement to the reference RF short. In 
this case, the drift in the measurement due to the low 
frequency circuitry within the network analyzer is 
minimized. 

4. The diode bias voltage is then increased to a 
value above the breakdown for which a suitable dc 
current is established. The RF power into the diode is 
varied between 23 dBm and -20 dBm (which is within the 
Capability of the Hp diodes used), and the return gain 
and phase are recorded. 

5. Steps 1 to 4 are repeated at the frequencies of 
interest. 

6% The results are then processed by a computer 
program (see Appendix 1) in order to evaluate the active 


admittance of the Impatt diode (Eq. (5.12)) and the 
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corresponding RF voltage amplitude across this diode 
ae loe ts): 

However, in order to evaluate the passive circuit 
admittance seen by the diode embedded in a practical 
oscillator, the experimental procedure (developed above 
to evaluate the parameters of the two networks) is 
repeated. But, in evaluating the insertion network 
parameters, the Jmpatt amplifier is to be replaced by 
the oscillator under consideration and the microwave 
measurements taken below breakdown voltages are to be 
performed at a power laevet -20 dBm. This power level 
is adequate, Since passive measurements are independent 


of the power level used. 


5.8 Discussion of results; examination of sources of 


Sgigene 


The junction conductances and susceptances for 
the two diodes as a function of the RF voltage amplitude 
and frequency are shown in Figs. 5.10(a) and 5.ll(a). 
This representation is useful in the characterization 
of Impatt diode joscillator) and amplifier circuits. 

These curves have shapes similar to those predicted by 
Been the theoretical studies [28,38] as well as the 
experimental results published for Si diodes [36,46]. 
These also resemble the curves obtained by Bracket [32] 


for the small-signal admittance of Germanium diodes 
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SUSCEPTANCE (mmhos) 


=e Sty ae | 
CONDUCTANCE (mmhos) 


Large signal admittance, Yp , of Impatt diode 
#1 as a function of RF voltage and 

frequency at constant de current of 25mA. 
Circus admittance, You as a function 

of frequency, as seen by the diode 

embedded in a 7 mm coaxial resonator. 
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harge siqdnatedons.tance, Yj, , Of impact 
diode #2 as a function of RF voltage and 
frequency at constant dc current of 25mA. 
CLrcniLte sachueeance, \Gpmrsmts ever PuUnCCLON On 
frequency, as seen by the diode embedded 


in a 7 mm coaxial resonator. 
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(Seerbag.  2.15)).) In Pigs. 5.10¢a) and 5.1(a), the wafer 
Susceptance is seen to increase with increasing RF 
voltage , whereas the negative conductance of the diode 
decreases. The latter effect can be observed on the 
characteristics showing the dependence of the diode 
conductance ony the RE “voltage, (Fig... 5.12) for drode #2). 
This criterion results in an observed decrease in the 
efficiency of RF power generation. In oscillator design, 
there is a trade off between the RF voltage level 
obtainable and the diode efficiency. It can also be 
noted that the decrease in negative conductance with 
the increase in RF voltage constitutes the basis of 
stapiLiity im-oscillators (refer “bo, Chapter-Ii3.. 
Figs.9.10(b) and 5.11(b)show the circuit admittance 
seen by the diode vs frequency when the diode #1 and #2 
respectively are embedded in the oscillator to be 
developed in Chapter VIII. 

It was mentioned in Section 5.3 that measurement 
errors arise in the course of microwave experiments. 
The errors associated with the measurement technique can 
be divided into two categories: 

1. The measured data 

2. The network solution, 

In the first category, the main errors are due to 
the inaccuracy in the measurement of the diode C-V 


characteristics. It has been estimated that the reflection 
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“Fig. 5.12 Diode negative conductance as a function of RF 
voltage amplitude for various fixed frequencies 
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measurement data are accurate to within 0.04 dB in 
return gain values and 0.25° in phase values. For 
typical diodes biased at oe =8-255mA@ (Table S2ce 
errors of this magnitude result in conductance and 
susceptance errors of about 2%. For each measurement 
of the return gain and phase, the error due to the 
instrument drift has been minimized by first making a 
measurement and then calibrating by switching over to 
an RF reference short circuit. Evaluation of the net- 
work parameters (AN,... and a,...) based on either the 


calculated reflection coefficients (T \erorweehe 


AA! 
evaluation of the error network, or impedances, (inser- 
tion network) results in an error in the value of these 
parame ters™ by4not® more-than’ 830t(Table 5412) ..oOTable BY 2 
shows Zn! the measured normalized impedance along with 
the corresponding values calculated by using the a, b 
and c parameters in Eq. (5.7). Thirteen microwave 
measurements together with the corresponding evaluated 
Impatt impedances below breakdown are used to obtain 
the optimized a, b and c parameters of the insertion 


network. When twelve or even eleven, rather than thir- 


teen microwave measurements were used to evaluate these 


parameters and then in turn used to calculate the microwave 
impedances at the drop out points, no effect could be detec- 


ted in the measured conductances and susceptances up to the 


seventh decimal place. 
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The effect of the error due to inaccurately 
measuring the diode C-V characteristics is demonstrated. 
by considering diode #1 as an example. The diode admit- 
tance ata 5.9 GHz is then evaluated, first by using the 
measured C-V characteristics, then for 95% of the 
ordinate values, and finally by considering 90% of the 
Ordinate, values “figs. 5.13, 5.14. and,.5.15. show the 
effect jor the error in the GeVscharacteristic on the 
measured diode conductance, susceptance and admittance, 
respectively. It can be noted from these figures that, 
generally, the values of susceptance and conductance 
decrease on reducing the capacitance values in the diode 
C-V characteristics. A computed error in the susceptance 
values of about 5% in the 95% case and 10% for the 903% 
case is obtained from the graphs of Fig. 5.14. The 
uniformity in the amount of error obtained for the span 
of the RF voltages considered is due to the small change 
in diode susceptances in this range (1.2 times). 

At this point it must be emphasized that the above 
example is for demonstration purposes only. In practice, 
a typical expected error in the C-V characteristics 
obtained through experimental measurements is not greater 
than 10):. OSipF #EfSuch tan€erron afises' due tothe *use'of 
an empty package to compensate for the value of re 
This error does not result in a corresponding error of 


more than 2% in the measured diode conductance and 
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Fig. 5.13 The effect of errors in the C-V measurements on the 
measured diode conductance vs. RF voltage curves at 
5.9GHz (diode #1). Deduced conductance based on: 
(A) measured diode capacitances, 

(B) on 95% of the measured diode capacitances, 
(C) on 90% of the measured diode capacitances. 


. 
——a eo ee 
} 
| 
i 
4 
4 % 
a 
a 
' 
‘ 
| _ ; 
— ee ee - a iene 
gs - Ot 
o is “ » ta 
. 


(aro) BIT TAIMA SQaTaov 


mnSivepen Vs 2 outs. at arora 8 


is gosdihon Ge .a9 sount outs Bok nsuae fr 
Divikion reise: hin Senwhed. .(L4% abotb] a 2 : 
| ah ~wedtiationge ebolh bs resem (A 
‘ 2 4nted Lopaen ‘BbOTS dex veeoet aria ie eee 


i aa 


a3 Rie f bd ib bawveser acy Rex Bates +f 


"i 


124 


SUSCEPTANCE (mmhos) 


2, 10 LS 20 25 


RF VOLTAGE AMPLITUDE (VOLTS) 


Fig. 5.14 The effect of errors in the C-V measurements on 
the diode susceptance vs. RF voltage curves at 
5.9GHz (diode #1). Deduced susceptance based on: 
(A). measured diode capacitances, 

(B) on 95% of the measured diode capacitances, 
(C) on 90% of the measured diode capacitances. 
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Fig. 5.15 The effect of errors in the C-V measurements on 
the measured diode admittance vs. RF voltage 
curves at 5.9GHz (diode #1). Deduced admittance 
based on: 
(A) measured diode capacitances, 
(B) 95% of the measured diode capacitances, 
(C) 90% of the measured diode capacitances. 
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susceptance. 
Since the errors in the method of measurement 
are independent of each other, the cumulative error is 
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CHAPTER VI 


INSTANTANEOUS POWER AND FREQUENCY MEASUREMENTS 


6.1 Basic principles of instantaneous power and 


frequency systems 


In order to compare the theoretically computed 
oscillator transients with those experimentally pre- 
dicted, an instantaneous frequency and amplitude mea- 
suring system is required. Since measurements are 
Carried out in the microwave frequency range, the 
system required must be capable of measuring signals 
with nanosecond changes of power and frequency. 

The present instantaneous frequency measuring 
systems [47] are more or less complex versions of phase 
discriminators [48-51]. In all microwave phase dis- 
criminatoxs! (48), Fig. 6.1, thel measured Signal is 
divided by a power splitter (either a 3-dB hybrid or 
hybrid tee) into two channels commonly consisting of 
two transmission lines of different electrical length 
BLy and BL, (8 represents the frequency-dependent phase 
constant of the lines). The signals from the outputs 
of the two channels are then combined in a power com- 
biner (usually similar sto ai power splittem). Because 
of the signal interference, the envelope of the output 
signal is proportional to the amplitude and frequency 


of the measured signal. Square-law detection, unity 
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VSWR (assumed for calculation convenience) results in 


a’net output voltage [48] 


2 
Ein ft) cos B (2, - Q 


Se seg 


where Ey {t) is the voltage amplitude of the measured 

Signal and Ya is the diode voltage sensitivity (matched 

diodes are assumed). 

However, these systems have several disadvantages 
which can be summarized as follows: 

(a) The need for cutting and mounting two rigid 
transmission lines with different accurately 
defined lengths. 

(2) A new piece of transmission line is required for 
any slight adjustments, or a change in the center 
frequency. 

35) When using two magic tees as power combiners or 
splitters a complicated geometry has to be used. 
This results in higher VSWR values which degrade 
the performance of the system, 

Instantaneous frequency measuring systems are 
usually complex and expensive. In commercial instan- 
.taneous frequency measuring systems, where cost and 
circuit, complexity. anreynory the Limiting factors,,sthe 
above mentioned problems are less important. If one 


does not need the sophisticated features of the expen- 


sive instantaneous frequency measuring systems, the novel 
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frequency discriminator, which has been developed during 
this study, gives all needed information with minimum 
requirements on microwave components. Principles of 
operation of the new discriminator together with its 
characteristics and calibration procedures are described 


in this chapter. 


6.2) Cingle hybrad tee discriminator: principle of 


operation 


The single hybrid tee discriminator [52] is schem- 
atuically shown in Fig: 682. It consists of "an ordinary 
hybrid tee, a transmission line terminated by a movable 
short circuit, and two crystal detectors which match and 
terminate the E- and H-arms of the tee, The principle 
of operation is also simple. The measured signal, re- 
presented by an incident normalized wave ays enters the 
hybrid tee via the 1 arm and is split, almost equally, 
into the H, E, and 2 arms giving rise to emerging waves 
b 


b b and bai: The wave b enters the auxiliary 


TRC Pan: ed be dy 
transmission line and travels a distance & before being 


reflected by the line short circuit. The incident wave 


returning to the hybrid tee via the 2 arm, ws given by 


9 12 expi(=25 2) (b.2) 


where by > is the total wave emerging from the 2 arm and 


B is the frequency~dependent phase of the transmission 


130 


) | ie | 

mest Leek ibaa 30 
ro ealeiloust pees 7 By shad 5 a 
ad} datw tedrapos ise. 

badicaal wre | sosniteatad ¢ 


ingitod ek (Se) sciences sat- absdes plgnte: edt 9 Vere 

WRELE MG, AA TO) AICABACD a7 Saf ‘at nunde eteokee = 

aldevomt « vd batantorges ots noteatepmaess 6, 002 bixd yt 
‘Seu docket dotitw EsOND eS He, igdaysn Owe bak, .2ae0aLD grote 
ted. HA) LO ene tea “S one ova Linaey 
ay ,leneat ortutenem wat ieLemnes als 2), nok sere@o 0 
eueae9 To Svev, eaibLewrson I gemhom- ca ‘ct beansaaaa 
‘3 rouse tiiga;ak hes rue i anit, phy 20g bindyd 


elgiontag eay 


aco 


éttaue 


rev palpramy of salt parviy ante &. boty + ates ‘odnt 
aah tbo orig Badin) cept BW sft 7 (BOs: ext eas? ‘ey 
‘eadeg oxdied i eyreJbld ploves: pas ahs nokeatmanaes 
owaw toebsonk, ett 


mia 7 sls +4 


| phipens soda ends ori? yd bedoatter ” > 


ee 


yd sievip wt eee, Cathy BEY vos Sisktyd 9M 09 meer 


: 
ea 


” - ro 
us) 7 


me 


al 


BSE 


“ro7eUTUTADSTp AOusnbezsZ 993 ptaqky etbhuts z°9 POT 


Yq 
Ea 
(WUW-H) ¥ 


LYOHS 
a TaVAOW 


—_—— 


asa 


(7) ANIT NOLISSINSNVUL 


QIYdAH 


IVNOIS 
GaynsvaN 


i, — 
a 
anh is oe he 
ales . 
: oes 
in * 
A 
beat me an 
cet i ag. om 
ae 
e 
4 ie go A 
~ & mie at ae 
a 4 
a&¢ 
t ; : 
ae ie 
Se) iste, 
” si 
: Lrg 
et ett 
Sain ey 3° 
oe . 
7 
+a 
ee ae ak 
p a Sate 
ey ee Cha 
S a om mitt nee] 2 
hi ) y ! een Hine le ee am 
Mal pa eae Dig ne Oe ie i Pe eat rt cic ie ey See =P ee ; 
ricer | ‘ ey ir f: . “. aa j 4 
\ ; ‘ ; © €sae 
“ 
om su) a, oe 
i Has ‘ a 
& i an aoe 
, 
ee ; 
ee ee 4) hon 
sein rbd Ss PRB ee Se BS 
te i 4 * pa 
sick ma d ‘ 
Sag Se Fede ates 4 
¥ 
¢ 
ny 
‘ 
Lom . i s bio 
Cee) Mole se 
i 
: a 
’ % 
rte ie 


g*g) Teepe waptry 


: t ; ae 
4 thy, bee agitated cia ae aN A Ae 
Pe ee Le Pen Piet ae 
i > + . 


i = 


a 


line. The a, wave undergoes similar splitting as the 
wave a), giving rise to emerging waves bo3 Dogs boy and 


b Since the phase difference of the waves b and 


neo 13 


b depends on the factor (282), the waves in the E~ and 


23 
H-arms interfere. The envelope of the signals emerging 
from the E- and H-arms thus depends on the amplitude and 
frequency of the measured signal. 

A scattering matrix representation is used for 
quantitative description of the discriminator proper- 


ties. Because of the isotropy and symmetry conditions, 


the scattering parameters of the hybrid tee satisfy [53] 
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Assuming an incident normalized wave ays match terminated 


E- and H-arms and match termination of the 2 arm for the 


wave aj, waves emerging from the l, 3 and 4 arms are 


Loe 


(6.3) 


g =~ 2 Fa 
Sie a, [Sj DS} exp (-j282) ] (6.4) 
b,=a,5,,[1+DS,, exp (-j282)] (6.5) 
b,=a,S,,[1-DS), exp{—j.262) ] (6.6) 
1 


where D = [1+S..exp(-j262)] «= |Dlexp jg =1. Assuming 
square-law detection, the discriminator net output vol- 


tage is 
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where Ya and Yq. are voltage sensitivities of the diodes 
3 4 


D. ana’ Dy, and 42 is the phase angle of thes 
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The input VSWR of the discriminator, which can be eval- 
Uetecerromend.» (6.4), varies on ther limits 
1+|8,,|+IDl |s, 41° 

dia eZ 


2 

| 
12 < VSWR < eae for Ee STO (6.8) 
|? 1-|S5,1> [Dh dsyod 


if -[s,, +10] Is 
Vows Stee hey 0s bli 


depending on the rare frequency and the length of the 
transmission line. 

The discriminator output voltage depends on the 
Si2 parameter according to_.Eq. (6.7). Since the amplitude 
of the Sio parameter for a lossless waveguide E-H tee 
Poure vacedeto. Ene. ret lection CoOertrtrelents: of the i.. 3, 
and 4 arms by [54] 
si)? 2 S31” : | 


2 
S44 
44 2 
oo biel hammer (6°19) 


only nonmatched hybrids, usually characterized by 


2.0.2 and |S > 0.2, will provide a useful sen- 
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sitivity. Scattering parameters of a typical hybrid tee 


depend very little on frequency in a sub-gigahertz 
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bandwidth (e.g., in the range of 5.7-6.5 GHz, the para- 
metersof a C622A Microlab hybrid tee are as follows: 


(554) -0.08-0.1; | = 0.65-0.60; 0.60-0.54; 


S13! 1Sy gl 
0.46-0.52) and the shape of the discrimina- 


S34 
and |S,5| 
tor response is determined by the cosine term. Thus the 
response curve of the single tee discriminator is very 
Similar tofthat of an ordinary phase discriminator. 

The results of experimental investigations ona 
Single hybrid tee discriminator are graphically shown in 
Figs. 6.3%and 6.4. To eliminate the erfects-om signal 


dependent sensitivities of detector diodes on the 


measured response curve, power meters were employed for 


the measurement of the output signals of the E- and H-arms. 


The polar display of the E- and H-arm output powers is 
shown in Fig. 6.3 for various frequencies. The x-y 
coordinates represent output power normalized to the in- 
put power. Triangles and dots indicate experimental and 
theoretical data, respectively. As can be seen, the 
agreement with respect to absolute values is very good 
iacine the phases do not agree as well. This is because 
of relatively poor accuracy in the measurement of the 
phase angles of the tee-scattering parameters. In the 
polar display of the discriminator outputs, the distance 
and phase angle of a displayed point uniquely determine 
the instantaneous power and frequency of the measured 


signal. The input VSWR of the single tee discriminator 
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Fig. 6.4 Input VSWR of the single hybrid tee 
discriminator. 
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is within the limits 1.35 < VSWR < 2.08 in the frequency 
range Of >.9-6.3 GHZ.— A 6 dB attenuator in the dis- 
Criminator input line will reduce the input VSWR to 
P08 < VSWR <2.20.5 The: discriminator Enput VSWR is 
shown in Fig. 6.4 for two different lengths of trans- 
mission line. The input VSWR lies within the predicted 


range. 


6.3 Instantaneous power and frequency measuring system 


The instantaneous power and frequency measuring 
system is essentially the single hybrid tee frequency 
discriminator combined with a directional coupler, to 
sample the input power together with an isolator to 
isolate the discriminator from the rest of the system. 
Prior to instantaneous power and frequency measurements, 
the measuring system should be calibrated using the pro- 
cedure discussed below. Fig. 6.5 shows a block diagram of 
the system used for calibration purposes. Added to the 
Single hybrid tee discriminator are a J-band (4-8 GHz) 
Signal generator, a frequency counter, two isolators, an 
attenuator (0-20 dB), a TWT amplifier and a power meter. 
The first isolator 1S usedhetovyisolate the signal generator 
beven the rest of the system and in this case the generator 
sees a constant load within the frequency range of inter- 
est. On the other hand the second isolator is used to 


isolate the TWT from the measuring system since reflected 
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Signals could cause frequency and power instabilities 
at the output of the TWT amplifier. The same system 
will be used, in the next chapter, to measure Impatt 
oscillator transients merely by replacing the signal 
generator and the TWT by the oscillator to be charac- 
terized. The crystal detectors used to record the dc 
voltage from the tee E~ and H-arms are HP 423A which 
will be called D, and D 


3 


power of as low as 130.uW for D 


4? respectively. A detected 


3 and 70 uW for D 


corresponds, to de) voltages of 0.4 mV and 0.1 mV, 


4 


respectively. Fig. 6.6 shows the detectors response 
to input signal levels at the frequencies of 5.8, 6.0 
and 6.2 GHZ in\the frequency range of interest. fit 
can be observed from these plots that the voltage sen- 
Sitivities of the two crystal detectors are different. 
For a given frequency, and with an unattenuated 
CW signal, the input power is adjusted to the maximum 
power required. Then, at the various power levels 
(5-90 mW), the D., and Dy responses are recorded 
sequentially. These steps are repeated for the fre- 
Guencies 207%, (0.6, oO y.0. 1d. ond Gwe GHZ. The 


constant power and constant frequency curves displayed 


in Fig. 6.7 will be of help for instantaneous power and 


frequency measurements. It can be seen that the curves 
are not exactly circular due to the different sensiti- 


vities of the crystal detectors. They are similar in 
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Fig. 6.7 Constant power and frequency response of the 
instantaneous measuring system. 
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shape to the curve derived theoretically using the 
scattering matrix concepts and the curve measured with 
the help of two power meters (Fig. 6.3). 

The calibration curves of Fig. 6.7 show the power 
and frequency characteristics of the measuring system 
under CW Signal operation. In order to use these curves 
for instantaneous (power and frequency) measurements of 
oserllator transtents, the icapabidji ty of the system to 
instantaneously measure the power and frequency of pulsed 
modulated signals has to be demonstrated. For this 
purpose, a Similar calibration system was set up for 
pulsed measurements (Fig. 6.8). The differences were 
Minor consisting essentially of the addition of a pulsed 
modulation stage in the RF path. This helped further to 
confirm the accuracy of the system in practical opera- 
L1ON. belie G.. Goa,;5 ie isolator following: the sagnal 
generator is essential in presenting constant load con- 
ditions to the generator that maintains a constant input 
RF signal (6 GHz), especially during the ON conditions 
of the PIN line. Under pulsed operations, the power 
meter is used to measure the average input power. Know- 
ing the duty cycle of the modulated signal, the peak 
‘power could be easily evaluated from these measurements. 

The RF signal was modulated with pulses of ampli- 
tude = 3. volts, duration =] %sec.,and repetition rate 


of 500 kHz. A PIN diode was the modulating element 
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The two RF envelopes recorded by the 

X- (top trace: 20 mV/cm) and Y-detector 
(bottom trace: 10 mV/cm). A modulating 
pulse of 3 V-amplitude, 2 us-period is 
used to modulate an RF source (6 GHz) by 
means of a PIN diode modulator. 
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(restricting the rise time of the resulting RF envelope 
EOmavoUtC® Fl oons The RF envelope recorded by the two 
detectors 1s shown vin Big... 6.9. ‘The ‘top trace corres- 
ponds to the D, detector and the bottom trace corresponds 
to the D, detector. In order to evaluate the instan- 
taneous frequency, at various time instants, starting 
from the leading edge of the two RF envelopes, the 
amplitude is recorded OFF the RF envelope of each 
detector. These values are plotted on the calibration 
curves of the system (X-Y polar display) and by inter- 
polation, the instantaneous frequency and power are 
determined for each time instant. These points are 
numbered (1 jr2,,.4 .gsecte.acorresponding to different time 
intervals along subsequent points on the RF envelope 
graphs. Here, the points of interest, are those on the 
top portion of the RF envelopes which are numbered 20, 22, 
24 and 26 (Fig. 6.10). The measured instantaneous fre- 
quency does not differ by more than 4 MHz from that 
measured on the counter (6 GHz). It can also be noted, 
that the measured instantaneous peak power never differs 
from the figure calculated from the measured average 


power by more than 4%. 
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CHAPTER VII 


THEORETICAL STUDY ON TRANSIENT AND 


STEADY STATES IN OSCILLATORS 


LaverlLneroductron 


The dynamic behaviour of a microwave oscillator, 
which is characterized by the time dependence of its 
Pistateaneous Premiendy and amplitude,is described by 
Eq. (2. 00)". two 1nlttral condrtr1ons are requrred “in 
order to find unique solutions for the amplitude A and 


the frequency w(t) related to the known x by 
w(t)=0 + dd/dt =Wtx. 


The physically accepted initial conditions are definable 
either by the conditions at the start of oscillations 
(i.e. the start of the transients) or by the conditions 
immediately following the transient stage, referred to 
as the "almost steady" state conditions. The solution 
based on initial conditions that characterize the start 
OL Osc#ilation bulld wp, will be called“the "forward 
recursive method, because it progresses along with the 
balla up Ol Isc rel atrous: wroltice roser. tations are 
initiated by thermal noise generated in lossy components 
of “the "circuit, an amplitude “value very close to zero 
may be used as one of the conditions for the forward 


recursive method. In this case, the starting value for 
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w(t) is determined mainly by the characteristics of the 
passive and active components of the oscillator. The 
alternative solution will be called the backward 
recursive method. Here, the initial conditions are 
taken at a point close to the steady state value. The 
forward recursive method can be easily implemented in 
the case of an oscillator with a simple pasSive circuit, 
e.g. a simple parallel GLC circuit. When applied to an 
actual oscillator circuit, the numerical solutions based 
on the forward recursive method very often diverge. 
Thisreouldeberduesetotthesfactathattthe buibdiup .of 
oscillations is exponential, and the numerical solution 
fails to converge to the steady state on account of 
cumulative round-off errors during successive iterations. 
On the other hand, the backward recursive method has proved 
to be a very effective method for the evaluation of the 
transient behaviour of practical oscillators. 
Differential equations are often numerically 
solved by transforming them to difference equations. 
Solving difference equations by backward schemes is 
well known in the theory of recursive computations. 
The first algorithm for such a backward scheme was 
Roe by Miller “in, 1954) to.,compute values, of the 
modified Bessel functions I, and the procedure is 
generally known as the Miller algorithm [55,56]. fThis 


algorithm has the distinct advantage over the forward 
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recursive scheme of being stable since cumulative round- 
off errors in forward computations seem to grow with the 
number of operations, as mentioned earlier. 

This chapter describes the methods for obtaining 
numerical solutions of the equations characterizing the 
steady and transient states of an oscillator. In the 
first section, equations and circuit parameters (evalua- 
tion of passive and active d-coefficients and w) common 
to both the transient and steady states are derived. 

The next section outlines the numerical steps followed 
in order to arrive at the steady state solution. Subse- 
quently, the numerical algorithms developed are applied 


for computing the oscillator transients using measured 


large signal admittance for the Impatt diode #2, and also 


by utilizing the data published by Scharfetter and 


Gummel [28]. 


7.2 Derivation of equations and parameters describing 


Steady apd tyansient. states in oscillators 


The steady state amplitude and frequency of the 
oscillator are needed to provide the proper starting 
conditions for the backward recursive method. Since the 
amplitude and frequency of oscillations are constant 
under steady state conditions, all jer ivaeices in the 
describing equation (2.66) are identically zero. Thus, 


the steady state operation is described by the following 
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set of two nonlinear algebraic equations 
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Here, the R's are nonlinear Lunctwons in A, x and the 
quasistationary frequency w. 

Generally; it is extremely difficult to solve 
the two nonlinear equations (Egs. (7.1)) in the two 
unknowns. On the other hand, one nonlinear equation 
in one unknown can easily be solved using available 
computer routines. It has been found that the proper 
solution can be obtained by reducing the problem to the 
solution of two nonlinear equations in a single unknown 
by treating the A variable as a known parameter. This 


requires arranging Res as follows: 
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where Ran are nonlinear functions in A and w. Therefore 


Eqs. (7.1) *avev’transformed to 


La . 
eb ee S20 (7. 3a) 
n=0 
aru = 
nex t=, 0 (7235) 
n 
n=0 


and the coefficients BH and Co are nonlinear functions in 


A and w as shown in Appendix It. 

Selecting values for A and using the predeter- 
mined value of W, to be evaluated in section 7.2.3, 

Eqs. (7.3a) and (7.3b) represent two equations in one un- 
known, x. Each equation can be then solved using known 
numerical routines. Generally, the sets of the eleven 
complex roots are different. The values of A have to be 
selected in such a way that the real roots of the two 
equations coincide. 

Most of the numerical routines used for the solu- 
tion of nonlinear algebraic equations, such as those of 
Eqs. (7.3), require analytically defined equations and 
thus the R coefficients have to be determined by means 
of the d-(spline) coefficients of the passive and active 
‘eircuittsy “hefquasistatiionary frequency) \v)Has %al'so rto 
be determined because the Be and Ce coefficients in 
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7.2.1 Evaluationsof the active device d-coefficients 


To evaluate the d-coefficients, the measured or 
theoretically evaluated device large-signal admittance 
data Yun havetto benexpressed,as/in®Chapter II, in the 


form 
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Next, the natural bicubic spline is fitted through the 


G eh e 5 


11°°°° Sun Pay comer 


corresponding to two sets of grid points, namely 


MeN functzonal values; e.g. G MN’ 


Ayreee sh The computer routine returns the 
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coefficients of the bicubic spline interpolate on the MN 


subrectangle [Ay jeitewas 4x4 matrixes Whe 


Avy) * [eye Oy 


spline interpolate is then given by 
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During the buildup of oscillations, the device 
conductance and susceptance for a given amplitude A and 
“frequency w, which lie between known amplitude and fre- 
quency values given by the device eRe characteristics, 
may be evaluated by the use of Eqs. (7.5) and (7.6). 


These can be transformed into the convenient forms 
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7.2.2 Passive circuit d-coefficients evaluation 


greeerFy (fox 


G(w) and B(w) /w) and the corresponding coordinates 


The set of N functional values FL,F 


Wi rWoreses Wy are fed into a natural-cubic spline sub- 


routine and an output array c of dimension (3, (N-1)) is 
obtained. In the i interval (Oo. , @ ) the spline 
N’ N+1 
interpolate, midaeon is given by 
4 Teal: 
ae 8 Dh, 9D 
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In order to evaluate d-coefficients of the passive 
Guecmmee conductance and susceptance Eq. (7.17) is: trans- 


formed to 
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where Gy and By are the passive circuit conductance and 
susceptance, respectively, and the d's are functions of 
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circuitry at the operating frequency. In practice, 
oscillations start with an amplitude very close to 
zero and the steady state frequency deviates only 
Slightly from | so that it Ws adequate to evaluate & 
at A= 0. 


Since the active device susceptance B(A“,w~) is 


aeeunctcLon of of and Ae and the pasSive circuit suscep- 


tance yan is asfunction ariel w may be evaluated as 
C 
e —2 —2 
BlA=0, Os), + Bl@e ea. (7. 29) 
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for Ghee =ntervals , N;NF1lPproh,ctorpy aswel? as “those ‘of 
the active device for the various subrectangles 

(MX N)'s, are evaluated in the previous section. Subse- 
quently Eq. (7.21) may be solved numerically in w using 
‘computer routines based on a Newton-Raphson iterative 
technique. However, under small-signal conditions (A=0) 
M=1 and there may be various values of N corresponding 


to the various frequency intervals. The routine can then 
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give different values of w for the various values of N. 
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calculated by the use of the corresponding N coefficients. 


the value w should be taken as that 


7.3 Steady state solutions 


Having evaluated the d-coefficients for active 
and. pasSive circuits, as well as determining 0, A is 
assigned a range of values lying within and conforming 
to the device characteristic (G and B for various A's) 
for thisevalnecofn ws Proceeding further, dor each of 
these chosen values of A, eleven complex roots of x are 
computed for each equation of (7.3) using the same 
YOULInNe wtliszedatonethesevalua&ionsof.@, eSance tthe 
steady state frequency is a real number, all complex 
roots are disregarded, and only physically acceptable 
values for x are considered in steady state calculations. 
Pipiacsbecnsobserved (Chat sLer a giver oscillators circuit, 
x has one negative value that leads to a converging 
solution. Also, it has been found that positive values 
of x cannot lead to a converging solution. The steady 
state solutions are the values of A and w+x which satisfy 
both the equations of (7.3) simultaneously. 

A flow chart demonstrating the different steps 
to be followed in order to evaluate the steady state 
values A and w+x is shown in Fig. 7.1. The calculated 


amplitude, A, and frequency W+ x may be used as the 
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Conductance & susceptance data as a 
function of A and w for the active 
device are specified. Conductance & 
susceptance data as a function of w for 
the passive circuit are specified. 


The active device b-coefficients are 
calculated using a bicubic spline 


computer routine. The passive b- 
coefficients are calculated using a 
cubic spline computer routine. 
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initial boundary conditions’for the start of the 


transient solution in the backward recursive approach. 


7.4 Transients in Impatt diode oscillators 


The transient response of an Impatt oscillator 
can be analyzed by solving Eqs. (2.66) which are non- 
linear differential equations in A and x. The time 
dependence of A and x (Ww taken to be constant) may be 
determined only by simultaneously solving Eqs. (2.66) 
using a digital computer. This has been accomplished 
through the use of a subroutine which utilizes a Runge- 
KUteorrOUGeEN OLder algoritnm, [$5/,50).. fe 1s often 
difficult and cumbersome to determine the proper size 
of the integration step in Runge-Kutta type methods for 
numerical integration of first order nonlinear differen- 
tial equations. The sub-routine used adopts a modified 
version of the Runge-Kutta method by Merson. This 
modification gives an automatic and rapid method for 
determining the interval (step length) to be used in 
order to obtain a predetermined accuracy. An accuracy 
close to that of the machine is used during the initial 
steps of integration in order to reduce accumulation of 
round—off errors. 

As pointed out before, there are two approaches 
for the evaluation of the dynamic solution, namely the 


forward recursive method and the backward recursive 
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d-coefficients of the active device 
and pasSive circuit are specified. 
Evaluate w, A and x as in the steady 
State case of bag. #.5. % (time). and 
H (initial integration step) are 
specified. | , 


Determine..M,...bland.accorndingly..the.cde 
coefficients. 


Evaluate Y_and Y 
Evaulate A” and x using Eqs (2.66) 


vt 


A & X are evaluated. 


es é [x (T+AT) ~x (T)l<jG? 
= A 0 — 
oo) : BIATHAN- ACTII<I1G-2 
7.2, Overall itterative-scheme for the calculation 


of the dynamic solution by using forward 
iterative scheme. (Appendix III). 
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d-coefficients of the active device 
and pasSive circuit are specified. 
Evaluate w, A and x as in the steady 
State Case, OF Fig oy. te come), and 
H (initial integration step) are 
specified. 


Time (T) and integration step (H) are 


Specified. 


Determine M, N for which d-coefficients 
are chosen. Evaluate c.9-coefficients. 


Evacuate Yana Y. 
Byvaluave fe tand x susing Ha. (2.66), 


A and x are evaluated. 


|x (T#AT) -x (T l<7g8 
is &IAT+AT)-A(T)I<IO 3 


Fig. 7.3 Overall iterative scheme for the calculation 


of the dynamic solution using backward 
iterative method. (Appendix III). 
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method. The steps followed in the two approaches to 
obtain the time dependence of A and x, are shown in 


Faqs fs 2eand’ 25. “(Appendix —£ Il)"; 


To demonstrate the applicability of the forward 
and backward recursive solutions, the Impatt oscillator 
transients are evaluated for a hypothetical oscillator. 
In these calculations, the Impatt diode, or more pre- 
cisely the semiconductor chip theoretically investigated 
by Scharfetter and Gummel, is assumed mounted in the 
Sy microwave package [59], ‘(contact area = Tose cm) 
and operated in a coaxial line-type resonator with a 
Single ‘tuning ’slug. The external circuit, consisting 
of the resonator and the lumped package parasitics, is 
tuned to present at the chip contacts a conductance of 
ine mho in parallel with a susceptance that tunes out 
the chip small-signal susceptance (corresponding to a 
bias current of 200 AVERT ) at / <4eGHZ . 

Fig. 7.4 shows the dependence of the active device 
negative conductance and susceptance upon the operating 
frequency and RF voltage as published by Scharfetter and 
Gummel [28]. It also shows the dependence of two passive 
Circuit admittances on the frequency as seen by the 
device. These admittances are calculated for two sets 


of parameters for a 7 mm coaxial type resonator shown 


in Fig. 7.5avand described an detail in Chapter VIII. 


ut) es 


a2 BS speorggs Ow? edt tg = 


ay 
ot nworts’ 61s yt ire, A 


wie? @d4 Io Y2tt beso Pbaai s va 
- tack 

tI 7 Qo 7: P| Ons ’ snccilicra | 

Jno . rorisouyel . ot bedi 


mie YO ab a6 sagt ant 


7 Bl J " ori3 
Lt Danwoa 
i : : rs 2 
| to 8 ‘opeabig | ove io 
f 4 melt tiehele a at beds: 
ye 7 
ite it. sQeil a eal ti 
| 
j cue Bag! oe 7s be bite idiancaes 
. - me, - fe = 
) Gino is 2348-232 


cys tats Ss 7eIGSsLDeUaA..& ae to GE at 


, 190) dan ephay swbane 5 : 


4 — oe 


yd botleiiduq —s 9Pps7to¥ +s bas” 
bnsash out WwOorls cade af “eer ifr 


: rar, S Pai 
Hoon ae VYoouscosrts iy ago eeomags bk fos “a i: 


ot Boteluoiss sé esonnautmiba ae mit .5o 


5 


; ar 7 Yt a + ; 
Ks ‘sIIQee 2 Says lelxsco eae tsa 


oe : 


in 


i = 
eee | 


i a etd ‘at ihaset : ah r scanned iin bis 


E*. oe a 


ss 
i 


161 


Also, indicated an Fig. 7.4 are.theM*xN values of 
the various subrectangles corresponding to each set 
of the d coefficients computed. It is clearly seen 
that for the five RF voltages and six frequencies 
available, there will be, respectively, four M and five 
N values. 

Ew Gee 1 DF Yo is the transformed circuit ad- 
mittance as seen by the packaged diode, and CS and Lo, 
are the package capacitance and inductance, respec- 
tively. A computer program is developed to calculate 
the circuit admittance as a function of frequency using 
the two sets of parameters presented in Table 7.1 and 
the equivalent circuit, of grig.) 7.5b as“an input. } in 
these calculations, the discontinuities due to the single 
slug are included as appropriate discontinuity capaci- 
tances [60]. The resonator parameters are adjusted to 
tune out the total small-signal oscillator susceptance 
at 7.4 GHz. As mentioned earlier, the forward recursive 
method sfaiis: to converge tou soractacal roscillator circuits. 
The results obtained on the basis of the backward recur- 
Sive method for these oscillators, are shown in Figs. 7.6 
and 7.7. The “amplitude origin ‘is chosen “tobe Meee volts 
which = about the same as the noise threshold level. 
Comparison of these figures shows that the frequency starts 
to increase appreciably only after the RF voltage has 
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Big. 7.4 (aie fImpatt diode’ admittancegas a functhon of 
RF voltage amplitude and frequency at a 
current density of 200A/cm (After 
Scharfetter tal / s[238)9. 

(b) Passive circuit admittance, Yo, as a 
funchion, omeerequiency for circuitvA 
(Tables 731) 

(c) Passive circuit admittance as a function 
of frequeney 1om circudt By (Table *7.1):. 
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Ry) Ror Rye XL5 and XL2 are the resonator parameters (Table 7.1). 


7.5a Cross-sectional view of a 7mm coaxial-type resonator embedding the Impatt 
diode. 
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Fig. “%.55 Equavalent circuit of an Impatt 
diode oscillator. 
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Table 7.1. The two sets of resonator parameters used 
to evaluate the passive circuit admittance 
as a function of frequency for circuits A 
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calculated steady state frequency is about the same 
(within 5%) for the two resonators used, even though the 
time needed to reach the steady state frequency is 
different. It is also obsesuvedgathat the oscillation 
rise time’ is different,for the two oscillators. This 
demonstrates that the interaction between the circuit 
admittance locus and the device characteristic curves 
plays an important role in the transient behaviour of the 
oseillater. Therefore, 1t is possible, by changingtthe 
resonator parameters, to optimize these transients. 

The characteristics of the Impatt diode oscilla- 
ton consisting of diode #2 mounted in a / mm .coaxial 
line resonator whose design is described in the next 
chapter, are used to calculate the transients. For this 
example, the backward recursive solution is obtained by 
using the admittance parameters measured in Chapter V 
(Fag; 95218). Imp Rig. 7.8 curve (a)shows the RF power 
build up for this arrangement. From this figure, it .can 
be deduced that the oscillation rise time is approximately 
5.4 ns. The frequency build up is shown by curve (b)in 
this figures "rigs 7. o formc™ tie™ baste .oOr “tie Cconmparison 
between the transients calculated by using the backward 
recursive solution and those measured experimentally for 
the Impatt diode oscillator developed by the author. 

In order to understand the behaviour of microwave 


oscillators, the evaluated instantaneous amplitude and 
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Fig. 7.8 Build up of power (a) and frequency (b) in the 
Impatt diode escillator,usingsdiode 2 (Ch.« Vids) . 
The measured oscillator admittances are shown in 
Bag ds Dol, 
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frequency are used to obtain the dependence of the device 
admittance on the RF voltage. This has been described 

in Chapter III as the device line [17]. The calculations 
are performed for the two hypothetical and the developed 
Impatt diode oscillators. A computer program is developed 
for this purpose in which the active and er eere d- 
eOctficients, A(t), x(t)» and w are fed as inputs. The 
program outputs both conductance and susceptance values 
asia Lunction of the RF voltage. These calculations are 
also useful, since they give the behaviour of the device 
admittance with time during the transient stage. Fig. 
7.9 shows the results obtained for the device line toge- 
ther with the admittance locus of the two hypothetical 
oscillators. The results obtained on the basis of the 
measured oscillator admittances are shown in Fig. 7.10. 
In both figures, the arrowheads attached to the locii 
indicate the\direction of increasing w in case of the 
admittance locus, and the/direction of increasing RF 
voltage as well as time for the device. The intersection 
between the device line and admittance locus determine 
the steady state frequency. The steady state frequency 
for the hypothetical oscillators is about, 7 GHz, and 

this can be easily seen as the intersection points which 
are, close PonymGHzan(® j=n7.4 GHz). It can’ be noted from 
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Fig. 7.9 Device line and impedance locus superimposed on a 
part of Impatt diode admittance characteristics 
(Fig. 4.9) published by Scharfetter and Gummel 
(28) aew = 7.4GHz). 
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tare the Same paths It can also be-noted, according to 
Pigs. o.oo ang 9/00) that at rtherseartiof, oscillation 
build up, the frequency changes slightly until the RF 
voltage builds up to more than one volt; thereafter a 
noticeable change in frequency occurs. These changes 
are different for various oscillators as shown in Figs. 
ieorand  /.10. ¢ Im case ob the oscillator developed, the 
steady state frequency is 5.804 GHz and w = 5.886 GHz, 
and the change occurring in frequency from that at the 
Start orgoscillationtbueld up is only about, $1.6 MHz. 
This is unlike the frequency change in the hypothetical 
oscillators, which is about 377 MHz. In general, the 
behavioureoh OScitiatiom build up is sSimilarjin alg 
oscillators. 

A stable oscillatory condition is possible pro- 
vided that, in addition to the conditions imposed upon 
the device and circuit admittances, the following 
requirement be fulfilled [61] at each intersection on 


the admittance plane plot: 
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Device line and impedance locus superimposed 
on the measured Impatt diode admittance 
characteristics, diode #2, assof Fig. +5<eil 
(Ww = 5.886GHz). 
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Using the parameters involved in Eq. (7.22), for 
the antersections of Figs. /.,9¢and./ 210, it can easily 
be shown that these operating points are stable. Appro- 
Ximate values for these parameters are given in Table 
eee Le dG) / 0A is positive, then stable oscillations 
can only be obtained if both dB /dw and dG /dw are nega- 
tive or else dB ./dw must be negative. It should be 
noted that although the transients for the two hypothe- 
tical oscillators are different, the angle between the 
device line and admittance locus is about the same, and 
equals tou. 07. This angle for the developed oscillator 
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CHAPTER VIII 


BAPE REMENTAL oO LUDTES*ONse INPATIODLODE OSCE LEATORS 


The emphasis so far has been on developing a 
nonlinear model to describe the behaviour of microwave 
active circuits embedding a certain class of negative 
resistance devices. This theory has been applied to 
Impatt diode oscillators in order to demonstrate how 
the general theory can be adopted to a particular device, 
and also to develop an insight into the characteristics 
of these oscillators under transient and steady state 
conditions. 

For experimental studies, an Impatt oscillator 
has been designed for pulsed operation. The oscillator 
design and experimental arrangement for measuring the 
transients is discussed in the first section of this 
chapter. In Section 8.2, steady state and transient 
studies are presented for oscillators using two Impatt 
diodes whose large signal admittances have been presented 
in Chapter V. The experimental results are compared with 
those predicted theoretically in the previous chapter. 
The operation of the bias modulated oscillators is 


diseussed in the Lact section. 


8247) Design Of Impatt, dicde oscillator and the measurement 
system 


The system to be used in the transient measurements 
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CONnS#sSeSeOGed=teSOnatolr.circult, the bias circuit and 
the instantaneous frequency and amplitude measuring 
arrangement. The latter has been discussed in Chapter 
VI. The measurement system is shown in Fig. 8.1 and 

an expanded cross sectional yviewsot the oscillator is 
given In -Fig.9 0,2. Thispescritatone 1s ‘equipped with 

a unique bias arrangement to make it suitable for pulsed 
operation mode. Two special terminals, one for monitor- 
inggthe diode currents and another for monitoring the 


modulating pulse have been provided. 


S.i41 §Thes resonator cilrecuit 


The Impatt diodes, whose specifications are given 


in Table 8.1, were separately end mounted in a single 


slug 7 mmicoaxtalvcavity (Pig wes 3).  Tolvensure #that only 


the TEM mode exists in the cavity, even at the second 
harmonic of the RF frequency, a 7 mm airline was used 
forithe cavity construction.) These™diode cavities’ are 
of the single transformer, fixed frequency type, so that 
ther could be represented by a simple equivalent circuit. 
A change of operating frequency requires a change of the 
center conductor; however, this was a straightforward 
procedure. The single transformer is designed to match 
the 50 ohm characteristic impedance of the transmission 
line to the impedance at the steady state operating 
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diode #2. The impedance transformer is designed using 
known techniques [61]. 

In addition, the discontinuities of the trans- 
mission line due to the single slug are included as 
appropriate discontinuity capacitances [62]. The 
transformed impedance can be put in parallel with the 
diode impedance, i.e. wafer and package, to obtain a 
complete model of the circuit. An approximate equivalent 
circuit representing this transformation is shown in 
Fig. 7.5b. The diode is biased by using a special arrange- 
ment to be discussed in the next section. The diode 
current is monitored by measuring the voltage across a 
15.0 seriesyresiscomun sche, biasing circuit. 

Tt snowldsber notedsthat’ ifsplanes A and|ZB .(Fig. 
8.3) coincided with each other, then one side of the 
diode would be effectively grounded. This configuration 
is suitable for CW operation or for pulsed operation 
with the pulse generator connected to the ground nega- 
Cive  temminalisyon thersdc power Supply. This last case 
is not acceptable because the pulse generator is then 
loaded by the capacitance of the power supply negative 
terminal, to’ ground: the value of this capacitance for 
the HP 6207B dc power supply used is about 1500 pF. 
Therefore, for pulsed operation, planes A and B are 
separated by a 2 mil thick teflon sheet and the pulse 


generator is connected between the heat sink and ground. 


7 


esate : 


iy) 


. = 
: Tis J = 
x . 

‘ 

7 

+ 


es ‘ . . pS aoris ss al .;) bom: if 2 
2 | ; a ome 
ta ms es ten 3 


a 


o 2 . tes j 
. 4 aes FLUO 
LEs840 yd uk 1otsiect -eohtee ” 
« Berrod ot BSt sorta sa i 
: alas a 7 a 
c pha, an i < SS Fg hg ve bebtontan) £ 


patwoxg Yovitsgthe ee bide 7 ‘sales 


} <3 
ba: hoitexssqo WO sot. eicied me 
. j ; J , % z 


be 4 aScrQO > fieasy salve sii 
e tawby _ ‘Re atestiniss avi: 
‘ue SAd Geusoed. | oidasgeond Age 

; a0 » hs gt ne ats . 
¢ 22. Fo con sa Lng edz | Sebso 
ren fr A? Jo ait av onc here ee 


sav seweq . ob avogi 


J a a = 
g bas A ganged F055 qo bosing 40% ‘ aie 4 


3 ee ae or 
4? bere feade noLies sotdy fie S * d bet, 
ine Seqt-ed? Asewzed | sesononn ae I 


1 


180 


ee ares NI] 
2 EROS CN TINE AR Boh aie ile aa) ET ae 


*apotp 4yedwuy pejuNnow pus 
uatm ARZTAeD TeTxeoo ww, ou FO MATA TeUOTIOSeS-ssoiIp ¢€°g °HTd 
qd ANWId 


d ANWVId 


adord 


is Com 
Si aa er ani 

fans wee 
Ps ae 
Rae 

Ke, sanite da sie aig 
Oi eisai Shas Seis ae 


—— 2 a 
| YHWHOASNVAL 


SoS SSeS OSOSSSS SSS BSS OSA TWNIWU SL 
eS es >Sese eS. <><] SVId-oOd 
= LATIOO PSOSeSOSOSOS OSES SS OS OS OSE 


GONVCUdWI 
Se 


aa 


ah 
arerts BEE ee nee oe te enw 


oVid od 


NOTdaL am 


LAdLAO — 


La VGTED Ya SPY mC wey 
9 i] . 7 


VAP 


a 


Rin? § es, 


ay wi 


472" 


iat 


£4 


a4 


e° 


e 


; JEVNELOMIES 


* THEROVECE 


rh 


oO? 2 
a ~ od 


Area 


In addition, there is a capacitance between planes C and 
D so that the total capacitance between the heat sink 
and ground is 60 pF as measured on a #71A Boonton meter. 
At 5.8 GHz this capacitance is equivalent to an impedance 
of about 1.44/7 ohms. The penalty for decreasing the 

RF impedance further is a prohibitively high time 
constant seen by the pulse generator. Noting that the 
output impedance of the pulse generator is 50 ohms and 
without including any other capacitances, the time 
eonstant tor the cCireurt 1873 ns P@here astardistinct 
trade off between the value of the RF short and the time 
constant associated with the circuit. The use of a 50 
ohm termination in parallel with the output of the pulse 
generator reduces the time constant of the circuit to 
ims ns. The value of 60 pF chosen for the RF blocking 
Capacitance (heat sink to ground) is considered near 


\ 


optimum for the circuit used. 


Se 120-2 the=bias Gireuit 


The impattreaL.oges 10 OSC Llatorwerrcurts are 
usually biased through a biasing element called a bias 
tee. This arrangement is adequate for CW operation. 
However, under pulsed operation the capacitance in the 
RF arm of the bias tee (= 12 pF) and the inductance in 
its ac arm (2 20% nH at’ 50. MHz) * @degrade 'the™“rise 


time’ of the diode bias Current which in turn affects 
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the rise time of the measured modulated signals (RF 
envelopes). To alleviate this problem, a different 
arrangement has been used as’ shown in Fig. OrnaeLn 
which the de voltage is supplied directly to the diode 
through a resistor. Two feed through capacitors 
(="12000pE.); acting as RF shorts, isolate the RF from 
the de. power supply...This, together with a 15 ohm 
resistor inserted in the de circuit, prevents insta- 
bilities which might result in bias oscillations [64]: 
The bias circuit as seen from the pulse generator 
terminals can be modelled as shown in Fig. 8.4. 

The major coneipeee ng factormaceetue, Length of 
the rise and decay times of the diode current is the 
time constant of the parallel RC circuit loading .the 
pulse generator. This is shown in Fig. 8.4 as a 
parallel combination of the source resistance Ro the 


terminating resistance R, and the heat sink to ground 


dk 


capacitance C The Ya se ‘time. OF stg Sane Cr rout can 


H° 
be decreased in one of the two ways: 
1. by decreasing the value of Ch 


2. by decreasing the value of R.- 


Cy acts as an RF bypass for the cavity. Decreasing this 
Capacitance below 60 pF results in an excessive amount 
of RF radiation losses through the back of the cavity. 


Decreasing the value of Re by further Toading the pulse 


generator results in a very low voltage pulse across the 
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PULSE 
GENERATOR 


where Ee diode resistance 
Ripete bias eurrentsmeasuring, resistor, +150 
Cyu= heat sink Capacitance to ground=60pt* 


C = feed through capacitor =12000pf 

Lp = DC supply input inductance =160nH at 
50MHz** 

Ro = pulse generator source resistance = 
502 

ies load resistance =50 

R,, = load termination = 502 


* measured with Boonton 71A inductance- 
capacitance meter 

** measured with Hewlett Packard HP8405A vector 
voltmeter 


Fig. 8.4 Bias circuit as seen from pulse 
generator terminals. 
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diode. There is. a trade off for high pulse repetition 
generators between short transients and high amplitudes. 
With these considerations, the compromise values of 

Re and Cr aresthose’ givenvinyFig. 8.4.7 LDtemayobe 
stressed) thati thesetresultsi|are\ restricted to the 
circuit used in the experiments reported here. 

It has been assumed in the above discussion that 
the RF oscillations start as soon as the pulsed current 
startsto flow. This assumption is seldom true as will 
be seen later in this chapter. 

The following are the various possibilities for 
biasing the diode under pulsed operation: 

iwebrasing at "ORV ido (1. e.,.nurGgc, bias voltage) 

2. biasing at the breakdown voltage 

3. biasing below the oscillation threshold but 

above the breakdown voltage 

4), biasing vat the oscillation threshold. 

A trade off exists in high-prf pulse generators between 
fast rise times and high pulse amplitudes. For the 
required rise ence Otelesso es than «2. anc eon yO] Vv. 
pulse generators are commercially available. Biasing 
below the breakdown voltage and pulsing the diode into 
avalanche breakdown results in fluctuations at the 
leading edge of the diode current pulse. This is due 
to the statistical nature of the avalanche process. 


The above considerations preclude biasing the diode 
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below the breakdown voltage. Biasing below the oscilla- 
tion threshold introduces some finite time delay before 
the current reaches the threshold level of oscillations. 
Therefore, a suitable biasing point for bias-modulation 
applications is at the oscillation threshold. This is 
in spite of the fact that there is always a time delay 
between the RF envelope and the diode current even 

at. tnizs Condition; biasing fat tne -oscrliation threshold 
at least minimizes this delay. This biasing point has 
the added advantage that the diode is already dissipat- 
ing one or two watts; hence the temperature fluctuations 


under pulsed conditions are minimized. 


8.2 Transient and CW measurements on the oscillators 


The two diodes whose large signal admittances 
have been measured and reported in Chapter V are 
embedded in the coaxial cavity discussed in Section 8.1. 
The oscillators are designed for frequencies around 
5.8 GHz. The diodes, under both CW and transient con- 
aeons operate well within the frequency range for 
which their large-Signal characteristics have been 
measured. However, depending upon the microwave cir-. 
cuitry, and functions involved, the diodes can operate 


over a wider frequency range (see Table 8.1). 
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Table 8.1 Specifications of HP 5082-0437 and HP 5082- 


0431 Impatt diodes used in experimental 


work [63] 


Diode parameters 


frequency range 

typical operating voltage 
ty picabpopenating acunnrent 
breakdown voltage (Vy) 
junction capacitance 
output power 

thermal resistance 

ae eps capacitance 


package inductance 


5-9. GHZ 
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To measure the CW frequency and power, the mea- 
surement set up shown in Fig. 8.1 is employed. The 
Current cource shown inekig. Gal tssvused in Order to 
maintain a constant current in the biasing circuit. 
Additional dc supply is connected in series with the 
constant current source because the latter has a maximum 
output voltage of 100 volts while the diode operating 
voltage is over 110 volts. For the steady state measure- 
ments, the Impatt diode is first reverse biased into 
avalanche breakdown with Tac = 25 mA and then left for 
about half an hour to obtain a stabilized condition. 

The value of this current is chosen to be the same as 
that for which the diode large-signal admittance has 
been measured. For each diode the steady state (CW) 


frequency and power are recorded and compared with those 


evaluated in the previous chapter, i.e. 


WO = wt x (Se) 


(8.22) 


2) 

H 
1S) oe 
) 
> 


Under CW conditions, both x and A are time independent 
quantities. The measured and theoretically evaluated CW 
power and the frequency of two oscillators are given in 


Table: 6:2 Zi. 
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Table 8.2. Theoretical and experimental steady state 
solutions (frequency and power) of the two 


Impatt diode oscillators used. 


CW operating CW operating 
frequency (GHz power (mW) 


5.8475 Delo, Ol 120 134 
Die Oo 5.804 106 ie Ona, 


The evaluated and measured steady state quanti- 


ties compare very well, in the case of the power the 
difference is about 12% and in the case of frequency 
the difference is only about 0.5%. Especially for 
diode #2, the measured and evaluated quantities are 
very close. Since both the cw power and frequency can 
BE measured with an accuracy better than a few percent, 
the above differences are mainly due to the inaccuracies 
associated with the admittance measurement, which are 
about 12%. Generally, in different cavities the oscil- 
lating diodes can see a little different environment at 
harmonic frequencies compared to that under which the 
diode admittanceswere measured. This fact eentecotn as 
for different inaccuracies characterizing the diodes 
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In order to compare the theoretical with the 
measured transients, the oscillator comprising diode 
#2 is considered under pulsed operation. In this case, 
contrary to the CW operation, the p-n junction 
temperature does not remain constant. The effect 
of junction temperature variation can be minimized by 
bias modulating the oscillator with pulses of high 
repetition frequency. With high pulse repetition fre- 
quency, on the other hand, the characteristics of the 
modulated RF will be dependent on the preceeding modul- 
ating pulses. This may not simulate well the conditions 
under which the transients of the oscillator are evaluated 
theoretically. However, decreasing the duty cycle by 
using pulses of smaller duration (about 30 ns) and thus 
increasing the off time, does not result in a noticeable 
change in the characteristics of the leading edge of the 
RF envelopes. Therefore a modulating pulse with an 
amplitude = 1.4 V, period = 30 ns with 10% ON time is 
ieee Fig. 8.5 shows the wave shapes of the modulating 
pulse, diode current and the RF outputs of the two 
detectors of Hola e ato The discriminator and 


microwave components inserted between the diode and the 


discriminator cause about 10 ns delay between the RF output 


Of sche (OSC lator sand) thew kr lcutput of thesdiscriminator. 
According to this figure, -theiapplied bias pulse 1s not 


completely constant during its flat portion and diode 


189 


abctb gntersande 5 


eoego azdd oF eee” 


jie ae Ronda 


rostonee md ony, nohse 


soohts edt _snenenee! noi 
yd Rhos imiteia ol righ. notary @ et: ' 
tpi to eselug Hiv 4036) 


«{yubom 


are Oi 


of 


tf Loe> 


Si] 


b> 


Pe, 


a) 


bf th 2 


Ve sis 


LD 


alee 


5 
a 


aint 


6 
Cl mek 


so¢sahasacetih sat peter mica eat adel 


} a. 


abort 


w3 ect Yo eye sete “te ei? bos a altel seers 


22 anda iasqes oalag pat (220 yoteupess 


evi.) 


toveiiiose dh Se esa is tanend ett 


7 


no’ @), As iw erp ets hot29q * dite 


sik 2d siiteiatincd avav etlt on b's & Pi) a 


y ’ 
ee ia ; 
mA T4 ve 


' : : : oy rs 


Ve 
stostate st? paw seidso 942.9. Rey 
anergy sth HS snebreash od ELiw wit | | 
‘tow etelamte: don gay eats 


1h ads yoless Ya2t nattowoity » 500 
i - - 
Of (awodey noiteareb iba 


Pt. J4Avers  20On eaaob a! sits 4 ato a2) © ni > ne] 
: A Lue " 
“ibsel sit Ae spideiresosseria: etm adi pared 


=f 


caine priseivbog s sactoged? ) «aes 


ed 


‘) One sheib,odd meewseg betzsani ed conogme: 
civ Gf .eta aeawced valet an 04 od aeyse: motes 
oath SCL a, Tey ENS es. arty baa sors: 
ovluq end Geliqge on sotepk? liek oo 


sboib Bas aobewog MER 624, pried ounbede 


7 


_— ye A ie 


190 


Fig. 8.5. From left to right; the modulating pulse 
(500 mV/cm), the diode current (100 mV/cm) 
and the two RF envelopes recorded by the 
X- and Y-detectors (20 mV/cm), respectively, 
for diode #2. Modulating pulse amplitude 
= 1.4 V, period = 300 ns with 10% ON time 
pulse. 
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current dOesnolexactly follow the bias pulse (Fig. 8.6). 
Small parasitic reactive elements in the current 
Measuring circuitry together with junction temperature 
effects may be responsible for small oscillatory fluc- 
tuations on the current waveform at its leading edge. 
The oscillation growth is taking place under almost 
isothermal conditions, note that the diode current before 
the pulse leading edge is constant. The oscillation decay, 
after Ghe™pulse is turned off, is significantly affected 
by the variable junction temperature. As can be seen, 
the diode current slowly decreases long after the bias 
pulse trailing edge and, because the junction temperature 
has strongly affected the diode RF characteristics, the 
trailing edge of the measured RF pulse reflects the 
transients caused by sudden change in the diode negative 
conductance which tssstullylarge enough to compensate 
for the circuit losses. Note that a relatively rapid 
transition of a trailing edge of the RF pulse is followed 
by a eae varying portion which characterizes the decay 
Of the small residual oscillation directly controlled® by 
the junction temperature of the diode. 

The two RF envelopes of Fig. 8.5 together with 
the “calibration curves in fig. 6.7 are used to Avabbetee 
the rise time of the oscillations. This has been 
achieved by calculating the instantaneous power at 


various time instants starting from the leading edge 
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Fig. "ots “A. reproduced "copy of Figuecvo with the timescale 
adjusted for a 10ns delay caused by the 
Measuremenc System 
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of either of the RF envelopes. The values are then 
plotted on the calibration curves of the system (x-y 
Dover Hsp love ioe 7 Jeane Dy Orerog lat Lom, othe 
instantaneous frequency and power are determined for 
each time instant. The points numbered 1,2,..., are 
corresponding to different time intervals along 
succesSive points on the RF envelope graphs. 

Fug. 8.8 shows the®growth ofthe oscillator RE 
power vs time as obtained from the two experimentally 
Gecermimed RP “envelopes an Taig.wiseosatterumul t ipd ying 
Dbyta Tactor "Or two tomadlow for thes 3) dbycoupker of 
Fig. 8.1. The figure also shows the growth as calcula- 
ted on the basis of the backward recursive method 
(Fig. 7.8). The measured power is shifted in time by 
8 ns in comparison to the time scale for the two RF 
envelopes (Fig. 8.5). The experimental curve indicates 
that the rise time of the RF oscillations is about 7 ns. 
In view of the discussion of the transients in the bias 
Circuitry in sectiongo.|l, this valtietcompares*tairly 
well with the value of 5.4 ns calculated on the basis of 
the backward recursive method (Fig. 7.8). It can also 
be noted that the measured power corresponding to the 
flat portion of the bias pulse is about 126 mW as | 
compared (to the calculated steady state value of 116.7 
mW. This difference can be attributed to the fact 


that the calculations are based on the diode admittance 
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which is measured under constant thermal conditions. 
However, in the pulsed mode of operation (due to the 
ON and OFF sequence of pulses) the diode junction tem- 
perature is slightly different from that at which the 
large signal admittance is measured. 

Fig. 8.9 shows the instantaneous frequency vs time 
behaviour (the crosses indicate measured values deduced 
from the tEwo IRF @envelones, of fig. 3.5 ,., and, the-dots 
show figures calculated from the theoretically deter- 
MineasvaluecswoLe. ancdeae(hig. 7. 6) mit mavealsowbe 
noted that in this case the time scale for the measured 
frequency is shifted by 8 ns in comparison to the time 
scale for the two RF envelopes. The two curves in Fig. 
$7.9 show a Similar trend. until about 33 ns. The frequency 
difference of about 80 MHz is observed at the leading 
edge of the bias pulse. For the flat portion of the 
pulse, a frequency difference of only about 10 MHz is 
observed. 

The Impatt diode oscillator used in these measure- 
ments is designed to oscillate at about 5.8 GHz if the 
Giode bias current is 25 mA. To avoid the RF pulse 
jitter and to operate at junction temperature close to 
that at which the large signal admittance of the diode 
was measured, the diode is pulsed from about 16 mA (OFF 
state) to 25 mA (ON state). The bias current reaches a 


value of 25 mA, and thus the design conditions are 
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Fig. 8.9 Measured and calculated time dependence 
of the frequency of the Impatt diode 
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aecaiioawauring the flat? portion of tire purse.” Inthe 
OFF state, no oscillations should be present in the 
circuit. When taking measurements, the output of the 
discriminator diodes was used to determine the OFF 
Pulses cumrente evel “Fig. 6. ce noOwc omOscCi | ati Ons 

in the OFF state. However, the diode sensitivity is 
very limited and any low power oscillations can hardly 
be detected. A spectrum analyzer should be used when 
selecting the OFF state current level. Fig. 8.9 
indicates that low power oscillations were present 

in the’ oscillator during the OFF state. The frequency 
of these oscillations, determined by the circuit 
impedance and the small signal impedance of the diode 
bias at 16 mA, seems to be around 6.0 GHz. When the 
diode is pulsed on, the oscillations grow from the low 
level oscillations, already present in the PaSenetee 
and not from the circuit noise whose frequency is deter- 
. Mined by the circuit impedance and the small signal 
impedance of the diode biased to 25 mA. Thus, at the 
beginning of the ON pulse, the frequency of the experi- 
mental oscillator changes rapidly and the conditions 
are very different from these used in the calculations. 
When approaching the flat portion of the pulse, the 
experimental conditions approach those of the conditions 


ysed, Inthe calculation stints, 16s Gerue for both. the 
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frequency and the rate of change of frequency. The 
undesired oscillations may be avoided if the oscilla- 
tOr 1S caretully tuned, not only for operation at 25 
MA OUP aLSOe LOL 16 mA. 

The time dependence of the oscillator frequency 
during and after the trailing edge of the bias pulse, 
i.e. for instants of time exceeding 33 ns, is deter- 
mined by the residual oscillations. Note that the 
measured RF envelopes indicate oscillations far into 
the OFF state. These oscillations could be prevented 
if lower bias current for the OFF state is selected or 


if the circuit 1s properly tuned, 


8.3 Bias modulated Impatt. diode oscillators 


Bias modulated Impatt diode oscillators can be 
used in transmitters of digital radio systems. These 
combine modulation and generation functions in a single 
stage and depending on output power requirements, RF 
amplification stages may not be needed. Hence such 
systems may be economically very attractive. The 
noise disadvantage of ASK modulation in such simple 
systems can be partially offset by the fact that 


higher pulsed power can be obtained from a given 
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osellilator: If a 50% modulation scheme is adopted, 
the oscillator pulsed output power can be almost 
doubled in comparison with its CW outputs. 

The spectral distribution of the energy of the 
transmitted signal is an important consideration ina 
modulation scheme. In an ideal ASK system, the RF 


spectrum, centered at the carrier frequency, is iden- 


tical in shape to the baseband spectrum of the modulat- 


ing signal. In this case, the percentage of the power 


in the main lobe of the RF spectrum is the same as the 


percentage of the power in the main lobe of the baseband 


spectrumpri tee, Othe *spectrum of ithe envelopey +For®an 


ideal ASK system, the RF spectrum is always symmetrical 


around the carrier frequency, regardless of the pulse 
shape. 

In a digitally modulated Impatt oscillator, the 
carrier frequency is not constant and an undesired FM 
component is present. This is due to the fact that 
the electronic susceptance is a function of the diode 
current, the amplitude of RF voltage, and the junction 
temperature. Hence when the current changes, the 
frequency will also change. Under pulsed conditions, 
the shape of the current waveform and that of the 


applied pulse are similar, resulting in the frequency 


being decreased during the build up of the oscillations. 


This observation can be verified from the change of the 
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instantaneous frequency during the pulse ON period as 
shown in Fig. 8.9. The frequency also changes during 
the main part of the pulse due to heating of the diode; 
this is a consequence of the temperature dependence of 
breakdown voltage as mentioned earlier. 

To give some insight into the limitations of 
bias modulated Impatt diode oscillators two sets of 
experiments are performed. The diode is dc biased at 
the threshold of oscillations and then modulated with 
pulses of either varying duty cycles or 50% duty cycle 
with different periods (Tables 8.3 and 8.4). 

For the operating conditions described in Tables 
8.3 and 8.4 the two RF envelopes detected by the X- and 
Y-detectors, and the diode current are recorded simul- 
taneously. The results corresponding to the duty cycles 
Of, 05. and 0.2 "arewimuustratved si twig See.o. | OmandpG clean. 
These figures show that the shapes of the RF envelopes 
are different and also the decay time increases with 
decreasing duty cycle. The latter observation is a 
consequence of temperature effects. The two instants 
at which maximum power and .87 maximum power occur 
are used as reference points for which the corresponding 
change in frequency is evaluated. These power points 
correspond approximately to the main portion of the RF 
envelopes. recorded for various: duty cycles. Fig. 8.12 


shows the measured effect of duty cycle variation on 
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From left to right; the diode current 
(100 mV/cm) and the two RF envelopes 
(20 mV/cm) recorded by the X- and Y- 
detectors,respectively, under the 
operating conditions stated in Table 
8.3 for a modulating pulse with duty 
cveLweumre 0.5. 
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From left to right; diode current 
(100 mV/cm) and the two RF envelopes 
(20 mV/com) recorded by the X- and Y- 
detectors, respectively, under the 
operating conditions stated in Table 
8.3 for a modulating pulse with duty 
cycler. U.2. 
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Table 8.3. Operating, conditions under changing duty 


cycle from 502 ‘to; 102. 


Characteristic 


Pulse duration 


DC bias current 


Applied voltage pulse amplitude] 2.1 V 


Table 8.4. Operating conditions under fixed (50%) duty 
cycle and changing pulse period from 60 ns 


tOeS0U0 Uns. 


Characteristic 


Pulse duration 


DC bias current 


Applied voltage pulse amplitude] 2.1 V 
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frequency during the main portion of the pulse. The 
dashed line for a given period reflects the frequency 
change between the two reference points in the main 
portion of the pulse. It is seen that, for a given 
pulse length, the frequency variation during the ON time 
increases as the duty cycle decreases. A similar 
behaviour is observed for the decay time of the RF 
envelopes which increases with the decrease in the 
duty cycle (Fig. 8.13). 

Below an OFF time of about 240 ns, the average 
diode junction temperature increases and there is a 
smaller temperature variation during the ON time. On 
the other hand keeping a 50% duty cycle and increasing 
the pulse period does not show any significant change 
in frequency because the temperature variation between 
ON and OFF periods of the pulse is smaller as compared 
to the smaller duty cycle pulses. A reference point is 
chosen at the maximum power level. The instantaneous 
frequency at this point is calculated and plotted against 
the corresponding pulse period (Fig. 8.14). It is 
observed that there is insignificant change in the 
instantaneous frequency with the pulse period employed. 
It was also noted that there was no significant change 
in the decay time. 

In realistic digital applications, the change of 


the carrier frequency is undesirable since it will 


: ‘ 

oft. .sedug an orn 
yorleaue att ‘expalts A 

geek a 

nian sty At 236: oe som ee ox" 


4 tot , dad + ape e at ee 


<y5 F- 


— 


it WO ety BRASCN F itg.iany¥ eels 
tmitawe . roan eta ove ¥ 
ia ae Ve nak Pye x02 


7.) » 
oh ae 


_—" sire 
Lan ae, oe “8 fo. 


ea 
“at ae, trodes ‘Io smig 190 ne work ty 


i] 5, 
l sage tons’ tae rogiios fx. shed seborek 


4 ‘ t¢ Atiw ascites 


= 


= « | iio MsSEL ASV = sqmet at 
paki, ‘bher ales 


Se! 
’ 
ey 
&32 
- 
— 
- 


wone Boa 2 meley ‘Bota hs ‘98 mpi 
») : : \ 
2) ; : ects ie 
ifetuay etuirdotsgmes ns yarsood ‘yore gest - 
: : Pec (4 By abe 
ering paioy> 4 
i P eae pt s2Wwoq mii 
Js6LuGI GE 2s LOG 
- rs 7 Fa ; < = 
pit) beitteq s2l:4q patbnagesiad 
: =e ber ¥ 
tY efit tngoiltepkes; 2. axed? shaied wane 
. ane Hailing oalua edt sy 2coupens avowtts 7 


“y 


eve 49-JneoMitne.s Off aaw aes TSH besos oot mm 5 


ar 


4 . “ - a — 
Ef 4 pons et sehen 2 s 


5.88 


7.86 


oo 


INSTANTANEOUS FREQUENCY CHANGE (GHz) 


§ 82 


Oat 0.2 Gus 0.4 0.5 
IDIOMINS (GNACIiaa 


Fig.8312 Frequency variation, for pulses of 30ns ON time, as 
a function of the duty cycle between two instances 
for which the power changes between its maximum 
value and 87% of this value. 
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Fig. 8.13 Decay time of the RF envelopes as a function of the 
aduty cycle due sto a, 30ns ON pulse. 
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INSTANTANEOUS FREQUENCY AT MAXIMUM 
POWER LEVEL OF THE RF ENVELOPES 


(GHz) 
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Fig. 8.14 Instantaneous frequency at maximum power 
level vs. pulse periods (50% ON time pulses). 
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SCALE : e8ns7daiv-: 


Modulation scheme for utilizing 
bias modulated Impatt diode 
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result in intolerable interference in adjacent channels. 
The experimental results shown above indicate that the 
most serious drawbacks of the simple ASK bias modulation 
associated with the changes in the carrier frequency 

can be minimized by employing coding schemes based on 
Dus Gury cycle pulses, =According to rig. 6.13, the 
narrowest pulse which can be handled by the bias modu- 
lated Impatt diode oscillator is 8 ns. Adopting a 
modulation scheme where binary zero is characterized by 
8 ns pulses and 50% duty cycle, about 40 M bits pulse 
stream (see Fig. 8.15) can be handled by the bias 
modulated Impatt diode oscillators without introducing 
unwanted interference. To decode such a signal, the 
decision time instant should be delayed behind the pulse 
beak by about 3 nS as Shown in fag. (8. tS: 

Since only a simple peak detector and a constant 
delay circuit are needed on the receiving side of the 
system, the above modulation scheme makes the bias 
modulated Impatt oscillators suitable for digital appli- 


cations handling bit rates up to 40 Mb/s. 
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CHAPTER IX 
SUMMARY AND CONCLUSIONS 


In this study a nonlinear model has been 
developed for the investigation of microwave circuits 
with embedded two terminal active devices. The general 
model is equally well suited for both discretely charac- 
terized and analytically described active and passive 
Circuits. It can be used in the design, study and 
optimization of microwave oscillators and amplifiers. 

By separating the passive microwave circuitry from 
the active device, the active microwave circuit can be 
well described by two nonlinear differential equations. 
Unlike the previous models, which take into account the 
RF voltage dependence of the active device admittance 
up to a maximum of the third order (e.g. Van der Pol's 
cubic nonlinearity), the model presented here considers 
the admittance dependence on both RF voltage and fre- 
quency up to the sixth order. Moreover, the previous 
models, based on single tuned GLC circuits, were res- 
tricted in their applicability to simple microwave 
circuits for narrow band operation. The model presented 
in this study, on the other hand, is of a more general 
nature which makes it suitable for any dynamic and steady 
state study of broad band oscillators and injection 
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has been illustrated with a transient and steady state 
StuaveOLe thea lmpattadiodesosei Maton. 

Numerical methods have been developed for the 
solution of nonlinear equations describing steady state 
and transient operation of Impatt diode oscillators. 

It has been found that the transient study based on a 
forward recursive method is suitable only for simple 
GLC circuits while the backward recursive method is 
found to be, so far, the only alternative for obtaining 
a solution for the more complex circuits encountered 

in? pracei ce: 

The determination of the model parameters requires 
the large-signal admittance values of the diode and also 
the admittance of the embedding circuitry. Therefore, 
an accurate method has been designed and tested for the 
measurement of both the large-signal admittance of 
an Impatt diode and the embedding circuitry. The mea- 
surement system is characterized by two 2-port networks 
whose parameters are determined with an error of less 
than 3%. The method is also readily adopted for measure- 
ments at various frequencies. One special feature of 
the measurement technigue is that it effectively 
eliminates the systematic errors inherent in any micro- 
wave measuring system. Two types of Impatt diodes were 
operated as reflection type amplifiers in a tunable 


coaxial resonator. The large-signal admittance of these 
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diodes and their associated circuitry, when embedded in 
a specially designed oscillator, were measured by the 
above mentioned technique at different power levels and 
atra single bias current’ of "25 mA. fhe overall! error 
in these measurements is estimated to be in the neigh- 
bourhood’ of 12%. This “accuracy is better than that 
obtained by all other methods except that of Van Iperen 
and Tjassens. These authors use a more complicated 
approach which limits the applicability of their system 
if measurements are sought at various frequencies. 

A coaxial type Impatt diode oscillator with a 
Single slug transformer has been designed. Special 
Dias circuitry optimizes the oscillator characteristics 
under pulsed operation. The rise time of the diode 
current pulse is about 3 ns or less which compares very 
well with the 2.0 ns rise time of the applied voltage 
pulse. The oscillator is capable of handling pulsed 
modulating signals up to 40 Mb/s. 

The two Impatt diodes characterized above were 
end mounted in two single slug transformers to form 
two oscillators of the developed type. The steady state 
power and frequency together with the transient response 
of the second diode were measured with the help of a 
frequency discriminator. 

In the course of this study, a novel frequency 


discriminator has been developed to measure the Impatt 
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diode transient and steady state behaviour. Its sim- 
plicity allows the discriminator to measure power and 
frequency changes in nanosecond time intervals. The 
discriminator can be calibrated using CW signals. Using 
a PIN modulator, the discriminator error in measuring 
the frequency of fast changing signals has been found 

to be better than 4 MHz, 

The measured steady state frequency and power 
together with the transient response of the two Impatt 
diode oscillators have been compared with the computed 
solutions obtained on the basis of the backward recur- 
sive method. Fair agreement between the theoretical and 
experimental data has been found. The steady state 
power differs by 12% and frequency by 0.5%. On the other 
hand, the measured rise time of 7 ns is about 1.6 ns 
longer than that evaluated theoretically. 

The backward recursive method has also been used 
to evaluate the response of the two hypothetical diodes 
for which the large-signal admittance has been published. 
These diodes are considered to be embedded in two 
resonators with different sets of parameters. It has 
been shown that although the two oscillators are tuned 
togtbe same. frequency ol oscillation, their calculated 
rise times differs, by about 203. .In addition, the: follow-— 
ing characteristics, havewbeon revealed Dby~the studies ‘on 


Impatt diode oscillators under steady state and transient 
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conditions: 
de The diode admittance dependence on RF voltage is 
not as simple as that expected by Van der Pol's cubic 
nonlinearity. 
7 ag The nonlinear dependence of the diode admittance 
on frequency cannot be ignored. 
Sn The calculated steady state and transient para- 
meters of the oscillators determined by using the present 
model depend on the accuracy by which the admittances 
of the diodes as well as the passive circuitry is 
measured. 
4, The value of the oscillator transients depends on 
the interaction between the device characteristics and 
Ener Cy CUlteELOCI . 

During the course of this work, certain aspects 
of the study arose which merit further investigation: 
a) The influence of the interaction between the 
device line and circuit locus in optimizing the transients 
Mieosoi Ulators. 
b) The stability anvoscillatorssand: amplifiers’. 
ec) Incorporation of the presented model into avail- 
able computer aided routines for design analysis and 
optimization of microwave, active circuits. 

The model developed here can Li be used in 
dynamic studies on broad band power amplifiers. The 


presented model has been successfully used in evaluating 
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distortion and intermodulation products in Impatt diode 
amplifiers [65]. 

The discriminator developed for the measurement 
of instantaneous power and frequency can be adapted, 
due torics simplicity and very fast sesponse time,, -to 


studying the transients in digital communication systems. 
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Apperdix I 


LARGE-SIGNAL ADMITTANCE OF IMPATT DIODES 


C THIS PROGRAMME EVALUATES THE ABCD PARAMETERS OF THE NETWORK BETWEEN THE 
MEASUREMENT REFERENCE PLANE AND THE ACTIVE IMPATT DIODE CHIP 

CE HK Xe He HE HE KE A OK Ok 

C THE USER MUST SPECIFY VALUES FOR THE FOLLOWING: 

AEPSsTHE ABSOLUTE DIELECTRIC CONSTANT (PF/CH) 


e 


Cc 


Cc 
- 
Cc 
e 
Cc 
Cc 
Cc 
Cc 
C 
6: 


=1.036 FOR SILICON 
=1.416 FOR GERMANION 


CONK:A CONSTANT 


=2.38D+17 FOR SILICON 
=1.74D+17 FOR GERHANIOM 


AJ: THE DIODE JONCTION AREA (CIRCULAR) (SQUARE CH) 
20: THE CHARACTERISTIC IMPEDANCE OF THE MEASURING SYSTEM (OHMS) 
NN: NUMBER OF MEASUREMENTS (20 MAX ) 


ES 
AeA We I a ae ai Hee eR oe IE 2k ak 3 ak 


FREQUENCY OF HICROWAVE IMPEDANCE DATA (GHZ) 


IMPLICIT REAL*8 (A-H,O-Z) 

REAL*8 LAMDC,LAMDG,LAMNDO 

COMPLEX*16 CHZ0,ZCR,YCR 

COMPLEX*16 GM,GE,ONE,GCAL,GAM (20) ,GAC (20) 

COMPLEX*16 DUMM, ¥D(20) , ZF (20) ,DET,AC,BC,CC,A,B 

COMPLEX*16 DELZD,DELZU,2ZN,AN,BN,CN,AF,BF,CF,DP,REFLC, DUMMY,DF2 , 
1 CHECK, ZD,ZMC (20) 

DIMENSION DELZD (20) , DELZN (20) ,PIF (20) , DIST (20) ,PIQ(20) 

DIMENSION V(20) ,CD(20),SHL(20), 

ICDAV (20) ,CONC (20) , RHO (2G) , DELRD (20) , DELXD (20) , DELCD (20) ,X (6), 


3 


FX (20) ,RD (20) ,XD(20) ,RM(20) , XM (20) , THEFF (20) 


DIMENSION RLF (20) ,THETF (20) ,PZN(20) ,VAC(20),R(20) , THEY (20) , RFF (20) 
COMMON/FUNC/ZM (20) ,Z2D(20),NN,NC 

COMMON/ALPHA/A (3,3) ,B(3,1) ,GM(20) ,GE(20) ,GCAL (20) , DET (20), 

*RL (20) , THETA (20) ,GAMNA (20) , THETAR (20) , WA (60) , PL 


DATA 
DATA 
DATA 
DATA 
Sid 
DATA 
DATA 


SHTIZ0LI0 At eR SIS, 535 a PS, 25, 55,206,265, 780007 
CL, LANDC f2e998D10,6.97D0/ 
AEPS,CONK,Ad,Z0/1.036,2.38D+17, 1. 0D-04,50.0/ 

CDVSUGL O30 10228670. 5, 6725682, 38, .35,+320, 2906, «289 


W201 6 95). 2 204 830.5) FORNS0 GOs TOLS0L 90. 100.7 
PIOQ/-20. Paces | US's Ss HOR gop c ples ll selec ge he plae/, 


EXTERNAL UOBJ 
Coksk sok tokio dk KAT NP OT DATA ROUTINE % 2% 4K 


NN=13 
NC=13 


14159265359D0 


DO 444 TII=1,2 
READ(5,441) F 
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NNM1=NN-1 

441 FORMAT (F10.5) 

LAMDO=CL/ (P*1.D9) 
LAMDG=LAMDO/DSORT (1-( (LANDO/LAMDC) **2) ) 
LANDG=LAMDG/2.54 

READ (5,4) (RL (Z) ,I=1, NC) 

READ (5,4) (THETA (I) ,L=1,NC) 

CALL ANCALI (X,LANDG,F,SHL) 

WRITE (6,940) 

940 FORMAT (1H0,3X,'MOVABLE',11X, "REFLECTION! , 22X, *REFLECTION',20X, 
x *REFLECTION',13X,'DIF IN CALC REPL. COEFP"'/4X,*SHORT',13X,'COEFFIC 
*¥LTENT®,21X,'COEFFICIENT', 19X, 'COBFFICIENT'/3X,"POSITION', 7X,? 
*THEORETICAL',21X, *MEASURED*,21X,* CALCULATED’) 

DOPE V3e1=1, NC ; 
WRITE (6,944) SHL(I),GE(I) ,GM(I) ,GCAL(I) ,DET(I) 

613 CONTINOGE 

944 FORMAT (2X,F7.3,3X, 1P8D15.5) 

NC=NN 
AC=DCHPLX (X (1) ,X(2)) 
BC=DCMPLX (X (3) ,X(4)) 
CC=DCNPLX (X (5) ,X(6)) 
ONE=DCNPLX (1.D0,0.D0) 
READ (5,4) (RL(I) ,I=1, NN) 
READ (5,4) (THETA (ZI) ,I=1, NN) 
WRITE (6,557) (RL (I) ,I=1, NN) 
WRITE (6,558) (THETA (2) ,I=1,NN) 
557 FORGAT (1HO,*RL="', 13 (F8.2,1X)) 
558 FORNAT(1HO, *THETA=',13(F8.2,1X)) 
4 FORMAT (8F10.4) 
DO 5 T=1,NN 
RL (I) =RL (I) /50.0 
THETA (1) =THETA (I) /10.0 
GAMMA (I) =10. ** (RL (I) /20.) 
THETAR (1) =THETA (I) *PI/180.0 
GCOS=GAMMA (I) *DCOS (THETAR (I) ) 
GSIN=GANMA (1) *DSIN (THETAR (I) ) 
GAM (I) =DCNPLX (GCOS, GSIN) 
GAC (1) = (GAM (I) -BC) / (AC- (CC*GAM (I) )) 
ZMC (I) = (ONE+GAN (I) ) / (CNE-GAM (I) ) 
ZM (1) = (ONE+GAC (I) ) / (ONE-GAC (I) ) 
5 CONTINGE 
8 C#R RAM AKRREKYRITE THE INPUT DATA *%*%% WK 

21 WRITE (6,22) NN,F,Ad 

22 FORMAT (1H1,10X,"NUMBER OF MEASUREMENTS: ',13/10X, "FREQUENCY OF MEAS 
JUREMENTS:',1X,F10.5,2X,"GHZ'*/i0X, ‘AREA OF DIODE JUNCTION: ',1X, 
2 1PD15.5,2X,*SQUARE CHt/) 

23 WRITE (6,24) 

24 FORMAT (1HO,4X,*BIAS*,2X,'DIODE CAPACITANCE',2X,"*NORMALIZED IMPEDAN 
1CE MEASUREMENT! ,5X, "NORMALIZED IMPEDANCE CALCULATED!,13X,'*REFLECTI 
SON*/ 3X, "VOLTAGE',4X,"AY 1 MHZ",16X,tAT FREQUENCY F',22X,"USING 
GANA-ABCt,21X, 'COEFFICIENT' /34X, "RESISTANCE", 8X, !REACTANCE', 8X, 
** RESISTANCE’, 8X,*REACTANCE', 8X,"MAGNITUDE', 6X,*PHASE'/3X,' (VOLTS 
*) ©,6X," (PF) *,100X,* (DEG) *) 
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DO 26 I=1,NN 
WRITE (6,25) V(I) ,CD(Z) , ZMC (I), ZN (I), GAMMA (ZI) , THETA (Z) 
25 FORMAT (2F10.3,5X, 1PD20.7, 1PD17.7, 12D20.7, 1PD17.7,D16.6,D13.3) 
26 CONTINUE 
CK HE 3K He Ae AE Be ee Hee ie A OK ak aK aE 3K ae OK ake a aK OK ae ah a ok Ye 
C CALCULATE IMPURITY CONCENTRATION FROM CAPACITANCE DATA 
CEM AE Me Be Hee OK De a AE Ae Ke He 6 Ke OK ake ak HK ah ae a aK ak OK 2 afc 2k 3K oe 
DIAM=DSQRT (AJ/PI) 2000. /2.54 
DO 32 I=1,NNN1 
CDAV (1) =0.5* (CD (I) +CD (I+1)) 
DELCD (I) =CD (I+1)-CD(I) 
DELV=V (I+14) -V (I) 
CONC (I) =- (CONK*DELV*CDAV (I) *3) / (DELCD (I) *DIAN**4) 
32 CONTINUE 
CE BE AK ee ae 3K a he Me 3 of aie a ke Ae OE oe 2h Hee ae fe 23e 9fe Yo at 9K Ke 
Cc CALCULATE RESISTIVITY 
Cc FOR N-TYPE SILICON: 
C 2K FE ie AE BK Oe KE KS Of aK as ap ae at oh fe ae ais ae ae oie fe 3K he Ne ie 
DO 34 I=1,NNN1 
RHO (1) =DEXP (28. 3637-. 7861*DLOG (CONC (I) )) 
34 CONTINUE 
CH HK A BENE Be He OK Fe ok Oe ie ie ke ae ae ae ee Be aK ah OK Oe Oi 3K OK oe 
C WRITE CONCENTRATIONS AND RESISTIVITIES IN YHE VOLTAGE INTERVALS 
CE AK DK AE 2 OE Fs ae AE Oe Se afc af 3s Dz ale Ae HE SHC TK Bie Ie DK Ic Acakx ec 3k 
WRITE (6,41) 
41 FORMAT (1H1,3X,'VOLTAGE', 11X, ‘IMPURITY! ,10X,"RESISTIVITY'/3X, 
1 "INTERVAL', 8X, ‘CONCENTRATION'/18X,'(PER CUBIC CH)',10X, 
2 ' (OHM-CH) t//) 
DO 43 I=1,NNM1 
WRITE (6,42) 1, CONC (I) , RHO (I) 
42 PORMAT (5X,13,4X,1P2E20.7) 
43 CONTINUE 
CttttrtTrrCtrrrrcLeeerrererrres ££. fe te 2 2 2. 
c CALCULATE CHANGES IN 1 MHZ AND VHF IMPEDANCE DATA 
c IN THE VOLTAGE INTERVALS 
COR I a Oo a ak te ak a) ok II IO I aa I ak ak ak ae ai ak 2k aie ak ok ake ak 2k ok ak ak ak ak ake 
DO 51 I=1,NNM1 
DELRD (1) =-AEPS*RHO (I) * (1. /CD (I+1) -1./CD (I) ) 
DELXD (1) = (500. *DELCD (i) ) / (PI*F#CD (I) ¥CD(I+1}) 
DELXD (1) =~ 1. *DELXD (I) 
DELZD (L) =DCNPLX (DELRD (I) , DELXD(I) ) 
DELZti (I) =ZM (I+ 1) -2M (1) 
51 CONTINOE 
CR He Fe He Se 3 Re ea Ok a aK ik ae Fi ie Tk aie aR aR it a AE aK Ks 9K 2 C2 2K 2 3h 2 2K 2K OK 
C WRITE THE CHANGES IN THE 1 MHZ AND VHF IMPEDANCE DATA 
c IN THE CORRESPONDING VOLTAGE INTERVALS 
CK He He Eee OK ae OE ote aR a 2 A RO Fe a OK I eK I ae A I IK Ae Oe He Ne I ie HIKE He AE IE OK 
WRITE (6,61) 
61 FORMAT (110, 3X," VOLTAGE',9X, "CHANGES IN 1 MHZ IMPEDANCE DATA, 
1 7X,"CHANGES IN MICROWAVE IMPEDANCE DATA'/3X, "INTERVAL"/19X, 
3 *RESISTANCE',11X, *REACTANCE', 10X, "RESISTANCE, 11X, * REACTANCE'/21X 
3,* (OHMS) ',14X,* (OHMS) ', 14X," (OHMS) ',14X,* (OHMS) '//) 
DO 64 I=1,NNM1 


Cen Oe i Tas 

“er i i i. Nia 
, , Patan 
j i a. 


; ae ha - a? 
| (3) ATENT STU) snag ¢t ity w (EpoRe, 
(4,0. 0%, F. Fae? .o.0SaaR, ronal 


coea i 


ig .t\. ongte 


PUES hia | my wre = 
(=) a 


(ee Tyas teaya cx ity waaay eae OBy Ke. 


jee CUE 7 
(sp eeheaines sain Ea 
cvTrexemn Wa WITAD: ae 
. | _ SOSER oe 
siete Srrrets hanes — 
9 ‘vo nut ur 0 
is C88 Bs) 110 . Pao 


40 SIM ATION KON vous 


Tepe « Wa 
ee in ied widanennen <p 
(CETVETESE RH OMe pene hii np jag be ina 
ee elie Dies dnansie we he Rah 
. (a apatem 
= Use, fT suv" tes (HPY EA eben 


ot J nical ; 
' py 1 i Ada” a, * Tava ot y 
as ANS atetito ry , 
. trypan, (1) 0% 
boa Me 
- + See wwe bo en 
a a" ! “GHA 2GA PF ee S306 US Ts 
ery GwaATaeov TAT 4 Pea 
was com 6 ew ee mee: sh sekebibhtememell 


Tat ted papa 
+L OOS wh) (Lp ORR et ag, a 
reve Tay i i or, iZuF « Ook) = (4) asd ‘ + : 

., 7 es, < ey okeede Ft Fy OKs: eS 
{ (rp@xuac, tz) as JES) KL MOGs Ly 
oS (PATy Mahe nas 


bd eaG itd bay oy iiaduiisguhinn 

0 SWAGHIAL GV Abs 2M 1 ART WS soolano' ave 

RIAY sx7es Beat PET mi TAN ORREEROD H4 ; 

weds et ee ee ; 

ele Se Seman 

Sun y he Re sha! te , torino Ey Re jour pray 

ttn ated eeandanT IVAPORD TA ME SHOUARD* 2 

A* by ete a iene". yor,” JoRAT: a u*. m3 Tota’ ene 
\* (TOD * ey ¢Sn40) abr ‘AO ah =f 


+) 
ie | 


WRITE (6,63) 1,DELRD (I) ,DELXD (I) ,DEL2M (Z) 
63 FORMAT (5X,13,4X, 1P4E20.7) 
64 COKTINUE 
CK 3K ae OK Ae Fe KD hs oe Oe Ae ae fe he he ok oO OC ae ok ok oe afk ok ae ok 2 ak ok ake ake a 
Cc CALCULATE ABSOLUTE VALUE OF 1MHZ IMPEDANCE 
RD (KN)=0.0 
C RK AE EK SE Be 2 Ne OE a OR ak Be BE a ak ee 2 3 oe oe aie of ale ak ke hk a ok ae of oe oe ok 
DO 71 I=1,NNM1 
J=NN-I 
RD (J)=RD (J+1)-DELRD (J) 
71 CONTINUE 
DOwI2r = HN 
XD (ZI) =-500. 0/ (PI¥F*CD (I)} 
ZD (I) =DCHPLX (RD (I) , XD (I)) 
RM (I) =DREAL (2M (Z) ) 
XM (I) =DIMAG (Z4 (I) ) 
72 CONTINGE 
CoE MK Ae aR Ae Ne fe ae aE ake OK A AE ake eH OK Be ole ok at 3h 2s 2k 3k afc 2 as aie ok oe ak ok ok ak af ake ok of 2K oe 
C ALL NECESSARY EXPERINENTAL DATA HAS NOW BEEN OBTAINED 
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Cc PROCEED WITH THE OPTINIZATION ROUTINE 
C DEFINE SCALAR VARIABLES REQUIRED BY ZXPOWL 
CK 2 oe Fe a ee OK SK Fe A ae aie He ak AN a OK aK 3K OK ok fe ie 3 3 ake aa Oe ik ais i 2 a ac aK ae 3K 2 
EPS=1.D-12 
ITNAX=10000 
N=NN 
N=6 
PEAK FE AE IK 2k Ms AEE I Ae OP aK ake Be ak Fe ae Es fs Bs OE I oe ie OK Ok is Os Ok OK OK 2 Ok ot ok of ale aK OK 2k 
Cc CALCULATE INITIAL VALUES BY SOLVING SYSTEM OF LINEAR EQUATIONS 
Cc RITH COMPLEX COFFFICIENTS OBTAINED FROM 3 SETS OF IMPEDANCE DATA 
C24 HK AEA Ae AK AK A Be SR A a Ts AS AER aK ie Ss IC OK a Sis ee a aE aie OK ah oe a XK 
AAR=1.0 
RRI-O.0 
DOs = 153 
J=3*T+4 


WO 7D (3) 
A (1,2) =DCHPLX (AAR, AAT) 
A (I, 3) =DCMPLX (XD (J) *XM(J)-RD (J) *RM (J) ,-RD (J) ¥XH (J) —-RY (J) ¥XD (J) ) 
B (I, 1)=Z8 (J) 
73 CONTINUE 
CALL LEQT1C (A,3,3,B,1,3,0,WA, IER) 
X (1) =DREAL (B(1,1)) 
X (2] =DINAG (B (1, 1)) 
X (3) =DREAL(B(2,1)) 
X (4) =DIUAG(B(2,1)) 
X (5) =DREAL (B (3,1) ) 
X (6) =DIMAG (B (3,1) ) 
WRITE (6,74) 
74 FORMAT (///, "THE STARTING ABCD PARAMETERS ARE:'/) 
WRITE (6,82) X(1) ,X (2) 
WRITE (6,83) X (3) ,X (4) 
WRITE (6,84) X(5) ,X (6) 
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FORMAT (1H0,4X%,°A =*,1PD20.7,2%,'+3 ', 1PD20.7) 

FORMAT (180,4X%;, &Be=",,F2D2097,2X, 143 -'A1PD20.7) 

BORMAT (100, 4% ,¢Ce='71PD20.7,2%,"+3 ©, 12D20.7) 

WRITE (6,85) 

BORMAT (AHO AIXeUD =",10X,°1.0", 9X,' +> *,11X%,'0.0'//) 
CALL ZXPOWL (UOBJ,EPS,N,X,FMIN, ITMAX,WA,IER) 

WRITE (6,81) 

FORMAT (1H1,4X,'THE NORMALIZED ABCD PARAMETERS ARE:'/) 
WRITE (6,82) X(1) ,X (2) 

WRITE (6,83) X (3) ,X (&) 

WRITE (6, €4) X(5) , X (6) 

WRITE (6,85) 


CR ROK ROOK IO kK cake ak a fe ae ake ic aie ak fe ai ak ofa ake ak a 
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CALCULATE THE ABSOLUTE ABCD PARAMETERS 
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AN=DCNPLX (X (1) ,X (2)) 
BN=DCMPLX (X (3) ,X (4) ) 

CN=DCMPLX (X (5) ,X (6)) 

DF2=1./ (Z0%* (AN-BN*CN) ) 

DX=DSQRT ( (DREAL (DF2) +CDABS (DF2) ) /2.0) 
DY=0.5*DINAG (DF2) /DX 

DF=DCMPLX (DX, DY) 

AF=ZO*AN*DF 

BF=Z0*BN*DF 

CP=CN*DF 

WRITE (6,87) | 

FORMAT (4X,'THE ABSOLUTE ABCD PARANETERS ARE:'/) 
WRITE (6,82) AF 

WRITE (6,83) BF 

WRITE (6,84) CF 

WRITE (6,88) DF 

FORMAT (150, 45 /UDe=%, 1PD2047,2%, "+0 ", 1PD20.77/) 
CHECK=AF*DF-BF*CF 

WRITE (6,89) CHECK 

FORMAT (10X,'THE VALUE OF AD-BC =',1PD20.7,2X,*+J ¢,1PD20.7) 


CHEK i 2 95 eB a ae ie OK IK Ok DK Be OK aK aK afc ae fs aie is 2K OK OB XE IKK OK OK 
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CHECK THE OPTIMUM NORMALIZED A,B,C,D 
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WRITE (6,92) 

FORMAT (1H0,*THE CALCULATED IMPEDANCES ARE:'*//) 

WRITE (6,93) 

FORMAT (6X,'BIAS*, 15X,'1 HHZ IMPEDANCE DATA',19X,'MICROWAVE IMPEDAN 
1CEB DATA',17X,'DIz IN CALC IMPEDANCES! 


4 A 5X (INOGPAGH GUI AVRESHSTANGE* ,11X, *REACTANCE', 10x, 
Ze tRSSLOCANCE © pit, GURBACTANCEY AG ker KR, 16x," DIFX'/ 


221%, (OHMS) ', 14X,* (OHMS) *) 

DO 95 I=1,NN 

ZHC (I) = (AN¥*ZD (I) +#BN) / (CN¥ZD (I) +1. 0) 
DET (I) = (2M (I) -ZuC (I) ) 

WRITE (6,94) V(I) ,ZD(I) , ZNC (I) , DET (1) 
PORNAT (2X,F 10.3, 1P6D20.7) 

CONTINUE 

ERR=0.0 
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DO 300 I=1,NN 
DUNM= (AN*ZD (I) +BN) / (CN*ZD(I)+1.0 )-2M (I) 
ERR=ERRtCDABS (DUNN) **2 
CONTINUE 
WRITE (6,301) ERR 
FORMAT (1HO, "ERROP=!,D20.7) 
FORMAT (1H0,2(19X,F10.3)) 
FORMAT (8F10.4) 
DO 251 I=1,NPON 
PIN (I} =PIQ (I) 
CONTINUE 
DO “1222 *10=1F2 
READ (5,500) (RLF (K) ,K=1,NPOR) 
READ (5,500) (THETF (K) ,K=1,NPOW) 
DO 600 I=1,NPOW 
WRITE (6,601) RLF (I) , THETF (I) 
CONTINUE : 
DO 501 K=1,NPOW 
PIF (K)=0.001* (10. ** (PIN (K) /10.)) 
RFF (K) =RLF (K) /50. 
THEFF (K) =THETF (K) /10. 
GAMA (K) =10. ** (REF (K) /20.) 
THETAR (K) =THEFF (K) *PI/180. 
GCOS=GANNA (K) *DCOS (THETAR (K) ) 
GSIN=GAMNA (K) *DSIN (EHETAR (K) ) 
ZF (K) =DCMPLX (1.+#GCOS,GSiN) /DCMPLK (1.-GCOS,-GSIN) 
YD (K) = (CN¥ZP (K) -AN) / (BN~-ZF (K) ) 
CONSV=CDABS ( (AN+CN) + ( (BN4+1.) *¥D (K) )) 
VAC (K) =DSQRT (8. PIF (K) *CDABS (AN- (BN*CN) )) /CONSV 
CONTINUE 
TY (O's5O. 1) WRITE (6,229) 
TF(LOssEO? 2) ONRETE (6,228) 
WRITE (6,502) 
DO 505 K=1,NPOH 
WRITE (6,503) PLF (K) , VAC (K) , YD (K) 
CONTINUE 
SMALL SIGNAL ADMITTANCE 
NPOR=7 


, dO 220 I=1,NPOW 


“ 1222 


PIN(I)=-20. 

CONTINUE 

CONTINUE 

COKTINUE 

FORMAT (1H1,'SMALL SIGNAL ADNITTANCES) 
FORMAT (1H1,'LARGE SIGNAL ADMITTANCE‘) 
FORMAT (1HO,6X,4 (D15.7, 10X) ) 


Zeal 


FORNA?T*(1H0,// 10X,* INPUT POWERS, 13%, “AC-VOLTAGE*, 20%, 'GDt, 24x, ' BD") 
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END 
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FUNCTION UOBJ(X) 
IMPLICIT REAL*8 (A-H, 0-2) 
COMMON/FUNC/Z4 (20) ,2D(20),NN,NC 
COMPLEX*16 AN,BN,CN,DELZD,DELZH,2N, DUMMY, 2D 
DIMENSION X (6) 
AN=DCHMPLX (X (1), X(2)) 
BN=DCMPLX (X (3} ,X(4)) 
CN=DCHPLK (X (5) ,X (6) ) 
UOBJ=0.0 
DO.100 Z=1,NC 
DUMMY= (AN*ZD (Z) +BN) / (CN*ZD (I) +1.0) - 24 (I) 
COBJ=UOBJ+CDABS (DUMMY) **2 
100 CONTINUE 
RETURN 
END 
SUBROUTINE ANCALI(X,LANDG,F,SHL) 
IMPLICIT REAL*8 (A~H,O-Z) 
REAL*8 LAMDG 
COMPLEX*16 GH,GE,AC,BC,CC,GCAL,ONE, UU,A,B 728, 2D, DET 
DIMENSION X(6),SHL(20) 
COMNON/FUNC/ZH (20) ,ZD (20) ,NN,NC 
COMMON/ALPHA/A (3,3) -B(3,1) ,GN(20) ,GE(20) ,GCAL (20) , DET (20), 
*RL (20) , THETA (20) ,GAHMA (20) , THETAR (20) , HA (GO) , PI 
EXTERNAL UOBJ 
N=6 
EPS=1.D-10 
WRITE (6,55) (RL(Z) ,I=1, NC) 
WRITE (6,56) (THETA (I) , I=1,NC) 
55 FORMAT (1H0,*RL=', 13 (F8.2, 1X) ) 
56 FORMAT (1H0,'THETA=",13(F8.2 ,1X)) 
DO 5 I=1,NC 
RL (I) =RL (I) /50.0 
THETA (I) =THETA (I) /10.0 
GAMMA (I) =10.** (RL (I) /20.) 
THETAR (I) =THETA (1) *PI/180.0 
GCOS=GANNA (1) *DCOS (THETAR (I) ) 
GSIN=GAUMA (I) *DSIN (THETAR (I) ) 
GM (I) =DCHPLX (GCOS,GSiN) 
5 CONTINUE 
DO 666 I=1,HNC 
COSV=DCOS (4. 0XPI*SHL (1) /LAMDG) 
SINV=DSIN (4. O¥PI*SHL (I) /LAMDG) 
GE (I) =DCNPLX (-COSV, SINV) 
666 CONTINUE 
: DCm2 0m les 
ie 
h (igs) =GE (J) 
A (I, 2) =DCMPLX (1.D0,0. D0) 
A (1,3) =-GH (J) *GE (J) 
B(1, 1) =Gi(J) 
20 CONTINUE 


{ ’ ; 
t2h2\h, iit be . 


THO .ca\ ee, ot: set oe | 


Or Us FOSS 
(OR) Sad, (OS; 2009, OS) Th, ORNS intl ae 
cal ae Lae ge et By. Ag FS) nt ae 


Sei as ries? a. 
BS 8 ie owe WS tee, so.) TOR 5 AAS. Fao! : oY ' 


a* 5 ates 
, F ee Te eee a Se 
f ae 
mf 
+ af #5 


b. been ae ty 
{iT} AALZKNS 


a ES aaa 


‘pice nase 
(mim ads (Tp enna 
af ao Sy aie) Lew 


ee 
~ roe . Z 
NiO Se. a 
+ Ae ° 
awe ees * ie eds 
i cS 
; we 9 
a - 
wnt id = AGM : Di 
: . 
‘ 
s % 4 
= ES SESE SAR PRM ART a ts Paci Ae oo ee a crete 
* u a 


ae ag h 4 ane Pa ath Sg cas 


as eae as 4 
so Duy ae a é van ee 


ie ade Lites 


CALL LEQTIC(A,3,3,B,1,3,0,WA, IER) 
WRITE (6,51) IER 
51 FORMAT (1HO,*AN-LEQTIC IER=',15) 
X (1) =DREAL(B(1,1)) 
X (2) =DIHAG(B(1,1)) 
X (3) =DREAL (B (2, 1)) 
X (4) =DINAG (B(2,1)) 
X (5) =DREAL (B (3, 1)) 
X (6) =DIHAG (B (3, 1)) 
WRITE (6,74) 


cag 


74 PORNAT (///, THE STARTING ABCD PARAMETERS FOR ANALIZER-CALI ARE:'/) 


WRITE (6,82) X(1) ,X (2) 
WRITE (6,83) X (3) ,X (4) 
WRITE (6,84) X (5) ,X (6) 
DO 220 I=1,KC 

ZD (1) =GE (I) 

2M (I) =GH (I) 

220 CONTINUE 4 
CALL ZXPOWL (UOBUJ,EPS,N,X,FMIN,20000,WA, IER) 
WRITE (6,52) LER 

52 FORMAT (1HO,*AN-ZXPOWL IER=",125) 
WRITE (6,81) 


81 FORMAT (1H1,4X,*THE OPTIMIZED ABC PARAMETERS FOR ANALYZER-CALI 


* RRE:'/) 
WRITE (6,82) ¥(1) ,X (2) 
WRITE (6,83) X (3) ,X (4) 
WRITE (6,84) X(5) ,X (6) 
CK KER Fe eR Oe He a OS OK Ne OK aK ie ts He oe Oe ae eK OK Ke ak ok 


S CHECK THE OPTIMUM NORMALIZED A,B,C FOR ANALYZER CALIBRATION 


CK Me WK He KK Ac eK RK HH Be SE HE 2 RK Oe RA HK At OK EK OK eK ORE OK 
ERR=0.0 
AC=DCHPLK (X (1) ,X(2)) 
BC=DCMPLX (X (3) ,X(4)) 
CC=DCMPLX (X (5) ,X(6)) 
ONE=DCHPLX (1-D0,0.D0) 
DO 600 I=1,NC 
GCAL (I) =(GE(Z) *AC+BC) / (GE (I) *CC+ONE) 
UOG=GH (1) ~GCAL (1) 
DET (Ij =UG 
ERR=ERRtCDABS (UU) **2 
600 CONTINUE 
WRITE (6,945) ERR 
945 FORMAT (1HO,'ERROP-AN=',D20.7) 
82 FORMAT (1HO,4X, AN=", 1PD20.7,2X,'tJ ',1PD20.7) 
83 FORMAT (1H0,4X,"BN=",1PD20.7,2X,'+d *,1PD20.7) 
84 FORMAT (1HO,4X,"CN=', 1PD20.7,2X,'+J ',1PD20.7) 
RETURN 
END 
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It is achieved by the following five programs. 


First: Evaluates the bicubic-spline coefficients 
g(i,j) or s(i,j) knowing conductance and 
susceptarce data of Inpatt diode. 

(Here FR = w and V = A) 
INPLICIT REAL*8 (A~H,0-Z) 
DIMENSION FK (5) ,V (6) ,FR(6),GR (6,6) ,C (4,4) , WK (340) , V1 (6) 
DATA W/0.0D0, .02Z5D0,. tD0 7.22500, -4D0;7. 625007 


DATA V1/ ~025D0,.1D0,.225D0,.4D0,.625D0/ 
DATA FR/.3249D0,.3481D0,. 3844D0,.4225D0/ 
DATA FK/ 3487D0,.3844D0,.4225D0/ 


DATA GR/12.441881D-3, 12. 731653D-3, 13.0336 75D-3, 13.347953D-3, 
*13.674482D-3, 14.013265D-3,12.625276D-3, 12. 874356D-3, 13.270122D-3, 
*13.7316540-3,14.178031D-3, 14. 52833 1D-3,13. 186755D-3, 13.458182D-3, 
*13.4984 15D-3, 13.73418D-3, 13.998617D-3, 14. 291723D-3, 14.613502D-3, 
*14.963952D-3/ 
SPLINE COEFFICIENTS FOR DIODE SUSCEPTANCE B(A2,82)/F 
SUN=0.0D0 
pO 2 K=1;5 
DOLL=A1 73 : 
CALL TGS2CU (GR, V, FR, 6,4, nlp C, WK LER) 
WRITE (6,10) K,L, IER 
10° FORMAT (1817, //5%, (R= 0/74 5G, (LH, 24 ASP BRE Th) 
WROTE (6 ,1.1)C (C (Dd a1 4). Et, 4) 
WROTE (12:1) C(G Gl ALY ai eg, T= 194) 
11 FORMAT(' ',////4D18.9,5X) 
D=V1 (K) -V¥ (K) 
B=FK (L) -FR (L) 
SUR=G (1,11). # (C (072) EB) (CV 7S) (BP 2y re (E0177, i (B473)i) # (C (272) FD) + 
1(C (3, 1) * (D**2)) + (C (4, 1) * (D¥*3) ) 
Do "9 u=2,4 
DO 7 N=2,4 
SUM =SOM +#C (M,N) * (D#* (4-1) ) * (B¥* (N=1)) 
7 CONTINUE 
9 CONTINUE 
WRITE (6,20) V1(K) ,SUM 

20 FORMAT (2(5X,D18.11)) 
1 CONTINUE 
2 CONTINUE 
1 FORMAT (4 (2X,D018.11)) 

STOP 

END 
RSS Po 10s Se eRe Os A) SER Th) os VE ES EM) SUPA ALI OD) OT IR. 
sy 


Subroutine ICS2CU was used from the International 
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SECOND: 
Evalutes the d-coefficierts of the diode krowirg 
the bicubic spline coefficients. 


THIS PROGRAM CALCULATES THE D-COEFFICIENTS OF GD(A2,W2) ,0 
BD (A2, 12) 

IMPLICIT REAL*8 (A-H,0-2Z) 

DIMENSION C(4,4),B(4,4),CN (4,4) ,A1(5,5,4,4) ,V (6) ,PR(6) 
DATE V0. 000, -02500,. 1D0,.225D0,.4D0,. 625D0/ 

DATA FR/.3249D-2, .3481D-2, .3844D-2, .4225D-2/ 
SUN=0.0D0 

DO '8 ‘K=1,5 

DO-e 1=1,3 

Remy Co, 21) 04 1 1(K, lu, 2,0) de 124), = 174) 
FORMAT (4 (2X,D18.11)) 


Sy Miegis eek 
CONTINUE 
B(1,1)=C (1,1) -C (1,2) *FR(L) +C (1,3) * (FR (L) **2) -C (1,4) * (FR (L) ¥*3) 

TEC OK2 tyme (2,2) PRL) Cll 3) Sieh (lle eos ee) ee (n) 

*+ (C (3,1) -C (3,2) *FR(L) tC (3, 3) * (FR(L) ¥*2) -C (3,4) * (PR (L) 83) ) ¥V (K) #2 

*4 (-C (4, 1) #C (4,2) ¥FR(L) -C (4,3) * (FR (L) ¥*2) #C (4, 4) ® (FR (LL) #3) ) *V (K) 2% 

*3 
B(1,2) =C (1,2) -2.*C (1,3) ¥FR(L) #3.*C (1,4) * (PR (L) ¥%2) 

K+ (=C (2,2) +2. *C (2,3) *FR(L) -3.*C (2,4) * (PR(L) *®2) ) *V (X) 

Kt (C (3,2) -2.*C (3,3) ¥PR(L) +3. *C (3,4) * (FR (L} *©*2) ) *V (K) #2 

Kt (-C (4,2) +2. *C (4, 3) XFR(L) -3.*C (4,4) * (FR(L) ¥*2) ) *V (K) ®¥3 
B (1,3) =C (1,3) -3.*C (1,4) *FR(L) + (-C (2, 3) +3. *C (2, 4) ¥FR(L) ) ¥¥ (RK) 

*4 (C (3,3) ~3.*C (3,4) *FR (L) ) ¥V (K) ¥*2=(C (4,3) -3.¥*C (4,4) ¥FR(L)) ¥V (K) ¥*3 
Big 4) =—C(1 7.4) =C (274) *V (K) #C (3,4) * (V (K) 42) -C (4,4) ¥ (V (RK) ¥*3) 
B(2,1)=C (2,1) -C (2, 2) *FR(L) +C (2,3) * (FR (L) **2) -C (2,4) * (FR (L) **3) 

*= 2% (C (3, 1) -C (3,2) *FR (L) #C (3,3) * (FR (L) **2) -C (3,4) * (FR (L) ¥*3) ) ¥V (K) 

*43.%(C (4,1) -C (4,2) *PR (L) #C (4,3) * (FR(L) **2) ~C (4,4) * (FR (L) **3) ) *V (K) 
Tx 2 
B (2,2) =C (2,2) -2.*C (2, 3) ¥PR(L) +3.*C (2,4) * (FR (L) **2) 

*4+2.% (-C (3,2) +2. *C (3,3) ¥FR(L) -3. *C (3,4) * (FR (L) ®*2) ) ¥V (K) 

*t3.%(C (4,2) -2.*C (4,3) ¥FR(L) +3.*C (4,4) * (FR (L) **2) ) ¥V (K) #*2 
B (2, 3) =C (2,3) -3-*FR (L) *C (2,4) +2.* (-C (3,3) #3. *C (3,4) *FR (L) ) *V (K) +3. 

1% (C (4,3) -3.¥*C (4,4) *FR(L) ) *V (K) **2 
B (2,4) =C (2,4) -2-#C (3,4) *¥ (K) #3. *C (4,4) * (V (K) ®*2) 

B (3, 1) =C (3,1) -C (3, 2) *FR (L) #C (3,3) ¥* (FR (L) *¥*2) -C (3,4) * (FR (L) **3) 
*43.%* (-C (4,1) +C (4,2) ¥FR(L) -C (4,3) * (FR (L) **2) +C (4,4) * (FR (L) **3) ) ¥V (K 
ok 

B (3,2) =C (3, 2) -2.*C (3, 3) *FR(L) +3.*C (3,4) * (FR(L) **2) 

#43. (-C (4,2) +2.¥*C (4,3) ¥FR(L) -3.*C (4,4) * (FR (L) ¥*2) ) ¥V (K) 
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B (3,3) =C (3,3) -3.*C (3,4) *FR (L) #3. * (-C (4,3) +3. *C (4,4) *FR (L) ) ¥V (K) 

B (3,4) =C (3,4) -3.*C (4,4) *V(K) 

B(4,1)=C (4,1) -C (4,2) *FR(L) +C (4,3) * (FR (L) ¥*2) -C (4,4) * (FR (L) **3) 
B(4,2)=C (4,2) -2.*C (4,3) *FR(L) +3. *C (4,4) * (FR (L) ¥*2) 

B (4,3) =C (4, 3) -3.*C (4,4) *FR(L) 

B(4,4)=C (4,4) 

DO 30 I=1,4 

DO 30 J=1,4 

FM= (FR (L+1) ** (J-1)) 

AM=(V (K+1) ** (I=1) ) 

SOM=SON+ (B (1,3) *FM*AM) 

CONTINUE 

FORMAT (1H ,5X,2D18.11) 

POGHALP (OTE ////4- DIG. 41759) 

WRITE (6,50) 

FORMAT (IY, “CN-COFFP!, 10%, *SUSCEPTANCE AS & EN OF W2, FOR TRANSIENT 
* CALCU. *) 
PRITECGpaa) Se (1,0), J=1,4),121,4) 

WREUE (7721) GB kid) 4J3=0, Why tad 

WRITE (6,60) SUM 

SUM=0.0D0 

CONTINOE 

CONTINGE 

STOP 

END 

THIRD: 
Evaluates the passive circuit cubic-spline coefficients 
as well as the corresponding d-coefficierts (Here x=w) 


IMPLICIT REAL*8 (A-H,O-2Z) 
DIMENSION GR(6),X(6),C (3,5) ,WK (70) , BK (5,4) 
DATA X/.390625D0, .4489D0, ~5476D0, -6084D0, -6724D0/ 
DOn2 yadda, 2 

READ (5,40) (GR(M) ,M=1,5) 

FORNAT (6 (D16.9,4X) ) 

WRITE (6,41) (GR(M) ,M=1,5) 

FORMAT (1H1,6 (4X,D16.9) } - 

CAuI Mees ICO(GR, X,5 ,C, WK, IER) 

WRITE (6,10) IER 

FORMAT (1HO, *IER=',110) 

wRIgEMe, 20) ((C (Ll, J) ,d=1,5 ),1=1,3) 
FORMAT(1H , 5(2X,D16.9)) 

DO eae d=1;4 

FA=X¢2) 

F2=P 1**2 

F3=F 1**3 

F4=X (I+1) 

F5=P 4% 2 

F6=F4**3 
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BK (I, 1) =GR(I) -C (1,1) ¥F14C (2,1) *F2-C (3,1) *F3 
BK (I,2)=C(1,1)-2.*C (2,1) *F14+3.*C (3,1) *F2 
BKC. 3) =C(2, 1) -3.*C (3, 1) ¥F 1 
BRAD Ly =64351) 
SUN=BK (I, 1) +BK (1,2) *F4+BK (I, 3) *P54BK (1,4) *F6 
WRITE (6, 30) SUM 

30 FORMAT (1H ,D16.9) 

1 CONTINUE 
WRITE(G6,15) ((BK (K,L) -L=1,4) ,K=1,4) 

15 FORMAT (1H0,4 (2X,D18.11)) 
WRITE (7,21) ( (BK (KL) ,L=1,4) ,K=1,4) 

2 CONTINGE 


21 FORMAT (4 (2X,D18.11)) 
STOP 
END 

FOURTH: 


Evaluates the quasistationary frequency w 
Having the active ard the passive d-coefficierts 
the value of w is determined by solving sixth 
order ron lizear algebriac equatiozr (Eq. 7.20). 
IMPLICIT REAL*8 (A-H,O0-2 
DEMENSZON. BK (4,4) GN (5,374,4) ,CN (5, 374,4) , GK (4,4) ,XCOP@ ),COFUT ) 
*,RR1(6 ),RIUI(6 ) 
DIMENSION GD (6) ,GT (6) 
DO 6 K=1,5 
DO 6 L=1,3 
GR EMO Ss oMG(GN (Kb, 0,0) d= 17,4) o= 174) 
DO 8 K=1,5 
DO 8 L=1,3 
SURED Ass )eChCN (kK, L, 1,0) ,0=1, 4) -2=1,4) 
READ (5 ,33)66(GKh(O, J)ezd=1 , def l= ta) 
READ (5,33) ( (BK (1,d) ,J=1,4) ,I=1,3) 
33 FORMAT (4(2X,D18.11)) 
D13=1./3. 
Di5=1./5. 
D17=1./7. 
A=0.0D0 
A2=A**2 
A4G=AX*Q 
A6=A**6 
K=1 
DO 994 HN=1,3 
DO 995 L=1,2 
CDO=GN (K,L, 1,1) +GN (K,L,2,1) *A2+#GN (K,L,3, 1) ¥A4+GN (K,L,4, 1) *46 
CD1=CN (K,L, 1,1) +3.*CN(K,L,2, 1) *A2*+5. *CN(K,L,3, 1) ¥AQ+7.*CN (K,L,4, 1) 
1*26 
CD2=-GN (K,L, 1,2) -GN(K,L,2, 2) ¥A2-GN (K,L, 3,2) *A4-GN (K,L,4,2) *A6 
CD3=-CN (K,L, 1,2) ~3.*CN (K,L,2,2) ¥A2~-5.¥*CN (K,L, 3,2) *AG-7.*CN(K,L,4, 2 
1) *A6 
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CD4=GN (K,L,1,3) +GN(K,L,2,3) ¥A2+GN (K, J, 3, 3) XAU+GN (K,L,4, 3) ¥A6 
CD5=CN (K,L, 1,3) +3. *CN (K,L,2,3) *A2¢5. *CN (K,L,3,3) ¥A4¢7.*CN (K,L,4, 3) 
1*A6 

CD6=-GN (K,L,1,4)-GN(K,L,2,4) *A2-GN(K,L,3, 4%) ¥A4-GN(K,L,4,4) ¥A6 
CD7=-CN (K,L, 1,4) -3.¥*CN (K,L,2,4) ¥A2-5.*CN (K,L, 3,4) ¥A4-7.¥CN (K,L,4,4 
1) ¥A6 

GOU2=CN(KeL, 171) + CNR, 2, 1) FAZFCN (K,L,3, 1) SA44CN (K,L,4, 1) *A6 
CD22=-GN(K,L,1,2) ~D13*GN(K,L,2,2) *22-D1I5*GN (K,L,3,2) *A4= 
1D17#GN (K,L,4,2) *A6 

CD32=-CN (K,L,1,2) -CN(K,L,2,2) *22-CN (K,L, 3,2) *A4-CN(K,L,4,2) ¥A6 
CD42=GN (K,L, 1,3) ¢D13*GN (K,L,2,3) £A2+D15¥GN (K,L,3,3) *A4+ 
1D17*GN (K,L,4, 3) *A6 

CD52=CN (K,L, 1,3) +CN (K,L,2,3) *A2+CN (K,L,3, 3) ¥A4*CN(K,L,4, 3} £6 
CD62=-GN (K,L,1,4) -D13*GN(K,L,2,4) ¥A2-D15*GN (K,1,3,4) *A4= 
1ID17*GN (K,L, 4,4) ¥46 

CD72=-CN(K,L, 1,4) -CN(K,L,2,4) ¥A2-CN (K,L,3,4) ¥A4-CN(K,L, 4,4) ¥A6 
WRITE (6,30) A,K,L, i 

30 FORMAT (1H ,*A=',D18.11,325) 

BO=BK (H, 1) #CD12 

Bi=0.0D0 

B2=BK (5, 2) -CD32 

B3=0.0D0 

B4=BK (4,3) +CD52 

B5=0.0D0 

BO=BK (Ni, 4) -CD72 

XCOF (1) =30 

XCOF (2) =B1 

XCOF (3) =B2 

XCOF (4) =B3 

XCOF(5 )=B4 

XCOF (6 )=B5 

XCOF(7 )=B6 + 

CALL POLRT(XCOF,COF,6 ,RR1,RIN1,1ER) 

WRITE (6,40) TER 

WRITE (6,50) 

WRITE (6,51) RR1 

WRITE (6,52) 

WRITE (6,51) RIN 

DO 200 KK1=1,6 

H1=RR1 (KK1) 

W2=W1**2 

W3=1**3 

W4=H1¥*4 

H5=W1**S 

WO=W1**6 

GD (KK1) =CD0-CD2*82+CD4*#4-CD6*WE 

GT (KK1) =CDO<GK (4, 1) + (GK (M1, 2) -CD2) *W2+ (GK (4,3) +CD4) *W4 
*+ (GK (NM, 4) -CD6) *W6 

200 CONTINUE 


T Subroutine POLRT was used from the Scientific Subroutine 


Package library (SSPLIB). 
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WRITE (6,800) 


WRITE (6,51) (GD (I) ,1=1,6) 

WRITE (6,801) , 
WRITE (6,51) (GT (1) ,I=1;6) 
CONTINUE 

CONTINUE 


FORMAT (1HO,*IER=",110) 

FORMAT (1H ,'REAL ROOTS') 

FORMAT (1H ,5(2X,D23.16)) 

FORMAT (1H ,*IM ROOTS*) 

FORMAT (1HO,*DIODE CONDUCTANCE!) 
FORMAT (1HO, ‘TOTAL CONDUCTANCE') 
STOP 

END 


FITTAs 
Havinrg w and d-coefficierts, this prograne evaluates 
the series state frequency vit x, assuming A, by 
solving eleventh order non linear algebriac equatiors 
(fqas3. 1.23) 
THIS ,PROGRAH EVALUATES THE ROOTS OF DA/DT=0 © DX/DT=0 AS FN OF 
SHERBOs), 2N 10 TH ORDER. IN, THE Two BONS.) X HAS TO BE THE SAE. 
EACH FN SET IN THE FORM: BOt+B1*X4+B2*X¥*24 53% X¥*3¢...4B1 14 X**11=0 
LUPLICIT REAL*S (Asi ,C=2) 
DIMENSIONS BK (4,4), GN lop 374 pdr Gl (043 pind), GK 4-4), XCOFP (12),,COF(12} 
*,RR1 (411) -RIM1(11) ,RR2(11) ,RIK2(11) 
#1=. 58864943 35896001D=1 
W2=H 1**2 
W3=W1**3 
WH=W1**4 
W5=H1*%*5 
W6=H 1**6 
W7=H1U**7 
D13=1./3. 
DAdS=1. 7/5. 
Di7=1./7. 
DG 6 K=1,5 
/DO 1 L=1,3 
PEADMOZ Sa) UGH (hy Loalb aii Cal ol ete eh 
CONTINUE 
CONTINUE 
DO 8 K=1,5 
DO, 2.L=1,3 
READ (5,33) ((CN(K,L,1,J9) ,J=1,4) ,I1=1,4) 
CONTINUE 
CONTINOE 
READ (5,33) ( (GK (1,J) ,J=1,4) ,i=1,3) 
READ (5,33) ((BK (I1-J) ,d=1,4) , I=1,3) 
FORHAT (4 (2X,D18.11) ) 
L=1 
i=2 
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READ (5,111) AA, DA,N 
WETETE (G, 012) AX, DAN 

TPVECENET (202735: 16 719) 

P12 GFORBAT (CtHO; AA=*7 D232. 16 ,o% fe DAK" /D 29. 16. N=", 15) 
DOCH 1 eC 1=t 8 
AA=AA+DA 
A=RA/ (10. DO*DSQRT (10. DO) ) 
iF (A piGire Us 0D eANDs 
IF (A ~Gke OS1S8774D0 Ir SAO 
IF (A sGue OF 31622500 SOAND. 
IF (A GE. 0.474342D0 .AND. 
IF (A e-GE. 0.632456D0 .AND. 

100 K=1 
GO TO 24 

200 K=2 
GO TO 24 

300 K=3 
GO TO 24 

301 K=4 
GO TO 24 

302 K=5 

24 A2=A**2 

AY=A%*Y 
AO=A¥*E : : 
CDO=GN (K,L, 1,1) +GN(K,L,2, 1) ¥A2#GN (K,L,3, 1) ¥A4#GN (K,L,4,1) *A6 
CDI=CN (K,L, 1,1) +3. *CN (K,L,2, 1) *A245.*CN(K,L,3, 1) ¥A4¢7.4CN (Ky L,4,1) 
T*A6 
CD2=-GN (K,L, 1,2) ~GN(K,L,2,2) *A2-GN(K,L,3, 2) ¥A4-GN(K,L,4, 2) *A6 
CD3=-CN (K,L, 1,2) -3.*CN(K,L,2, 2) *A2-5.¥CN (K,L, 3,2) *A4-7.¥*CN (K,L,4,2 
1) *A6 
CD4=GN (K,L, 1,3) +GN(K,L, 2,3) *A2tGN (K,L, 3, 3) *A44+GN (K,L, 4, 3) ¥A6 
CD5=CN (K,L, 1,3) +3.*CN (K,L, 2,3) #2245. ¥*CN (K,L,3, 3) ¥A4G#7 2 *CN (K,L,4, 3) 
1*A6 
CD6=-GK (K,L, 1,4) -GN(K,L,2,4) *A2-GN (K,L,3,4) *A4-GN(K,L,4,4) *A6 
CD7=-CN (K,L, 1,4) -3.*CN(K,L, 2,4) *A2-5.¥*CN (K,L,3,4) ¥A4-7.¥*CN (K,L,4,4 
1) *A6 
CD 122=CH (KeL: 171) tCN 1K, 1,2, 1) AZ CN (Ke, TPA ECON (Kh, Lb, 4, 1)6A6 
CD22=-GN (K,L, 1,2) ~DI3*GN (K,L,2,2) *A2-DI5S*GN (K,L,3,2) *A4- 
1D17*GN(K,L,4,2) *A6 
CD32=-CN (K,L, 1,2) -CN (K,L,2,2) *A2-CN (K,L,3,2) *A4-CN(K,L,4, 2) *A6 
CD42=GN (K,L, 1,3) +D13*GN (K,L,2,3) *A2+DIS*GN (K,L,3,3) *A4+t 
1D17¥*GN (K,L,4, 3) *A6 
CD52=CH (KGL, 1,3) FCN (KZ, 2,3) Ret CNA GS, oy) AatTCN( Kh, Lb, 4,3) "AG 
CD62=-GN(K,L,1,4) -D13*GN(K,L,2,4) *¥A2-D15¥*GN (K,L,3,4) *A4- 
1D17*GN (K,L,4,4) *A6 


a 


~LT. .158114D0) GO TO 100 
-LT. .316228D0 ) GO TO 200 
.LT. .474342D0) GO TO 300 
.LT. .632456D0) GO TO 301 
-LE. .790569D0) GO 70 302 
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CD72=-CN (K,L, 1,4) -CN(K,L,2,4) ¥A2-CN(K,L,3 74) KAK~CN(K, Lp 4,4) *26 
R11=CD1-CD3#3. *N2+CD5*5. ¥H4—-CD7*7. #HO*BK (H, 1) +BK (M,2) *3.*W2+ 

1 BK (H,3)*5.*84¢BK (M,4) ¥7. *R6 

R12=-CD3*6.*W1+CD5*20. *W3-CD7¥*42. ¥HS*BK (MH, 2) *6.*W1+BK (M,3) *20. *W3 
1#BK (4,4) ¥42.*H85 

R13= (BK ({1,2) -CD3) *3. DO+ (BK (M, 3) +CD5) ¥30. DO*W2+4 (BK (M,4) -CD7) ¥105. DO 
TeW4 

R14= (BK (4,3) CDS) *20. DO*¥W14 (BK (4,4) -CD7) *140. DO*N3 

R15= (BK (41,3) +CD5) *5.D0+ (BK (M,4) -CD7) *105. DO*W2 

R16= (BK (1,4) -CD7) *42. D0*W1 

R17= (BK (M,4) -CD7) *7. D0 

R21=GK (i, 1) *CD0+ (GK (1,2) -CD2) *W2+ (GK (M, 3) +CD4) *W44 (GK (11,4) -CD6) *16 
R22= (GK (M,2) -CD2) *2. DOFR14 (GK (1,3) #CD4) *4, DOXW3+ (GK (M, 4) -CL6) *6. DO 
1*u5 

R23=GK (N, 2) -CD2+4 (GK (M, 3) #CD4) #6. DO*W2+ (GK (MH, 4) -CD6) *15.D0*H4 

R24= (GK (i, 3) +CD4) #4. DOX R14 (GK (11,4) -CD6) #20. D0*H3 

R25=GK (ii, 3) #CD4+ (GE (4,4) -CD6) *15. DO*W2 

R26= (GK (fi, 4) -CD6) #6.D0*81 

R27=GK (N, 4) -CD6 

R3 1= (BK (M, 2) -CD3) *3. DO*# 1+ (BK (NM, 3) +CD5) *10. DO*W3+4 (BK (M,4) -CD7) *21. 
*DO*xRS 

R32= (BK (1, 2) -CD3) *3. LO+ (BK (4,3) #CD5) #30. DO*R2+ (BK (HM, 4) -CD7) *105. D0 
Twa 

R33= (BK (1,3) +CD5) ¥30. DOXH14 (BK (M, 4) -CD7) *210. DO*H3 

R34= (BK (M,3) +CD5) #10. DO+ (BK (M,4) -CD7) #210. DO*W2 

R35= (BK (M,4) -CD7} *105. DO*W1 

R36= (BK (M, 4) -CD7) #21. 50 

R41= (CD22-GK (M, 2) ) *2.DO*W1- (CD42+GK (4,3) ) *4. DOXN3+ (CD62-GK (MN, 4))* 
16.D0*W5 

R42= (CD22-GK (MH, 2) ) *2.DO- (CD42+#GK (H, 3) ) *12. DOXH2+ (CD62-GK (M, 4) ) *30. 
1DO*HY 

RY 3=- (CD&42+GK (M,3) ) *12.D0*W1+ (CD62-GK (M,4) ) *60.D0*H3 

R4G=- (CD42+GK (M,3))*4.D0 + (CD62-GK(M,4)) *60.D0*W2 

R45==- (CD42+GK (4,3) ) + (CD62-GK (11,4) } *30.D0*H1 

R4&6= (CD62-GK (M, 4) ) *6. DO*H1 

R51=CD12*W#1-CD32*W3+CDS2*W5-CD7I2* W7+BK (HM, 1) €W14BK (11, 2) *W3+ 

1BK (4,3) *W5+BK (4,4) ¥N7 

R52=CD12-CD32*3.*R2+CD52*5. *W4-CD72*7. *N64+BK (M, 1) +BK (HM, 2) *3.* N24 
1BK (4,3) *5.*H44+BK (M4) *7. *W6 

R53=-CD32*3. *¥H14CD52* 10. *83-CD72*21. *85+BK (M, 2) *3.*W14 

1BK (M,3) *10.*034BK (N, 4) #21. *H5 

R54=BK (N, 2) -CD32+ (BK (1,3) +CD52) *10.DO*H2+ (BK (M,4) -CD72) #35. DO*NA 
R55= (BK (4,3) +CD52) ¥5. DO*W14 (BK (HM, 4) CD72) *35. DO*W3 

R56=BK (i$, 3) +CD52+ (BK (M, 4) ~CD72) *21. DO*H2 

R57= (BK (UM, 4) CD72) 7. DO*H1 

R6 1=CD22-GK (HM, 2) — (CD42+GK (HK, 3) ) *¥6.DO*W2+ (CD62-GK (HM, 4) ) *15.DO*H4 
R62=~ (CD424+GK (M,3)) *12.D0*W 1+ (CD62-GK (4,4) ) *60.D0*W3 

R6 3=- (CD42+GK (4,3) ) *6.D0+ (CD62-GK (MN, 4) ) *90.DO*W2 

R64= (CD62-GK (4,4) ) *60.D0¥*H1 

R65= (CD62-GK (M, 4) ) *15.D0 

WRITE (6,30) AA,K,L,M 

30 FORMAT (1H ,'AA=',D23. 16,315) 
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BO=R31*R51-R51*R21 

B1=R31*R52+R32*R5 1-RG 1*R22~ RE2¥R21 
B2=R31*R53¢N3 2*R52+R33* RS 1- RE I*XR23— RSH 2*R22~ NS I*R21 

B3=R3 1¥R544R3 2*RS3+R33*RS2+4RIA*RS I= RH TFRQU=RH 2¥R2ZI—RGO 3B*R2Z2—ROY*R21 
B4=R31¥R554+R3 2*RSGtR33*R5 34+ R34¥RS 2+ R35S¥*RS 1- RO 1*R25- RG 2* RIG— RE 3¥ R23 
*= R6G*R22-N65*R21 

B5=R3 1*R56+R32*R55¢R3 3% RS44 RIG*R5D 3+ RI5S¥*R52+RIG*RS 1- RS 1*¥ R2G—- RO 2W*RIS 
#*— R63* R24 RHUXRZ3-RE5¥R22 

BO=R3 1¥ R574 23 2* R564 RS 3¥RH55¢RIGERS AE RZ5S¥RS 34 R36¥R52— RG 1*R27I-NGC 2*R2G 
*— R6 3% R25~R64*R24~ B65%*R23 

B7=R32*R5 7+ R33¥*R5G6¢R34*R554R35 *RS44+R36¥R5S 3— RG 2Q2*RIT—-RGOZB¥*R2G 

1- REG R25-RO5¥R24 

B8=R33*R5 7+ RBG R5G4R35*R55+ PB 36*%R5 4 RE 3¥R27— R64 *XR2G6— RE S¥*R2IH 
BO=R3K*RS 74+ R35* R56 R36¥R55—-ROY*R2 T= ROHS¥*RIG 
B10=R35*R57+R36*R56~-RO5*R2Q7 

B11=R36*R57 

XCOF (1) =BO 

XCOF (2) =B1 

XCOF (3) =B2 

XCOF (4) =B3 

XCOF(5 )=B4 

XCOF (6 )=B5 

XCOF (7 )=B6 

XCOF (8 )=B7 

XCOF(9 )=B8 

XCOF (10) =B9 

XCOF (11) =B10 

XCOF (12) =B11 

CALL POLRT (XCOF,COF,11,RR1,R1M1,IER) 

WRITE (6,40) LER 

WRITE (6,50) 

WRITE (6,51) RR1 

WRITE (6,52) 

WRITE (6,51) RIM1 

CO=R41*221-R11*R51 

C1=R41*¥R224+R42*R21-R11*R52—R12*R51 

C2=R41* R234 RU2* R2 2+ RY 3*¥R2Z1-R1 1*RS3—R I2*R5S2-RU3*RS1 

C3=R4 1*KR24+ BU 2% R23+ RU 3*R22+ RUG*ERZ 1- RIV IFRS 4- RY 2*R53-R1Z¥RS2-RIG*RS 1 
C4=RY1¥R254 RY2¥R2G4 RY 3¥ R234 R4G*RD 2¢ R4US*R21- R11 ¥R55~R12* RO4=R13*RS3 
*= R1G*R52-R15¥*R51 

C5=R4 1*R264 RA QERZ54 RY 3X R244 RAG*RIZ34 R45S*RI24ERGSH*R2I-R11* RS 6—RU2KRS5 
*—R13%R54-R1G*¥R5 3-R1S*R52-R1G*RS1 

C6=R41*R274RU2* R264 RA Z*R25+ RYY* R244 RUS*R2 34+ RYGXRQ2—R1U*¥RST—-R1U2* R56 
%— R1Z*R55-RIG*RSY-R1S*R5 3-R1G¥RS2- R1IT*RS1 

CT=R42*R27+ R4 3*RIGERAU*RIZS+RGS*RZU 4 RUG¥*R2IZ—RIQ*RS T-R13*¥R56—-RV4*R55 
*— R15*R54-R16¥R5 3-RIT*R52 
C8=R43¥RZ7+R44E*R264R45*¥R254 RUG¥R24— RUZ¥*R5S 7— RY4*R5S6-R1S*R55-R16*R54 
*~R17*R53 
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C9O=R44*R27+ RYS5S*R26+ RYG6*R25—RUG*RS 7-R15*R56—R16¥*R55-B17* RSM 


C10=R45*R27+R46*R26-RI5S¥*R57—RIG¥R56—-RIT*R55 
C11=R46*R27-R16*R57-R17* R56 ; 
XCOF (1) =C0 

XCOF(2)=C1 

XCOF (3) =C2 

XCOF (4 )=C3 

XCOF(5 )=C4 

XCOF (6 )=C5 

XCOF (7 )=C6 

XCOF (8 )=C7 

XCOP (9 )=C8 

XCOF (10) =C9 

XCOF (11) =C10 

XCOF (12) =C11 

CALL POLRT (XCOF,COF,11,RR2,RIN2, IER) 
WRITE (6,40) IER 
FORMAT (1HO, *IER=',110) 

WRITE (6,50) 

FORMAT(1H ,*REAL ROOTS") 

WRITE (6,51) RR2 

FORMAT (1H ,5(2X,D23.16)) 

WRITE (6,52) 

FORMAT (iH ,'IM ROOTS') 

WRITE (6,51) RIN2 

CONTINUE 

STOP 

END 
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Appendix III 


6 


OSCILLATOR TRANSIENT SOLUTIONS A(t) & x(t) 


It is obtained by solvirg two ron linear differential 
Sgustions is Aguex and wame(tje. 
The user must supply the following parameters. 
d-ccoefficierts of the diode 
d-coefficients of the passive circuit 
W qusistatiornary frequency 
The initial values of A and x 


IMPLICIT REAL*8 (A-H,O0-Z) 
DIHENSION X0(2),#K(15) 
CONNON /ONE/GN (5,3,4,4) ,CN(5,3,4,4) ,BK (4,4) ,GK (4,4) 
COMMON /TRO/G,H1 
EXTERNAL F 
PI=3.1415927D0 
H=-10.D0 
AA=-49500.D0 
XO (1) =.9624812145501785D-2/ (10. D0*DSORT(10.D0}) 
XO (2) =—. 334 1535974404524D-04 
W1=. 58864 94335896001D-1 
N=2 
DO 6 K=1,5 
DO 60L=1, 3 
6 herp (Syocid (CN (Kp hl, 3) ,0= t,4)),2=1 78) 
DO 80K=1,5 
HO38 ‘L=1,3 
6 BEA0G5 233) ((CN(K,L, 1,9) ,J=1,4) =1,4) 
REED (5, 33) C(GK (1, 3)ida1, 8) = 10} 
READ(5,33) ((BK (I,J) ,J=1,4) , I=1,3) 
33 FORMAT (4(2X,D18.71)) 
BB=-100.D0 
DO t10081=500, 120075 
B=BB*I 4 
CALL DASCRU(F,AA,B,H,N, X0,WK) 
AA=B 
AMPL=10.DO*DSQRT (10. DO) *X0 (1) 
FREQ=X0 (2) «1.011 
TIME=B/ (2.D0*PI*1.D11) 
HRITE (6,200) AUPL, FREQ, TIME,H,B 
100 CONTINUE 
200 FORMAT (1H , ‘AMPL=",D23. 16, 1X, *FREQ=',D23. 16,3X,'T=',D15.7, 1X, "H=", 
“043.7, 1%, (B=",D1567) 
STOP 
END 


t Subroutine DASCRU was used from the International 


Mathematical and Statistical library (IMSL). 
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200 
300 
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21 
400 
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22 
150 
250 
350 


SUBROUTINE F(X0,2,N,YP) 


JHPLICIT REAL*8 (h~H, 0-2) 
DIHENSION X0(2),Y¥2(2) 
COMMON /ONE/GN (5,3,4,4),CN(5,3,4,4) ,BK (4,4), GK (lt, th) 
COMHON /THO/G, #1 


X=K0K2) 
K2=KHK2 
X3=K%**3 
X4=X%x*Q 
X5=XeKS 
X6=X**6 
A=X0 (1) 
A2=A%*2 
AG=A%RG 
A6=A**6 


H2=W1**2 
W3=4 14%3 
W4=W 15 *4 
H5=8 1**5 
W6=W1x*6 
W7=W Tk *7 


D13=1./3. 
D15=1./5. 
D17=1./7. 


LF (A 
IF (A 
LF (A 
IF(A 
IF (A 
K=4 


GOSTO 124 


K=2 


GOsTOe2 1 


K=3 


GOg70 721 


K=4 


GO TO 21 


Ks 
IF (8 1#+X 
IF (W1+4K 
IF (W1+X 
L=1 


IF (W1+X 
IF (W1+X 
IF (H1+X 
M=7 


GO TO 23 


M=2 


GO IO 23 


eG. 
«GE. 
GE, 
«GE. 
GE. 


«GE. 
-Gk. 
GE. 


GE. 
GE. 
-GE. 


0.0D0 
0. 155%7148D0 
0.316228D0 
O.474342D9 
0.63245 


~C57 «DO 
- 0590D0 
- 0620D0 


057 «DO 
-0575D0 


2 0582550 


6D0 


~ AND. 
AND. 
AND. 
APS RODS 


~ AND. 


eA ID < 
e AND. 
e AND. 


be fo cet fe 


Wi+X 
W1+k 
W14+X 


W1+X 
W1+X 
WitX 


~ 


-LT. 
~LT. 
«LT. 
LT. 
LT. 


LT. 
-LT. 
e LE. 


-LT. 
-LT. 
~ LE. 


- 1581718D0) 


~21522800 } 


-474 34200) 
-632456D0) 
- 79056900) 


-059D0 )GO 
.062D0 )GO 
- 06500 


were) 


- 057550} GO 
.05825D0) Go 
. 059250036 


TO 
TO 
TO 


?O 
To 
LO 


9 


TO 
TO 
TO 
Hh) 


400 
509 
600 


150 
250 
359 


100 
206 
200 
301 
202 
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23 CDO=GN (K,L,1,1) +GN(K,L,2, 1) *A2*GN (K,L,3, 1) *A44+GN (K,L,4,1) *A6 
CDI=CN(K,L, 1,1) #3. *CN (K,L,2,1)*A245. *CN (Ky 1,3, 1) AUT. £CN (Ky Ly 4,1) 
1*A6 
CD2=-GN (K,L, 1,2) -GN(K,L,2, 2} *A2-GN (K,L,3, 2) *A4- GN (K, L,4,2) *A6 
Chs=-CN (K,1L, 1,2) -3. xCN(K,L, 2, 2) FA 2=55 *CN(K, L,3,2) *A4- We *CN(K,L,4,2 
1) *A6 
CD4=GN (K,L,1,3) +GN(K,L,2,3) *A2€GN (K,L, 3,3) *A44GN (K,L,4, 3) ¥A6 
CD5=CN (K,L, 1,3) #3. ¥CN (K,L, 2,3) ¥A2*5.*CN (K,L,3, 3) *AG+7.*CN (K,L,4, 3) 
1*2hO 
CD6=-GN(K,L, 1,4) -GN(K,L,2,4) *A2-GN(K,L,3, 4) ¥AU-GN(K,L, 4,4) *A6 
CD7=-CN (K,L, 1,4) -3.¥*CH (K,L, 2,4) *A2-5. #CN (Ky L, 3,4) *AG-7.*CH (K,L, 4,4 
1) ¥AS 
CD1I2Z=CN (K,L, 1,1) +CN (K,L,2, 1) ¥A2+CN (K,L,3, 1) ¥AG#CN (K,L,4, 1) ¥A6 
CD22=-GN (K,L, 1,2) ~D13#GN (K,L,2,2) #A2-D15¥GN (K,L, 3,2) *A4~— 
1D17¥GN (K,L,4,2) ¥6 
CD32=-CN (K,L, 1,2) -CN (K,L,2, 2) ¥A2-CN (K,L,3,2) ¥A4-CN(K,L,4,2) ¥A6 
CD42=GN (K,L, 1,3) +D13*GN (K,L,2,3) *A2+D15*GN (K,L,3,3) ¥A4+ 
1D17*#GN (K,L,4, 3) *A6 
CD52=CN (K,L,1,3 ) #CH (Kel 2,3) ¥A2+CN (K,L,3,3) ¥AG4CN (K,L,4,3) *A6 
CD62=-GN (K,L, 1,4) -D13¥*GN(K,L,2,4) *A2-D15*GN (K,L,3,4) ¥A4- 
IDI7=GN (K,L,4,4) *A6 
CD72=-CN (K,L, 1,4) -CN (K,L,2, 4) *A2-CN (K,L, 3,4) ¥AU-CN(K,L, 4,4) ¥A6 
R17 1=CD1-CD3*3 .*W24+CD5*5.*N4-CDI*7. *H6+BK (NM, 1) +BK (Kh, 2) *3.*82+ 
1 BK(M,3) *5.*W44BK (4,4) *7. *W6 
Ri2=-CD3*6.*W14+CD5*20. *W3-CD7*42. *WStBK (M, 2) #6. *WN1+BK (A, 3) *20.*W3 
Be 4) *42.*W5 
R13= (BK (MN, 2) -CD3) *3.DO+ (BK (M, 3) #CD5) *30. DOXH2Z+ (BK (MH, 4) -CD7) *105. D0 
1*H4 
R14= (BK (1,3) #CD5) ¥20. DOXH 14 (BK (M, 4) -CD7) *140. DOXW3 
R15= (BK (4,3) +CD5) *5.DO+ (BK (4,4) -CD7) #105. DO*N2 
R16= (BK (M, 4) -CD7) #42. DO*H1 
Ri7= (BK (4,4) -CD7) ¥7.D0 
R21=GK (MN, 1) +CDO+ (GK (M, 2) -CD2) ¥W2+ (GK (M, 3) +CD4) *H4+ (GK (M,4) -CD6) *R6 
R22= (GK (M, 2) -CD2) *2. DOXW1+ (GK (NM, 3) #CD4) #4. DO*W3+ (GK (M, 4) -CD6) *6. D6 


R23=GK (N, 2) -CD2+ (GK (MN, 3) #CD4) *6.DO*W2+ (GK (4,4) -CD6) *15. DO* HY 

R24= (GK (11, 3) +CD4) *4. DO*W1+4 (GK (M,4) -CD6) #20. DO*W3 

R25=GK (I, 3) #CD4+ (GK (H,4) -CD6) *15. DO*H2 

R26= (GK (MN, 4) -CD6) *6. DO*HT 

R27=GK (i, 4) -CD6 

R3 1= (BK (MH, 2) -CD3) *3.DO*W1+ (BK (M,3) +CD5) *10. DO*N3+ (BK (M,4) -CD7) *21. 
*DOxHS 

R32= (BK (1, 2) -CD3) *3. DO+ (BK (N,3) +CD5) *30.D0*W2+ (BK (M,4) -CD7) *105. D0 
1*84 

R33= (BK (1,3) #CD5) ¥30. DO*H14# (BK (IM, 4) -CD7) #210. DO*H 

R34= (BK (f1, 3) +CD5) #10. DO+ (BK (M,4) -CD7) #210. DO*W2 

R35= (BK (1,4) -CD7) *105. DO*wW1 
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R36= (BK (4,4) ~CD7) *21. D0 

R4 1= (CD22-GK (MN, 2) ) *2. DO*H1- (CD42+GK (MN, 3) ) *4.D0*W3+ (CD62~-GK (NM, 4} ) * 
16.D0*K5 ; 

RY 2= (CD22-GK (MN, 2) ) *2. DO- (CD42+GK (H, 3) ) *12. DO*H2+ (CD62-GK (MH, 4) ) *30. 
1DO*W4 

R43=— (CD42+GK (N,3)) *12.D0*W1+4+ (CD62-GK (M,4) ) ¥60.D0*W83 

R44=~ (CD42+GK (M,3)) *4.D0 + (CD62-GK (N,4) ) *60.D0*H2 

R45==— (CD42+GK (M,3) ) + (CD62-GK (M,4) ) *30.D0*H1 

R46= (CD62-GK (f, 4) ) *6. DO*H1 

R5 1=CD12*01-CD32*W34+CD52*N5S-CDI2*WH7+BK (M, 1) *HI*#BK (M, 2) *W3+ 

1BK (2, 3) #H54BK (4,4) *W7 

RO2=CD t2-CDS 273, .* R2eCDS ZFS. 2UN=-CD7T2*7 SHNGtBK(N, U) BK (1,2) 43. th Ze 
1IBK (4,3) *5.*W4eBK (4,4) *7.* 96 
R53=—-CD32*3.*K14+CD52* 10. *W3-CD72%21. *H5¢4+BK (M, 2) *3,*H1+ 

1BK (4,3) *10.*H3+BK (4,4) *21.*85 

R54=BK (M, 2) -CD32+ (BK (4, 3) #+CD52) *10.D0*H2+ (BK (4,4) -CD72) *35. DOFWY 
R55= (BK (M,3) +CD52) *5. DO¥H1+ (BK (HM, 4) -CD72) *35.D0*H3 

R56=BK (4,3) +CD52+ (BK (M, 4) -CD72) *21.D0*W2 

R57= (BK (4,4) -CD72) *7.D0*W1 

R6 1=CD22-GK (4,2) — (CD42¢GK (11,3) ) #6. DO*W2+ (CD62~-GK (KM, 4) ) *15.DOFRE 
R62=— (CD42+GK (NM, 3) ) *12.D0*W1+ (CD62-GK (N,4) ) *60.D0*83 

R63=— (CD424+GK (MW, 3) ) *6.D0+ (CD62-GK (NH, 4) )¥90.D0*H2 

R6 4= (CN62-GK (4, 4) ) *60.D0*#1 

R65= (CD62-GK (4,4) ) *715.D0 

R= RV14R12*X4+R1IZB*¥X2¢RI4*XZ34R1S*¥X44R16¥*X54+R17% X6 

R2=A* (R274 R22*X*tR2ZeX 2+ R24¥*X3+R25*X4+R26*X5S¢R2Z7I* XG) 

R3=A* (R31*+K32*X4+R33*X2ERIG*XZERIS¥*¥KN+tRZO*X5) 

R4= R4W4+R4E2*X+R43*X2¢4+ R4U*XZ4R4ES*XUtRAGKKS 

R5=A* (R514 R52*X+R53*X2+R5SU*XS4RS55*X44+R56*X54R57* XG) 

R6=A* (RO14+RO2*X+ROB*X2+RO4*XI+ROS* KU) 

YP (1) = (R3*R5-R6*R2) / (RO*R1I-RI¥RG) 

YP (2) = (R4*R2-Ri*RK5) / (RO*¥RI-R3*RY} 

RETURN 

END 
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